Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,276,887 T→C 11.9% intergenic (+34/‑115) garD → / → prlF D‑galactarate dehydrogenase/antitoxin of the SohA(PrlF)‑YhaV toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,276,8870TC11.9% 278.0 / 26.1 109intergenic (+34/‑115)garD/prlFD‑galactarate dehydrogenase/antitoxin of the SohA(PrlF)‑YhaV toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (47/49);  new base C (11/2);  total (58/51)
Fisher's exact test for biased strand distribution p-value = 1.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 110 total aligned reads displayed.

GCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTGA  >  NC_000913/3276738‑3277035
                                                                                                                                                     |                                                                                                                                                    
gCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCt                                                                                                                                                      >  2:1227941/1‑150 (MQ=255)
   ggAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAg                                                                                                                                                   >  1:496450/1‑150 (MQ=255)
     aaGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGTAGAAGAAAACCTTCTCGGATCTATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTg                                                                                                                                                 <  1:1347006/150‑1 (MQ=255)
     aaGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTg                                                                                                                                                 <  2:1403370/150‑1 (MQ=255)
        ttgttCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCa                                                                                                                                              <  1:1446681/150‑1 (MQ=255)
          gttCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATg                                                                                                                                            <  1:842728/150‑1 (MQ=255)
            tCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGtt                                                                                                                                          <  2:552278/150‑1 (MQ=255)
              cACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCACTAAGGCCATGtttt                                                                                                                                        <  2:1188249/150‑1 (MQ=255)
                cTTTATTCTCGACGTCGCCAGCGGCAAGAAGAAAACCTTCTCGGATCAATGGGGCCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTCATTTCTTTAGCAACGACCCACGCAAAACATGCCCTAATTGCCAAGttttaa                                                                                                                                      >  2:694885/1‑148 (MQ=255)
                 tttATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACATCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTAt                                                                                                                                     <  2:204955/150‑1 (MQ=255)
                 tttATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCGGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGATCGTCCCTCGCAAAACATGGCACTAAGGCCATGttttgcg                                                                                                                                     <  1:506492/150‑4 (MQ=255)
                 attATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGGTTAACCCGGCTCCGGTGACCTGATTTCGTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCTTGTTTTTAt                                                                                                                                     >  1:193529/2‑150 (MQ=255)
                  ttATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTAtt                                                                                                                                    >  1:1097930/1‑150 (MQ=255)
                   tATTCTCGCCGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGACCGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTg                                                                                                                                   >  1:127418/1‑150 (MQ=255)
                     ttCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGtt                                                                                                                                 <  2:925373/150‑1 (MQ=255)
                     ttCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGtt                                                                                                                                 <  2:500263/150‑1 (MQ=255)
                      tctcGACGTCGCCAGCGGGAAGAAGAACACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGttt                                                                                                                                <  1:638205/150‑1 (MQ=255)
                       ctcGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACAAGGCCTTAGTGCCATGTTTTTATTGTTTa                                                                                                                               >  1:1525594/1‑150 (MQ=255)
                        tcGCCGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTaa                                                                                                                              <  2:127418/150‑1 (MQ=255)
                          gACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAg                                                                                                                            <  2:1526334/150‑1 (MQ=255)
                               cGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATATCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGccccc                                                                                                                       >  2:1517319/1‑150 (MQ=255)
                                gCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCa                                                                                                                      <  1:1227941/150‑1 (MQ=255)
                                    gCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTc                                                                                                                  <  1:635599/150‑1 (MQ=255)
                                    gCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTc                                                                                                                  <  1:514675/150‑1 (MQ=255)
                                      gggAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTcc                                                                                                                 >  1:819048/1‑149 (MQ=255)
                                      gggAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTcc                                                                                                                 <  2:819048/149‑1 (MQ=255)
                                       ggAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCAt                                                                                                               >  2:1318738/1‑150 (MQ=255)
                                        gaagaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCAtt                                                                                                              <  1:1530914/150‑1 (MQ=255)
                                          agaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaa                                                                                                            >  1:193733/1‑150 (MQ=255)
                                          agaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACAGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGGTTAAAGCCCCCCCGTCCATtaa                                                                                                            >  1:999412/1‑150 (MQ=255)
                                          agaagaaAACCTTCTCGGATCAATGGGGGCTGCATAAACAGCTGGCGGTGTTTAACCCGGCCCCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAATCCCACACGTCCATtaa                                                                                                            >  1:819435/1‑150 (MQ=255)
                                           gaagaaTACCTTCTCGGTGCACTGTGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACAGGTGACCTGATTTGTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGACATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaat                                                                                                           <  1:224612/150‑1 (MQ=255)
                                            aagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaata                                                                                                          >  1:1068133/1‑150 (MQ=255)
                                             agaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaataa                                                                                                         >  2:1383491/1‑150 (MQ=255)
                                               aaatCCTTCTCGGATCAATTGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCTGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtg                                                                                                       <  2:57498/150‑1 (MQ=255)
                                               aaaaCCTTCTCGGATGAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtg                                                                                                       <  1:685297/150‑1 (MQ=255)
                                                   ccTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatt                                                                                                   >  2:1489677/1‑150 (MQ=255)
                                                   ccTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatt                                                                                                   >  1:491015/1‑150 (MQ=255)
                                                     ttCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatttg                                                                                                 <  1:999454/150‑1 (MQ=255)
                                                      tctcGGATCAATGTGTGCTGCATATCCATCTGGCGGTGTTTAACCCGGAACTGTTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgctaatgc                                                                                                <  2:1165844/150‑7 (MQ=255)
                                                         cGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatttgcatt                                                                                             >  2:718949/1‑150 (MQ=255)
                                                           gATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTAc                                                                                           >  2:501382/1‑150 (MQ=255)
                                                           gATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTAc                                                                                           >  2:1392535/1‑150 (MQ=255)
                                                            aTCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTAcc                                                                                          <  2:491015/150‑1 (MQ=255)
                                                              cAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTCAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCtt                                                                                        <  2:1133759/150‑1 (MQ=255)
                                                              cAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCtt                                                                                        <  1:1396067/150‑1 (MQ=255)
                                                                   ggggCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTc                                                                                   >  1:1482717/1‑150 (MQ=255)
                                                                    gggCTGCATAACCAGCTGGCGGTGGTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCa                                                                                  <  2:1565465/150‑1 (MQ=255)
                                                                      gCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAg                                                                                <  2:819131/150‑1 (MQ=255)
                                                                        tGCTTAACCAGCTGGCGGTGATTAACCCGGCACCGGTGACCTGATTTCTTTTGTAACGTCCCTCGCAAAACATGGCTTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCt                                                                              <  1:62027/150‑1 (MQ=255)
                                                                        tGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCt                                                                              <  2:1023768/150‑1 (MQ=255)
                                                                           aTAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTtat                                                                           <  1:402014/150‑1 (MQ=255)
                                                                               ccAGCTGGCGGTGTTTAACCAGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTATTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAAtt                                                                       <  1:348983/150‑1 (MQ=255)
                                                                                  gCTGGCGGTGTTTACCCCGGCACCGGTGACCTGATTTCGTGTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAGAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTGAGTTCACGCTTATAATTTgc                                                                    >  1:1068410/1‑149 (MQ=255)
                                                                                  gCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGATCGTCCATCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCCCGTCCATTAAAAAAGCATTTGCATTACCTTAGTACAAGCTTATAATTTGa                                                                    >  1:1062720/1‑150 (MQ=255)
                                                                                  gCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATATAGCAATTGCATTACCTTAGTTCCAGCGTATAATTAGa                                                                    >  2:1426405/1‑150 (MQ=255)
                                                                                  gCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTTCCTTAGTTCAAGCTTATAATTTGa                                                                    >  1:235648/1‑150 (MQ=255)
                                                                                  gCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCAACGTCCATAAATAATGCATTAGAATTAACTAAATTCAATCTtatattataa                                                                    >  2:1055473/1‑144 (MQ=255)
                                                                                   cTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAg                                                                   <  1:540058/150‑1 (MQ=255)
                                                                                    tGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGc                                                                  <  2:605958/150‑1 (MQ=255)
                                                                                    tGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGc                                                                  >  1:1437203/1‑150 (MQ=255)
                                                                                     ggaggTGTTTAACCCGGCACCGGTGACCTTATTACTTTTGGAACGTCCCTGGCAAAAAATTGCCTTAGTGCCATGTTTTTATTGTTAAAAGCGCCCACGTCCATTAATAAAGCATATGCATTACCTTAGTTCAAGCTTATAATTTGAGCa                                                                 <  2:92675/147‑1 (MQ=255)
                                                                                       cGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCACCCACGTCCATTTATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAgc                                                               >  2:1265548/1‑149 (MQ=255)
                                                                                          tgtTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGaaaa                                                            >  2:1175749/1‑150 (MQ=255)
                                                                                           gtTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAAc                                                           >  1:1141219/1‑150 (MQ=255)
                                                                                           gtTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAAc                                                           <  2:1653549/150‑1 (MQ=255)
                                                                                           gtTTAACCCGGCACCGGTGACCTGATATCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAAc                                                           <  1:1517319/150‑1 (MQ=255)
                                                                                             ttAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAACCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAg                                                         <  1:308877/150‑1 (MQ=255)
                                                                                               aaCCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGc                                                       <  2:496450/150‑1 (MQ=255)
                                                                                                aCCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCt                                                      >  2:587195/1‑150 (MQ=255)
                                                                                                 cccGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTg                                                     <  2:1097930/150‑1 (MQ=255)
                                                                                                 cccGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGGTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTAAAGCTGTTAATTTGAGCAGAAAACAGGATg                                                     >  1:320059/1‑150 (MQ=255)
                                                                                                  ccGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGt                                                    >  1:148673/1‑150 (MQ=255)
                                                                                                       aCCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAg                                               <  2:193733/150‑1 (MQ=255)
                                                                                                        ccGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTTTAATTTGAGCAGAAAACAGGCTGTAAAAgg                                              >  2:187160/1‑150 (MQ=255)
                                                                                                          ggTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGAc                                            >  2:812575/1‑150 (MQ=255)
                                                                                                           gTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACa                                           <  1:379653/150‑1 (MQ=255)
                                                                                                               ccTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAATCAGGCTGTATAAGGACAGTGa                                       >  2:206265/1‑150 (MQ=255)
                                                                                                                   aTTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCa                                   <  2:1068531/150‑1 (MQ=255)
                                                                                                                    tttCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTAAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGCCAGTGAATCAt                                  >  2:835289/1‑150 (MQ=255)
                                                                                                                        ttttGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGccc                              >  1:912914/1‑150 (MQ=255)
                                                                                                                           tGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCt                           <  1:1392535/150‑1 (MQ=255)
                                                                                                                               aCGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAATCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGATCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATg                       <  2:819435/150‑1 (MQ=255)
                                                                                                                                 gTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgct                     >  2:480910/1‑150 (MQ=255)
                                                                                                                                   aacaataaaaaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCAGTAGCTTAGTTCAAGCTTATTATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATTCCCGCTTATgctcg                   >  1:879025/8‑150 (MQ=255)
                                                                                                                                     caataaaaaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATTATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcgct                 >  2:206939/6‑150 (MQ=255)
                                                                                                                                     cTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGTTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGGGAATCATGCCCGCTAATgcacgct                 >  2:57056/1‑150 (MQ=255)
                                                                                                                                       ataaaaaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGctct               >  1:1332515/4‑150 (MQ=255)
                                                                                                                                       ataaaaaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGctct               >  2:1216936/4‑150 (MQ=255)
                                                                                                                                        taaaaaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCtctc              >  2:1172346/3‑150 (MQ=255)
                                                                                                                                        gCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCtctc              >  1:1284912/1‑150 (MQ=255)
                                                                                                                                         cAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCa             >  2:1353313/1‑150 (MQ=255)
                                                                                                                                         aaaaaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAATGACAGTGAATCATGCCCGCTAATGCTCGCTCTCa             >  2:674139/2‑150 (MQ=255)
                                                                                                                                          aaaaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGATTCATGCCCGCTAATGCTAGCTCTCAc            >  1:207409/1‑150 (MQ=255)
                                                                                                                                          aaaaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCAc            >  2:1310236/1‑150 (MQ=255)
                                                                                                                                          aaaaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCAc            >  2:1089952/1‑150 (MQ=255)
                                                                                                                                          aaaaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCACTAATGCTCGCTCTCAc            >  2:1077699/1‑150 (MQ=255)
                                                                                                                                           aaaCATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACg           >  1:639487/1‑150 (MQ=255)
                                                                                                                                                tGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCACCCACGTCCATTAATAATGCATTTGCATTACATTATTTCAAGCTTATAATTTGAGCAGAAAACAGGATGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTa      >  1:969302/1‑150 (MQ=255)
                                                                                                                                                    cTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTGa  >  1:165945/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTGA  >  NC_000913/3276738‑3277035

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