Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,659,467 A→C 7.2% E79A (GAA→GCA) ‡ mdtE → anaerobic multidrug efflux transporter, ArcA‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,659,4670AC 7.2% 405.3 / 16.2 138E79A (GAA→GCA) ‡mdtEanaerobic multidrug efflux transporter, ArcA‑regulated
Reads supporting (aligned to +/- strand):  ref base A (62/65);  new base C (5/5);  total (67/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 139 total aligned reads displayed.

CCGCCATGACGCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCAGGCATCGTTGC  >  NC_000913/3659320‑3659616
                                                                                                                                                   |                                                                                                                                                     
ccgccATGACGCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGCGCGGTATTATCATTAAACTCAACTTTATCGAAg                                                                                                                                                     >  1:1342171/1‑150 (MQ=255)
   ccATTACGCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTAGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCTTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGGGGGCGGTATTATCATTAACCGCAACTTTATCGAAGGCg                                                                                                                                                  >  2:1481484/1‑150 (MQ=255)
    cATGACGCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGa                                                                                                                                                 >  2:932359/1‑150 (MQ=255)
      tGACGCCTGAGGTCTGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCAACAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATa                                                                                                                                               <  2:1838997/150‑1 (MQ=255)
       gACGCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATaa                                                                                                                                              >  2:581845/1‑150 (MQ=255)
        aCGCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCATCTTTATCGAAGGCGATaaa                                                                                                                                             >  1:102701/1‑150 (MQ=255)
          gCCTGAGCGCGGTGTCGGCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGt                                                                                                                                           <  1:336520/150‑1 (MQ=255)
           ccTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGGCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGCCTTCGATAAAGtt                                                                                                                                          <  2:1683164/150‑2 (MQ=255)
            cTGAGGTCGGTGTCGTCACACTCTCCCCGGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATTCGTCCCCAGGTTGTCGGTATTATCATTAAACGCAACTTTATCGCCTTCGATAAAGttg                                                                                                                                         <  2:900714/150‑3 (MQ=255)
            cTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGCCTTCGATAAAGttg                                                                                                                                         <  1:2043069/150‑3 (MQ=255)
            cTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGGGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGCCTTCGATAAAGttg                                                                                                                                         <  1:110284/150‑3 (MQ=255)
                ggTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCa                                                                                                                                     <  2:1098868/150‑1 (MQ=255)
                   cGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAggg                                                                                                                                  >  2:2019159/1‑150 (MQ=255)
                    ggTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGc                                                                                                                                 <  1:1097023/150‑1 (MQ=255)
                      tgtCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGa                                                                                                                               >  2:905463/1‑150 (MQ=255)
                           tcACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCga                                                                                                                          >  2:1451728/1‑149 (MQ=255)
                           tcACACTCTCCCCCGGGTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGACGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGc                                                                                                                          <  2:770491/150‑1 (MQ=255)
                            cacaCTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCt                                                                                                                         <  1:1631412/150‑1 (MQ=255)
                              cacTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGt                                                                                                                       <  1:100463/150‑1 (MQ=255)
                               acTCTCGCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCAGGTAGAACCGTTCCTTATGAAGTTGTCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTa                                                                                                                      <  2:1950890/150‑1 (MQ=255)
                               acTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTa                                                                                                                      <  2:1083605/150‑1 (MQ=255)
                               acTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTa                                                                                                                      >  1:1153155/1‑150 (MQ=255)
                               acTCTCCCCCGGTGCGGGCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTa                                                                                                                      <  1:543443/150‑1 (MQ=255)
                                ctctCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTAt                                                                                                                     >  1:1483424/1‑150 (MQ=255)
                                 tctcCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATc                                                                                                                    >  1:2049035/1‑150 (MQ=255)
                                  ctcCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCa                                                                                                                   >  1:48617/1‑150 (MQ=255)
                                      cccGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGGGGGCGGTATTATCATTAAACTCAACTTTATCGAAGGCTATAAAGTGAACCAGGGCGATTCGCTGTATCAGAtt                                                                                                               <  2:2078928/150‑1 (MQ=255)
                                         ggTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGTTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGAt                                                                                                            <  1:1473368/150‑1 (MQ=255)
                                         ggTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGAt                                                                                                            <  1:1888624/150‑1 (MQ=255)
                                           ttCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATcc                                                                                                          <  1:522306/150‑1 (MQ=255)
                                           gtCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATcc                                                                                                          <  1:1671822/149‑1 (MQ=255)
                                            tCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCt                                                                                                         <  1:1632470/150‑1 (MQ=255)
                                               gTCAATGTGTGGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATGCGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGACGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCa                                                                                                      <  1:755591/150‑1 (MQ=255)
                                                 cAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCAcc                                                                                                    <  2:467660/150‑1 (MQ=255)
                                                 cAATGTGTGGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCAcc                                                                                                    <  2:541252/150‑1 (MQ=255)
                                                        ttGAGCGAATTACCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAg                                                                                             <  1:1478562/150‑1 (MQ=255)
                                                         tGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAgg                                                                                            >  2:1891874/1‑150 (MQ=255)
                                                          gAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGc                                                                                           >  1:1260771/1‑150 (MQ=255)
                                                            gCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCg                                                                                         >  2:2019193/1‑150 (MQ=255)
                                                            gCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCg                                                                                         <  1:1910989/150‑1 (MQ=255)
                                                             cGAATTGCCCGTTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCAAGGGGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGa                                                                                        <  2:391835/150‑1 (MQ=255)
                                                                 ttGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTa                                                                                    <  1:1800859/150‑1 (MQ=255)
                                                                 ttGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGAACCTTTACAGGCCGAGCTa                                                                                    >  1:1347111/1‑150 (MQ=255)
                                                                   gCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTaaa                                                                                  >  1:766040/1‑150 (MQ=255)
                                                                    cccTGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCGAGGTGGGCGGTTTTATCATTAAACTCAACTTTTTCGAAGGCGTTAAATTGAACCAGGGGGATTCTCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAAc                                                                                 <  1:1525718/150‑1 (MQ=255)
                                                                        gTAGAACCGGTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCCGCACCTTTACAGGCCGAGCTAAACTCCg                                                                             <  2:71716/150‑1 (MQ=255)
                                                                         tAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATATAGTGAACCAGGGCGCTTCGGTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAACCTCCGc                                                                            >  1:53372/1‑150 (MQ=255)
                                                                           gAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCa                                                                          <  1:1776574/150‑1 (MQ=255)
                                                                            aaCCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCaa                                                                         >  2:2069701/1‑150 (MQ=255)
                                                                             aCCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCaaa                                                                        >  2:1420369/1‑150 (MQ=255)
                                                                                gTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGc                                                                     >  1:2092416/1‑150 (MQ=255)
                                                                                    cTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGc                                                                 >  1:1832818/1‑150 (MQ=255)
                                                                                    cTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGc                                                                 >  1:1446563/1‑150 (MQ=255)
                                                                                     ttATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCt                                                                <  1:184487/150‑1 (MQ=255)
                                                                                       aTGAGGTTGCCGAGATACGTCCCCAGGTGTGCTGTATTATCATTAAACGCAACTTTATCGAAGGAGATAAAGTGAAGCAGGTCGAGTCGCTGTATCAGATTGATCCTGCACCTTTTCAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTgg                                                              <  1:343168/150‑1 (MQ=255)
                                                                                        tGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTAGCTGGc                                                             >  1:184588/1‑150 (MQ=255)
                                                                                        tGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAAACAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGATCGCTGGc                                                             >  1:20817/1‑150 (MQ=255)
                                                                                         gAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCg                                                            <  2:48617/150‑1 (MQ=255)
                                                                                         gAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCg                                                            <  1:1015960/150‑1 (MQ=255)
                                                                                             ttGCCGTGATACGTCGCCAGGTGGGCGGTTTTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACGTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAg                                                        <  2:380625/150‑1 (MQ=255)
                                                                                              tGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAgc                                                       <  2:1317015/150‑1 (MQ=255)
                                                                                               tcagagaTAAGTCCGCAGTTGGGCTGTATTATCATTAAACGCAACTTTATCGTATGCTATAAAGTGAACCAGTGCTATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAgcg                                                      <  1:163558/147‑1 (MQ=255)
                                                                                               gCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAgcg                                                      >  2:2102184/1‑150 (MQ=255)
                                                                                                 cGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGCGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGct                                                    >  2:544004/1‑150 (MQ=255)
                                                                                                 cGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGct                                                    >  1:952158/1‑150 (MQ=255)
                                                                                                   agaTACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGctct                                                  <  1:888474/150‑1 (MQ=255)
                                                                                                    gaTACGGCCCCAGGTGGGCGGTATTATCATTAACCGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGctctc                                                 <  2:687885/150‑1 (MQ=255)
                                                                                                      tACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTa                                               >  2:1429867/1‑150 (MQ=255)
                                                                                                      tACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTa                                               >  1:1217977/1‑150 (MQ=255)
                                                                                                         gTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCg                                            <  2:132020/150‑1 (MQ=255)
                                                                                                          tCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGc                                           >  2:803120/1‑150 (MQ=255)
                                                                                                           ccccAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGcc                                          >  2:793652/1‑150 (MQ=255)
                                                                                                            cccAGGGGGGCGGTATTATCATTAAACGCAACATTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCa                                         <  2:175740/150‑1 (MQ=255)
                                                                                                                  tGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATg                                   <  1:1645232/150‑1 (MQ=255)
                                                                                                                    ggCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGcc                                 >  2:1472821/1‑150 (MQ=255)
                                                                                                                            tatCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCAc                         >  2:1651319/1‑150 (MQ=255)
                                                                                                                             atCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCAcc                        <  2:2067810/150‑1 (MQ=255)
                                                                                                                               cATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCTGATTGATCCTGCACCTTTACAGGCCGCGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCtt                      >  2:1274112/1‑150 (MQ=255)
                                                                                                                               cATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCtt                      <  2:1124890/150‑1 (MQ=255)
                                                                                                                               cATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCtt                      <  1:2053459/150‑1 (MQ=255)
                                                                                                                               cATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCtt                      >  2:1748127/1‑150 (MQ=255)
                                                                                                                                aTTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCttt                     >  2:1469898/1‑150 (MQ=255)
                                                                                                                                 ttAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATGCGCGGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTa                    <  1:391606/150‑1 (MQ=255)
                                                                                                                                   aaaCGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCTATTCGCTGTATCAGATTGATCCTGCACCTTTAAAGGCCGAGCAAAACTCCGCCAAACGCTCGATGGCGAAAGCGCGCTCTACCGCCAGCAATGCCCGCATCACCTTTacc                  >  2:1241795/1‑148 (MQ=255)
                                                                                                                                   aaaCGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAAc                  >  2:954272/1‑150 (MQ=255)
                                                                                                                                   aaaCGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAAc                  <  1:1573719/150‑1 (MQ=255)
                                                                                                                                    tggttcaCTTTATCGCCTTCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTTCAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAAcc                 >  2:82550/7‑150 (MQ=255)
                                                                                                                                    tggttcaCTTTATCGCCTTCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAACGGCTCGCTGGCGAACGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAAcc                 >  2:1085460/7‑150 (MQ=255)
                                                                                                                                     ggttcaCTTTATCGCCTTCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCg                >  2:681102/6‑150 (MQ=255)
                                                                                                                                     aCGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCAGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCg                >  1:729460/1‑150 (MQ=255)
                                                                                                                                       gCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGcc              <  2:2077651/150‑1 (MQ=255)
                                                                                                                                        tcaCTTTATCGCCTTCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGTCCGACCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAAACGCCa             >  2:878374/3‑150 (MQ=255)
                                                                                                                                        cAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGGGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCa             >  2:1207817/1‑150 (MQ=255)
                                                                                                                                        cAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCa             >  1:1738247/1‑150 (MQ=255)
                                                                                                                                             ttATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAAAGCCCGCATCACCTTTACCCGCCCGGCAt        >  2:1266117/1‑150 (MQ=255)
                                                                                                                                              tATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGAATAAACCTGCACCTTTACAGGCCGAGCTGAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTCCCGCCACAAATGCCCGCCTCACCTTTAACCGCCAGGCATc       >  1:109139/1‑150 (MQ=255)
                                                                                                                                               aTCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCAGGCATCg      >  1:1482547/1‑150 (MQ=255)
                                                                                                                                                 cGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTAACAGATTGATCCTGCACCTTTACACACCGCGCTAAACTCCGCCAAATGCTCGCTGGCGAACGCGCTCTCGACAGCCAGCAATGCCCGCAGCACCTTTAACCGCCAGGCAACGtt    >  2:329706/1‑150 (MQ=255)
                                                                                                                                                  gAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCAGGCATCGTTg   >  2:1654405/1‑150 (MQ=255)
                                                                                                                                                  gAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCAGGCATCGTTg   >  1:470673/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCGCCATGACGCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCAGGCATCGTTGC  >  NC_000913/3659320‑3659616

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