Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,978 A→C 27.4% intergenic (‑188/+120) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9780AC27.4% 368.6 / 151.0 204intergenic (‑188/+120)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base A (61/75);  new base C (42/3);  total (126/78)
Fisher's exact test for biased strand distribution p-value = 2.12e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-03

Only 100 of 208 total aligned reads displayed.

GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705879‑3706077
                                                                                                   |                                                                                                           
   gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                          <  2:104716/100‑1 (MQ=255)
     aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTCTACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                        <  1:858909/100‑1 (MQ=255)
      ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                       <  2:1414241/100‑1 (MQ=255)
        gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                     <  2:1552556/100‑1 (MQ=255)
        gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                     >  2:627113/1‑100 (MQ=255)
         ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                    >  1:2711336/1‑100 (MQ=255)
          tggccgATGGGGCTGCACGCAGAGGCACAGCGCCAGCGGGTATGCAGATGGGTGACTAAGTGTGCTTTGTTTATGCCGGATAAGGGGTTAACGccgca                                                                                                     <  1:1204262/93‑4 (MQ=19)
          tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTCGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                    >  1:2530029/1‑99 (MQ=255)
              ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                               <  2:812152/100‑1 (MQ=255)
               gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                              <  2:2627334/100‑1 (MQ=255)
               gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                              <  2:2511255/100‑1 (MQ=255)
                aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                             <  2:2556677/100‑3 (MQ=37)
                aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                             <  2:687116/100‑3 (MQ=37)
                 tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATa                                                                                            <  1:2392408/100‑1 (MQ=37)
                   gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                         >  1:298616/1‑100 (MQ=255)
                        gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                                    >  1:2114859/1‑100 (MQ=255)
                            gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                <  1:1628590/100‑1 (MQ=255)
                            gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                <  2:723945/100‑1 (MQ=255)
                            gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                <  2:156242/100‑1 (MQ=255)
                             cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                               >  2:1863234/1‑100 (MQ=255)
                             cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                               <  1:1998917/100‑1 (MQ=255)
                             cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAAGACTTATCCGGCCTAC‑AAAATCGTGc                                                                                <  1:266293/99‑1 (MQ=255)
                                    acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCaa                                                                        <  1:748110/100‑1 (MQ=255)
                                    acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCaa                                                                        >  2:1107768/1‑100 (MQ=255)
                                     caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                       >  2:75776/1‑100 (MQ=255)
                                           cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                 >  1:1714925/1‑100 (MQ=255)
                                            aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                                >  2:213659/1‑100 (MQ=255)
                                            aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                                <  2:2680935/100‑1 (MQ=255)
                                             ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCa                                                               <  1:2901766/100‑1 (MQ=255)
                                                 gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAt                                                          <  1:2240157/100‑1 (MQ=255)
                                                    taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAt‑t‑‑‑tt                                                   >  2:1833524/1‑100 (MQ=255)
                                                    taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAt‑t‑‑‑tt                                                   >  1:2285169/1‑100 (MQ=255)
                                                      cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑g                                               <  1:767779/100‑1 (MQ=255)
                                                       aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑Gt                                              <  2:1101458/100‑1 (MQ=255)
                                                        gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTa                                             <  1:1066102/100‑1 (MQ=255)
                                                         aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAg                                            >  2:2234618/1‑100 (MQ=255)
                                                         aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAg                                            <  2:1931441/100‑1 (MQ=255)
                                                          ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAgg                                           <  1:2278694/100‑1 (MQ=255)
                                                            gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGcc                                         >  1:2088129/1‑100 (MQ=255)
                                                                aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGAt                                     >  2:1373498/1‑100 (MQ=255)
                                                                aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGAt                                     <  1:1023825/100‑1 (MQ=255)
                                                                 cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATa                                    >  1:2423367/1‑100 (MQ=255)
                                                                   aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAg                                  <  2:505799/100‑1 (MQ=255)
                                                                     gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCg                                >  1:8102/1‑100 (MQ=255)
                                                                       gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                              >  1:1134815/1‑100 (MQ=255)
                                                                       gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                              >  2:1992020/1‑100 (MQ=255)
                                                                        cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                             <  1:105123/100‑1 (MQ=255)
                                                                        cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                             <  1:2875506/100‑1 (MQ=255)
                                                                         gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                            >  1:458707/1‑100 (MQ=255)
                                                                          cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                           <  1:349317/100‑1 (MQ=255)
                                                                           attttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcgc                          >  2:1839977/10‑100 (MQ=18)
                                                                            ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCa                         <  2:2088136/100‑1 (MQ=255)
                                                                             ttcgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCGTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GGCT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCAt                        >  2:2311783/8‑100 (MQ=18)
                                                                              ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCA‑GGGA‑C‑‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATc                       >  1:2351461/7‑100 (MQ=18)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                       <  1:2447633/100‑1 (MQ=255)
                                                                               tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                      <  1:1512983/100‑1 (MQ=255)
                                                                                 taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                    >  2:2291973/4‑100 (MQ=21)
                                                                                 taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                    <  1:1396944/97‑1 (MQ=21)
                                                                                 taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                    >  2:1396946/4‑100 (MQ=21)
                                                                                 taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                    >  2:2519674/4‑100 (MQ=21)
                                                                                 taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GGAG‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                    >  2:1634854/4‑100 (MQ=18)
                                                                                 taggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                    >  1:1123923/4‑100 (MQ=18)
                                                                                  aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  2:2882801/3‑100 (MQ=18)
                                                                                  aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATGGGC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  1:2129657/3‑100 (MQ=17)
                                                                                  aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  1:1492369/3‑100 (MQ=17)
                                                                                  aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                   >  2:885089/3‑100 (MQ=18)
                                                                                   ggctggataaggcgTTCACGCCGCATCCAGCCTAC‑AAAATTGTGCAAAATCAATGGATTGCA‑CAGC‑T‑‑‑AGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  1:1699792/11‑100 (MQ=25)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:2741719/2‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑AAT‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAgatc                  >  1:1402919/2‑97 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGA‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:297573/2‑100 (MQ=17)
                                                                                    gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                 >  2:60616/1‑100 (MQ=255)
                                                                                    gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                 >  2:1317900/1‑100 (MQ=255)
                                                                                    gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATGCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑CTAGGCCTGATAAGCGGTGCGCATCAGGCAt                 >  2:1143506/1‑100 (MQ=255)
                                                                                    gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                >  1:1882686/1‑100 (MQ=18)
                                                                                    gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                 >  2:308360/1‑99 (MQ=17)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                 >  1:1597626/1‑99 (MQ=21)
                                                                                     ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                >  2:191299/9‑98 (MQ=17)
                                                                                     ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TGAG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                >  1:561257/9‑98 (MQ=17)
                                                                                     ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                >  1:1658508/9‑98 (MQ=21)
                                                                                     ccggataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                >  2:278676/9‑100 (MQ=18)
                                                                                     ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGCGCGCATCAGGCAt                 >  2:1286012/1‑99 (MQ=255)
                                                                                     ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                >  1:1448217/1‑100 (MQ=255)
                                                                                      cggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtt               >  2:2710258/8‑97 (MQ=18)
                                                                                      cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GGAG‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtt               >  2:2153981/8‑97 (MQ=17)
                                                                                       ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑AAT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt               >  1:2149795/7‑97 (MQ=21)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGAA‑GAGA‑T‑‑‑AAT‑‑GTAGGCCTGATAAGCGTAGCGAATCAGGCAattt              >  1:2273190/7‑100 (MQ=14)
                                                                                       ggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCA‑GGAA‑C‑‑‑GAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta              >  1:675479/7‑96 (MQ=25)
                                                                                        gataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt             >  1:1187139/6‑100 (MQ=17)
                                                                                        gataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑AAT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCAattt              >  2:2020099/6‑100 (MQ=21)
                                                                                        gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCA‑GGGC‑C‑‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagcttt             >  1:1840469/6‑95 (MQ=18)
                                                                                          tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttttt           >  1:1390426/1‑98 (MQ=255)
                                                                                           gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc          <  2:1448220/100‑1 (MQ=255)
                                                                                           gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc          >  1:1199982/1‑100 (MQ=255)
                                                                                            cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct         >  2:390814/1‑100 (MQ=255)
                                                                                            cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct         >  1:1906928/1‑100 (MQ=255)
                                                                                               cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg      <  1:1933571/100‑1 (MQ=255)
                                                                                                 ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc    <  2:1209265/100‑1 (MQ=255)
                                                                                                 ttAACGCCTTACCCGGCCTAC‑AAAATCGTGCATATGCAAGACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGTGCGCATCAGGCATTtttgcctttgcc    <  2:174846/100‑1 (MQ=21)
                                                                                                  tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa   <  2:463066/100‑1 (MQ=255)
                                                                                                   aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt  <  2:2408959/100‑1 (MQ=255)
                                                                                                   |                                                                                                           
GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705879‑3706077

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