Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,276,886 C→A 12.0% intergenic (+33/‑116) 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,8860CA12.0% 295.7 / 24.9 109intergenic (+33/‑116)garD/prlFD‑galactarate dehydrogenase/antitoxin of the SohA(PrlF)‑YhaV toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base C (47/48);  new base A (11/2);  total (58/51)
Fisher's exact test for biased strand distribution p-value = 1.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

Only 100 of 111 total aligned reads displayed.

GGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTGA  >  NC_000913/3276737‑3277035
                                                                                                                                                     |                                                                                                                                                     
ggCTGGAAGGTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCATTGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGcc                                                                                                                                                       <  2:83691/150‑1 (MQ=255)
 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)
                  tttATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTAt                                                                                                                                     <  2:384117/150‑1 (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)
                        ctcGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTa                                                                                                                               >  1:695530/1‑150 (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                                                                                                                 <  2:819048/149‑1 (MQ=255)
                                       gggAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTcc                                                                                                                 >  1:819048/1‑149 (MQ=255)
                                        ggAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGATCTGATTTCTTTTGGTACGTAGCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCAt                                                                                                               <  2:1226448/150‑1 (MQ=255)
                                         gaagaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCAtt                                                                                                              <  1:1530914/150‑1 (MQ=255)
                                           agaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaa                                                                                                            >  1:193733/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)
                                                        ctcGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatttgca                                                                                               <  2:1445427/150‑1 (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)
                                                                   gggggCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGAATTTGCATTACCTTAGtt                                                                                    >  2:618542/1‑150 (MQ=255)
                                                                    ggggCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTc                                                                                   >  1:1482717/1‑150 (MQ=255)
                                                                       gCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAg                                                                                <  2:819131/150‑1 (MQ=255)
                                                                         tGCTTAACCAGCTGGCGGTGATTAACCCGGCACCGGTGACCTGATTTCTTTTGTAACGTCCCTCGCAAAACATGGCTTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCt                                                                              <  1:62027/150‑1 (MQ=255)
                                                                         tGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCt                                                                              <  2:1068133/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)
                                                                                    cTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAAGGTCCCTCGCAAAACATGGCCTTAGTTCCATGTTTTTATTGTTTAAAGCCCCCACTTCCATTTATAATCCATTTGCAATACCTTAGTTCAAGCTTATAATTTgtg                                                                   <  1:43287/150‑3 (MQ=255)
                                                                                    cTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAg                                                                   <  1:540058/150‑1 (MQ=255)
                                                                                     tGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGc                                                                  >  1:1437203/1‑150 (MQ=255)
                                                                                     tGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGc                                                                  <  2:605958/150‑1 (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)
                                                                                                 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)
                                                                                                                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)
                                                                                                                                      caataaaaaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcgct                 >  1:1668738/6‑150 (MQ=255)
                                                                                                                                      cTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGTTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGGGAATCATGCCCGCTAATgcacgct                 >  2:57056/1‑150 (MQ=255)
                                                                                                                                        ataaaaaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGctct               >  2:1216936/4‑150 (MQ=255)
                                                                                                                                        ataaaaaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGctct               >  1:1332515/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:1089952/1‑150 (MQ=255)
                                                                                                                                           aaaaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCAc            >  2:1310236/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)
                                                                                                                                                     |                                                                                                                                                     
GGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTGA  >  NC_000913/3276737‑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: