Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,752,107:1 +G 15.4% coding (116/975 nt) yiaU → putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,752,1071.G15.4% 585.5 / 116.4 175coding (116/975 nt)yiaUputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (68/80);  new base G (9/18);  total (77/98)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

TTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTA  >  NC_000913/3752008‑3752206
                                                                                                    |                                                                                                   
ttAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCaa                                                                                                      <  2:3583811/100‑1 (MQ=255)
   aGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ata                                                                                                  <  1:3948357/100‑1 (MQ=255)
    gTACCGGGAGTTAAAAATTATCTCGGTGATCGCTGCCAGTGCAAATAACAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCCACCTCCGCAA‑ctat                                                                                                 >  1:1084796/1‑100 (MQ=255)
      acagggAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ata                                                                                                  <  1:1424998/94‑1 (MQ=255)
      aCCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCt                                                                                               <  2:887850/100‑1 (MQ=255)
        cGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTg                                                                                             >  2:4054298/1‑100 (MQ=255)
           gAGTTAAAAATTATCTCGGTAATCTCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGc                                                                                            <  1:481080/98‑1 (MQ=255)
           gAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTg                                                                                          <  1:3372593/100‑1 (MQ=255)
            aGTTAAAAATTATCTCGGTAATCGCTGCGTGTGAAAATAACAGCCATGCCGCGACGGTACTTGGCATCGCAAAGGCCAACGTCAGCAA‑ATATCTTGc                                                                                            <  1:1853843/97‑1 (MQ=255)
             gTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATAT‑TTGCTGac                                                                                        <  2:4296305/100‑2 (MQ=255)
             gTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATAT‑TTGCTGac                                                                                        <  2:4316244/100‑2 (MQ=255)
             gTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGACCAACGTCAGCAA‑ATATCTTGCTGAt                                                                                        >  1:2184105/1‑100 (MQ=255)
               tAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATAT‑TTGCTGacgt                                                                                      <  2:4074396/100‑4 (MQ=255)
               tAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGAttt                                                                                      >  2:3660999/1‑100 (MQ=255)
                aaaaaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATAT‑TTGCTGacgtt                                                                                     <  1:1281205/100‑5 (MQ=255)
                aaaaaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGAtttt                                                                                     <  2:2690403/100‑1 (MQ=255)
                 aaatTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATAT‑TTGCTGacgttg                                                                                    <  1:4230385/100‑6 (MQ=255)
                 aaaaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTg                                                                                    <  1:2617046/100‑1 (MQ=255)
                 aaaaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTg                                                                                    <  2:2263658/100‑1 (MQ=255)
                 aaaaTTATCTCGGTAATCGCTGCCAGTGAAAATATAAGCAATGCCGCGACTGTACTTGGCATCGCACAGTCCAACGTCAGCAAGATAT‑TTGCTGacgttg                                                                                    <  1:3563355/100‑6 (MQ=255)
                  aaaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATAT‑TTGCTGacgttgg                                                                                   <  2:2334328/100‑7 (MQ=255)
                  aaaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATAT‑TTGCTGacgttgg                                                                                   <  2:3393473/100‑7 (MQ=255)
                   aaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATAT‑TTGCTGacgttggc                                                                                  <  1:3576582/100‑8 (MQ=38)
                   aaTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGaa                                                                                  <  2:1058507/100‑1 (MQ=255)
                    aTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTATTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAAt                                                                                 <  1:432462/100‑1 (MQ=255)
                    aTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATAT‑TTGCTGacgttggcc                                                                                 <  2:1857516/100‑9 (MQ=38)
                    aTTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAAt                                                                                 <  2:3754602/100‑1 (MQ=255)
                      tATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACATCAGCAAGATAT‑TTGCTGacgttggcctg                                                                               <  2:2126332/100‑11 (MQ=37)
                       aTCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCaa                                                                              <  1:170016/100‑1 (MQ=255)
                          tcGGTAATCGCTGCCAGTTAAAATATCAGCCATTCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAg                                                                           <  1:3445501/100‑1 (MQ=255)
                           cGGTAATCGCTGCCAGTTAAAATAACAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAg                                                                           >  1:3125963/1‑99 (MQ=255)
                           cGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGt                                                                          >  1:2900706/1‑100 (MQ=255)
                             ttAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTgg                                                                        <  1:2289759/99‑1 (MQ=255)
                               aaTCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGt                                                                      <  1:3180548/100‑1 (MQ=255)
                               aaTCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGt                                                                      <  1:3155145/100‑1 (MQ=255)
                               aaTCACAGCCAGTGAATATATCAGTCACGCCGCAACTGTACTTGGTTTCGGATAGGCCAACGTCAGCAA‑GTATCTTGCTGCTTTTGAATCAAACGCGGGt                                                                      <  1:2802295/100‑1 (MQ=255)
                                aTCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGACTCAAAAGTggggt                                                                     <  2:2983574/100‑3 (MQ=255)
                                   gCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGttt                                                                    <  1:4192133/98‑1 (MQ=255)
                                   gCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCGTGCTGATTTTGAATCAAAAGTGGGTTTaa                                                                  >  1:3458546/1‑100 (MQ=255)
                                    cTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTaaa                                                                 >  2:3083590/1‑100 (MQ=255)
                                    cTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTaaa                                                                 >  1:4099759/1‑100 (MQ=255)
                                    cTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTaaa                                                                 <  2:3311681/100‑1 (MQ=255)
                                     tGCCAGTGTAAATATTTGCTATGCTGCGATTGTACTTGGTATTGAAGAGGAAATCGTCTGCAA‑ATATCTTGCTGATTTTGTATCAAAAGTGGGTTTa                                                                   <  1:3083893/97‑1 (MQ=255)
                                      gCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAg                                                               >  2:4240713/1‑100 (MQ=255)
                                        cAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTc                                                             >  2:2329492/1‑100 (MQ=255)
                                         aGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTGTGAATCAAAAGTGGGTTTAAAAGTCt                                                            >  2:2187282/1‑100 (MQ=255)
                                          gTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCtt                                                           >  2:182026/1‑100 (MQ=255)
                                           tGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGTTTTTGAATCAAAAGTGGGTTTAAAAGTCtt                                                           <  1:2891888/99‑1 (MQ=255)
                                           tGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCttt                                                          <  2:2749662/100‑1 (MQ=255)
                                            gAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTg                                                         <  1:4091087/100‑1 (MQ=255)
                                            gAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTg                                                         <  2:155108/100‑1 (MQ=255)
                                            gAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTg                                                         >  1:680602/1‑100 (MQ=255)
                                            gAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTg                                                         <  2:666000/100‑1 (MQ=255)
                                               aaTATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGAcc                                                      <  2:3255597/100‑1 (MQ=255)
                                                 tatCAGCCATGCCGCGACTGTACTTGGCATCGCACATGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGa                                                    >  2:3099311/1‑100 (MQ=255)
                                                   tCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGaan                                                  <  2:3874436/100‑2 (MQ=255)
                                                    cAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAacc                                                 <  1:2230303/100‑1 (MQ=255)
                                                    cAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAacc                                                 <  2:4382157/100‑1 (MQ=255)
                                                     aGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAacca                                                >  2:857235/1‑100 (MQ=255)
                                                     aGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAacca                                                >  1:3003774/1‑100 (MQ=255)
                                                      gCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAaccac                                               >  2:83485/1‑100 (MQ=255)
                                                        cATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCAccc                                             >  1:3847171/1‑100 (MQ=255)
                                                        cATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCAccc                                             >  2:130136/1‑100 (MQ=255)
                                                             cgcgACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCcagcag                                        >  1:886386/1‑100 (MQ=255)
                                                              gcgACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGc                                       >  2:4291398/1‑100 (MQ=255)
                                                              gcgACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGc                                       <  2:3803992/100‑1 (MQ=255)
                                                              gcgACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGc                                       >  2:3736908/1‑100 (MQ=255)
                                                                   tGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATg                                  <  2:3907694/100‑1 (MQ=255)
                                                                    gTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGc                                 <  2:3342455/100‑1 (MQ=255)
                                                                     tACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCt                                >  2:1733057/1‑100 (MQ=255)
                                                                      aCTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTc                               <  1:1133508/100‑1 (MQ=255)
                                                                      aCTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTc                               >  1:241217/1‑100 (MQ=255)
                                                                      aCTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTc                               >  2:2352400/1‑100 (MQ=255)
                                                                      aCTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTc                               <  2:1234309/100‑1 (MQ=255)
                                                                       cTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTca                              >  1:637860/1‑100 (MQ=255)
                                                                         tGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTcaca                            >  2:1614948/1‑100 (MQ=255)
                                                                            cATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCt                         >  2:966928/1‑100 (MQ=255)
                                                                               cGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCtttt                      <  1:4373230/100‑1 (MQ=255)
                                                                                 cacaGGCCAACGTCAGCAC‑ATATCTTGCTGCTTTTGACTCACAAGTGGGTTTAACAGTCTTTGTCCGAACCACCCGGCAGCTGATGCTCACGCCTTTTgg                    <  2:4126587/100‑1 (MQ=255)
                                                                                 cacaGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACTGAACCACCCGGCAGCTGATGCTCACACCTTTTgg                    <  2:3925579/100‑1 (MQ=255)
                                                                                  acaGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGc                   <  2:2459944/100‑1 (MQ=255)
                                                                                  acaGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGc                   <  1:391216/100‑1 (MQ=255)
                                                                                    aGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCAc                 <  2:3637207/100‑1 (MQ=255)
                                                                                     attcaaaatCAGCAAGATAT‑TTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGCCCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCAc                 >  2:40294/9‑99 (MQ=37)
                                                                                      ttcaaaatCAGCAAGATAT‑TTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACcg               >  1:773050/8‑100 (MQ=38)
                                                                                       tcaaaatCAGCAAGATAT‑TTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACcgc              >  1:2468643/7‑100 (MQ=255)
                                                                                       ccAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACcgc              >  1:3143515/1‑100 (MQ=255)
                                                                                           aatCAGCAAGATAT‑TTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGACACCGCGCTg          >  1:4123628/3‑100 (MQ=255)
                                                                                             tCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGtt        <  2:1664544/100‑1 (MQ=255)
                                                                                              cAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTa       >  2:3579465/1‑100 (MQ=255)
                                                                                               aGCAA‑ATATCTTGCTGATTTTTAATCTAAAGTGTGTTTAAAAGTCTTTGATCGAACCAGCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTAc      <  1:3021115/100‑1 (MQ=255)
                                                                                               aGCAA‑ATATCTTGCTGATTTTAAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTAc      <  1:2769675/100‑1 (MQ=255)
                                                                                                gCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTAcc     >  2:328486/1‑100 (MQ=255)
                                                                                                gCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTAcc     >  1:2678708/1‑100 (MQ=255)
                                                                                                gCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTAcc     >  1:3193177/1‑100 (MQ=255)
                                                                                                gCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTAcc     >  2:295573/1‑100 (MQ=255)
                                                                                                 cAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCg    <  2:2052534/100‑1 (MQ=255)
                                                                                                  aa‑aTATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGt   >  1:3850399/1‑100 (MQ=255)
                                                                                                  aa‑aTATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGt   >  1:916876/1‑100 (MQ=255)
                                                                                                   a‑aTATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTa  <  1:4031184/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
TTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAA‑ATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACCGAACCACCCGGCAGCTGATGCTCACACCTTTTGGCACCGCGCTGTTACCGTA  >  NC_000913/3752008‑3752206

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑