Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,276,886 C→A 5.3% 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,8860CA 5.3% 770.7 / 20.3 218intergenic (+33/‑116)garD/prlFD‑galactarate dehydrogenase/antitoxin of the SohA(PrlF)‑YhaV toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base C (102/104);  new base A (4/8);  total (106/112)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 218 total aligned reads displayed.

TGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTAC  >  NC_000913/3276740‑3277032
                                                                                                                                                   |                                                                                                                                                  
  gAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAg                                                                                                                                                >  4:1549901/1‑149 (MQ=255)
  gAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAg                                                                                                                                                >  3:916487/1‑149 (MQ=255)
     gttgttCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGc                                                                                                                                             >  4:1709147/1‑149 (MQ=255)
      ttgttCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGcc                                                                                                                                            <  1:1203734/149‑1 (MQ=255)
          tCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGt                                                                                                                                        >  1:426462/1‑149 (MQ=255)
            cACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGttt                                                                                                                                      >  2:1442122/1‑149 (MQ=255)
            cACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGttt                                                                                                                                      <  2:1712778/149‑1 (MQ=255)
            cACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑ACTAAGGCCATGttt                                                                                                                                      <  2:1263556/149‑1 (MQ=255)
             aCTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGtttt                                                                                                                                     >  2:1131471/1‑149 (MQ=255)
               tttATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTa                                                                                                                                   <  3:84110/149‑1 (MQ=255)
                 tATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTAtt                                                                                                                                 <  3:316428/149‑1 (MQ=255)
                 tATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTAtt                                                                                                                                 <  1:315129/149‑1 (MQ=255)
                  aTTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTg                                                                                                                                >  2:1100153/1‑149 (MQ=255)
                  aTTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑ACTAAGGCCATGttttgcgag                                                                                                                                <  3:831397/149‑6 (MQ=255)
                      tcGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑Ct                                                                                                                                                   >  2:1530543/1‑126 (MQ=255)
                      tcGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑Ct                                                                                                                                                   <  1:1530543/126‑1 (MQ=255)
                      tcGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTa                                                                                                                            >  4:647764/1‑149 (MQ=255)
                        gACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTaaa                                                                                                                          <  1:1122811/149‑1 (MQ=255)
                        gACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTaaa                                                                                                                          <  3:1254886/149‑1 (MQ=255)
                           gTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGcc                                                                                                                       <  1:1337/149‑1 (MQ=255)
                              gCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTATGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGcccac                                                                                                                    >  4:1105042/1‑147 (MQ=255)
                                   cGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTc                                                                                                               <  1:1442122/149‑1 (MQ=255)
                                      gaagaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCAt                                                                                                            >  1:892409/1‑149 (MQ=255)
                                        agaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATta                                                                                                          >  1:1648843/1‑149 (MQ=255)
                                             aaaaCCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTaataat                                                                                                     <  2:676924/149‑1 (MQ=255)
                                                 ccTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcat                                                                                                 >  1:488610/1‑149 (MQ=255)
                                                  cTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcant                                                                                                <  3:645601/149‑1 (MQ=255)
                                                   ttCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcattt                                                                                               <  4:1515867/149‑1 (MQ=255)
                                                   ttCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcattt                                                                                               <  3:1436877/149‑1 (MQ=255)
                                                   ttCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGCGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcattt                                                                                               <  4:293944/149‑1 (MQ=255)
                                                   ttCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑ACTAAGGCCATGttt                                                                                                                                      <  3:534849/110‑1 (MQ=255)
                                                   ttCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑ACTAAGGCCATGttt                                                                                                                                      >  4:534849/1‑110 (MQ=255)
                                                    tctcGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatttg                                                                                              <  2:565089/149‑1 (MQ=255)
                                                     ctcGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatttgc                                                                                             <  2:1339156/149‑1 (MQ=255)
                                                      tcGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGc                                                                                                                                             <  1:367255/100‑1 (MQ=255)
                                                       cGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatttgcat                                                                                           <  2:1231783/149‑1 (MQ=255)
                                                               tGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCAAATGCATTattaatg                                                                                   <  1:785822/149‑7 (MQ=255)
                                                                   ggCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCa                                                                               >  1:205613/1‑149 (MQ=255)
                                                                    gCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCaa                                                                              <  1:56442/149‑1 (MQ=255)
                                                                      tGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTaataat                                                                                                     >  3:1285230/1‑124 (MQ=255)
                                                                      tGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTaataat                                                                                                     <  4:1285230/124‑1 (MQ=255)
                                                                        cATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTTAAGCtt                                                                          <  2:190362/149‑1 (MQ=255)
                                                                            aCCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATaa                                                                      <  2:1215162/149‑1 (MQ=255)
                                                                              cAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAAtt                                                                    >  1:1101793/1‑149 (MQ=255)
                                                                               aGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAAttt                                                                   <  3:136467/149‑1 (MQ=255)
                                                                                 cTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGa                                                                 <  4:916487/149‑1 (MQ=255)
                                                                                  tGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAg                                                                <  4:784449/149‑1 (MQ=255)
                                                                                      ggTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGa                                                            >  3:978207/1‑149 (MQ=255)
                                                                                          tttAACCCGGCACCGATGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAAc                                                        >  2:638252/1‑149 (MQ=255)
                                                                                           ttAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACa                                                       <  1:705799/149‑1 (MQ=255)
                                                                                             aaCCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAgg                                                     <  2:1689597/149‑1 (MQ=255)
                                                                                             aaCCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTCATAATTTGAGCAGAAAACAgg                                                     <  1:955950/149‑1 (MQ=255)
                                                                                              aCCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGa                                                                 <  2:833608/136‑1 (MQ=255)
                                                                                              aCCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGc                                                    <  2:576432/149‑1 (MQ=255)
                                                                                               cccGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAAttt                                                                   >  3:719166/1‑133 (MQ=255)
                                                                                               cccGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCt                                                   >  1:1705771/1‑149 (MQ=255)
                                                                                                ccGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTg                                                  >  1:150680/1‑149 (MQ=255)
                                                                                                          tGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACa                                        <  2:365082/149‑1 (MQ=255)
                                                                                                          tGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACa                                        <  3:59622/149‑1 (MQ=255)
                                                                                                            aCCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAg                                                      >  4:189518/1‑133 (MQ=255)
                                                                                                            aCCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAg                                                      <  3:189518/133‑1 (MQ=255)
                                                                                                            aCCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGt                                      <  1:1291838/149‑1 (MQ=255)
                                                                                                              cTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTta                                                                         <  1:855723/112‑1 (MQ=255)
                                                                                                              cTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGa                                    >  2:862320/1‑149 (MQ=255)
                                                                                                              cTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGa                                    >  1:240886/1‑149 (MQ=255)
                                                                                                               tGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTaaa                                                                                                                          >  3:1244486/1‑62 (MQ=255)
                                                                                                               tGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTaaa                                                                                                                          <  4:1244486/62‑1 (MQ=255)
                                                                                                               tGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGaa                                   <  3:1103421/149‑1 (MQ=255)
                                                                                                                gATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAAt                                  >  1:664407/1‑149 (MQ=255)
                                                                                                                 aTTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATc                                 >  3:1111805/1‑149 (MQ=255)
                                                                                                                 aTTGCTTTTGGAACGTCCCTCGCAAAACATGGC‑ATTAGTGCCATGTTTTTATTCTTTAAAGCACCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATc                                 <  2:296988/149‑1 (MQ=255)
                                                                                                                   ttCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGc                                                                            >  2:39891/1‑104 (MQ=255)
                                                                                                                   ttCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCAt                               <  1:1216828/149‑1 (MQ=255)
                                                                                                                   ttCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCAt                               <  1:1012226/149‑1 (MQ=255)
                                                                                                                      ttttGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGcc                            >  3:1196702/1‑149 (MQ=255)
                                                                                                                      ttttGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGcc                            >  2:236235/1‑149 (MQ=255)
                                                                                                                      ttttGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGcc                            >  1:51160/1‑149 (MQ=255)
                                                                                                                       tttGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGt                                                                                  >  3:479011/1‑94 (MQ=255)
                                                                                                                       tttGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGt                                                                                  <  4:479011/94‑1 (MQ=255)
                                                                                                                       tttGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATAccc                           >  3:1455114/1‑149 (MQ=255)
                                                                                                                           gAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTa                                                                                                                                   <  1:1396580/41‑1 (MQ=255)
                                                                                                                           gAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTa                                                                                                                                   >  2:1396580/1‑41 (MQ=255)
                                                                                                                              cGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATg                    <  2:1032141/149‑1 (MQ=255)
                                                                                                                                tCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgct                  <  1:1752579/149‑1 (MQ=255)
                                                                                                                                 cccTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCACAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctc                 >  3:1473188/1‑149 (MQ=255)
                                                                                                                                   cTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcgc               >  3:456369/1‑149 (MQ=255)
                                                                                                                                   cTCGCAAAACATGGC‑CTTAGTGCCATGGTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcgc               >  3:1645094/1‑149 (MQ=255)
                                                                                                                                    tCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcgct              >  4:577040/1‑149 (MQ=255)
                                                                                                                                     cGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcgctc             <  3:1549901/149‑1 (MQ=255)
                                                                                                                                     ataaaaaCATGGCACTAAG‑GCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTa                                                                                    <  2:1323711/75‑1 (MQ=25)
                                                                                                                                     ataaaaaCATGGC‑ACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTa                                                                                    >  1:1323711/4‑78 (MQ=25)
                                                                                                                                      gCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGctct            >  4:1381844/1‑149 (MQ=255)
                                                                                                                                        aaaaCATGGCACTAAG‑GCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCa          >  1:1227991/1‑149 (MQ=255)
                                                                                                                                        aaaaCATGGCACTAAG‑GCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCa          >  4:816624/1‑149 (MQ=255)
                                                                                                                                        aaaaCATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAg                                                                                   >  3:1055816/1‑76 (MQ=255)
                                                                                                                                        aaaaCATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAg                                                                                   <  4:1055816/76‑1 (MQ=255)
                                                                                                                                        aaaaCATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCa          <  3:345605/149‑1 (MQ=255)
                                                                                                                                         aaaCATGGCACTAAG‑GCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCAc         >  4:1187232/1‑149 (MQ=255)
                                                                                                                                         aaaCATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCAc         <  1:106635/149‑1 (MQ=255)
                                                                                                                                                gC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGc                                                    >  2:49591/1‑99 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTAC  >  NC_000913/3276740‑3277032

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑