Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,276,890 G→A 8.7% intergenic (+37/‑112) 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,8900GA 8.7% 429.3 / 25.0 138intergenic (+37/‑112)garD/prlFD‑galactarate dehydrogenase/antitoxin of the SohA(PrlF)‑YhaV toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (66/60);  new base A (5/7);  total (71/67)
Fisher's exact test for biased strand distribution p-value = 5.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

Only 100 of 138 total aligned reads displayed.

AAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACT  >  NC_000913/3276743‑3277033
                                                                                                                                                   |                                                                                                                                               
aaGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGt                                                                                                                                                <  1:157391/149‑1 (MQ=255)
  gttgttCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGt                                                                                                                                                <  1:862503/147‑1 (MQ=255)
       tCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGt                                                                                                                                         >  4:1436233/1‑149 (MQ=255)
       tCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGt                                                                                                                                         >  2:812237/1‑149 (MQ=255)
         cACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCACTAAGGCCATGttt                                                                                                                                       <  2:695222/149‑1 (MQ=255)
            tttATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCACTAAGGCCATGt                                                                                                                                         <  2:247826/144‑1 (MQ=255)
            tttATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCACTAAGGCCATGt                                                                                                                                         >  1:247826/1‑144 (MQ=255)
                   tcGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGttt                                                                                                                                       >  2:330510/1‑139 (MQ=255)
                   tcGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTa                                                                                                                             >  3:856819/1‑149 (MQ=255)
                   tcGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTa                                                                                                                             >  4:327891/1‑149 (MQ=255)
                     gACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTaaa                                                                                                                           <  3:1436233/149‑1 (MQ=255)
                     gACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTaaa                                                                                                                           >  1:611445/1‑149 (MQ=255)
                              aGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACg                                                                                                                  >  2:721025/1‑149 (MQ=255)
                               gCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACGTGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGt                                                                                                                 >  4:1547064/1‑149 (MQ=255)
                                  ggAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGt                                                                                                                                <  3:1579689/131‑1 (MQ=255)
                                    aagaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCa                                                                                                                                            <  1:481218/117‑1 (MQ=255)
                                    aagaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCa                                                                                                                                            >  2:481218/1‑117 (MQ=255)
                                         gaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGCCCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCACTAAGGCCATGtt                                                                                                                                        <  2:1787571/116‑1 (MQ=255)
                                         gaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGc                                                                                                                                              <  2:1216667/110‑1 (MQ=255)
                                         gaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaataa                                                                                                       >  1:1355459/1‑149 (MQ=255)
                                         gaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCACTAAGGCCATGtt                                                                                                                                        >  1:1787571/1‑116 (MQ=255)
                                          aaaaCCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTaataat                                                                                                      >  3:1195136/1‑149 (MQ=255)
                                            aaCCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTaaa                                                                                                                           >  4:1325835/1‑126 (MQ=255)
                                              ccTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCAc                                                                                                                   <  4:1174198/132‑1 (MQ=255)
                                                ttCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcattt                                                                                                >  3:1891664/1‑149 (MQ=255)
                                                   tcGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTATTTAAAGcccc                                                                                                                      <  3:1084618/124‑1 (MQ=255)
                                                   tcGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCACTAAGGCCATGtt                                                                                                                                        >  3:941120/1‑106 (MQ=255)
                                                     ggATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatttgcatt                                                                                           <  1:334084/149‑1 (MQ=255)
                                                     ggATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatttgcatt                                                                                           >  3:434849/1‑149 (MQ=255)
                                                      gATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTa                                                                                          >  3:688366/1‑149 (MQ=255)
                                                       aTCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTAc                                                                                         <  2:611445/149‑1 (MQ=255)
                                                       aTCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTAc                                                                                         <  3:327891/149‑1 (MQ=255)
                                                        tCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTa                                                                                          <  2:1377503/147‑1 (MQ=255)
                                                        tCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTa                                                                                          >  1:1377503/1‑147 (MQ=255)
                                                        tCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTAcc                                                                                        <  4:856819/149‑1 (MQ=255)
                                                        tCAATGGGGGCAGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTAcc                                                                                        >  3:1111915/1‑149 (MQ=255)
                                                          aaTGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAGTAATGCAAATGCATTattaa                                                                                      <  1:721025/149‑5 (MQ=255)
                                                           aTGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTa                                                                                     <  2:1355459/149‑1 (MQ=255)
                                                                       tAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatttgc                                                                                              >  4:613798/1‑128 (MQ=255)
                                                                       tAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatttgc                                                                                              >  1:68473/1‑128 (MQ=255)
                                                                       tAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatttgc                                                                                              <  2:68473/128‑1 (MQ=255)
                                                                       tAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatttgc                                                                                              <  3:613798/128‑1 (MQ=255)
                                                                         aCCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCAt                                                                                                                                           >  3:1241895/1‑81 (MQ=255)
                                                                         aCCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCAt                                                                                                                                           <  4:1241895/81‑1 (MQ=255)
                                                                         aCCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATaa                                                                       >  4:1558284/1‑149 (MQ=255)
                                                                          ccAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAAt                                                                      >  1:1023654/1‑149 (MQ=255)
                                                                           cAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAAtt                                                                     >  1:1551724/1‑149 (MQ=255)
                                                                              cTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGa                                                                  <  1:527423/149‑1 (MQ=255)
                                                                                ggCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGc                                                                <  4:688366/149‑1 (MQ=255)
                                                                                 gCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCACACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCa                                                               <  2:1023654/149‑1 (MQ=255)
                                                                                       tttAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAAc                                                         >  1:838896/1‑149 (MQ=255)
                                                                                       tttAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAAc                                                         >  1:1453369/1‑149 (MQ=255)
                                                                                        ttAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACa                                                        >  1:784839/1‑149 (MQ=255)
                                                                                          aaCCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAgg                                                      <  2:1551724/149‑1 (MQ=255)
                                                                                                  aCCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCCTGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATGTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTaaaa                                              >  1:1213254/1‑149 (MQ=255)
                                                                                                  aCCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTaaaa                                              <  4:434849/149‑1 (MQ=255)
                                                                                                      gTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGAc                                          <  3:1558284/149‑1 (MQ=255)
                                                                                                      gTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGAc                                          <  1:1688987/149‑1 (MQ=255)
                                                                                                       tGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACa                                         <  2:1621092/149‑1 (MQ=255)
                                                                                                           cTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGa                                           >  4:1362924/1‑143 (MQ=255)
                                                                                                           cTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGa                                           <  3:1362924/143‑1 (MQ=255)
                                                                                                           cTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCACTAAGGCCATGt                                                                                                                                         >  1:366048/1‑49 (MQ=18)
                                                                                                            tGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCa                                                                                >  3:1473824/1‑105 (MQ=255)
                                                                                                            tGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGa                                                                  <  1:1019243/119‑1 (MQ=255)
                                                                                                            tGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGa                                                                  >  2:1019243/1‑119 (MQ=255)
                                                                                                            tGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAgg                                                      <  3:1610637/131‑1 (MQ=255)
                                                                                                            tGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAgg                                                      >  4:1610637/1‑131 (MQ=255)
                                                                                                            tGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGaa                                    >  2:718421/1‑149 (MQ=255)
                                                                                                              aTTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAg                                                       >  3:322555/1‑128 (MQ=255)
                                                                                                              aTTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAg                                                       <  4:322555/128‑1 (MQ=255)
                                                                                                              aTTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATc                                  <  2:838896/149‑1 (MQ=255)
                                                                                                              aTTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATc                                  >  1:1689195/1‑149 (MQ=255)
                                                                                                              aTTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATc                                  >  2:405574/1‑149 (MQ=255)
                                                                                                               tttCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGaa                                    <  2:1500502/146‑1 (MQ=255)
                                                                                                                ttCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGACATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCAt                                <  2:1689195/149‑1 (MQ=255)
                                                                                                                   ttttGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATg                               >  4:1735999/1‑147 (MQ=255)
                                                                                                                   ttttGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGc                              <  4:867015/148‑1 (MQ=255)
                                                                                                                   ttttGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGc                              >  3:867015/1‑148 (MQ=255)
                                                                                                                    tttGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaataa                                                                                                       >  3:1781261/1‑74 (MQ=255)
                                                                                                                    tttGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaataa                                                                                                       <  4:1781261/74‑1 (MQ=255)
                                                                                                                    tttGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGccc                            >  4:1310790/1‑149 (MQ=255)
                                                                                                                      tGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaat                                                                                                         <  4:450090/70‑1 (MQ=255)
                                                                                                                      tGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaat                                                                                                         >  3:450090/1‑70 (MQ=255)
                                                                                                                         aaCGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGa                                                                  <  4:1461386/106‑1 (MQ=255)
                                                                                                                         aaCGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGa                                                                  >  3:1461386/1‑106 (MQ=255)
                                                                                                                               ccTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTta                                                                          <  2:384204/92‑1 (MQ=255)
                                                                                                                                cTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGaaaa                                                          <  3:1326592/107‑1 (MQ=255)
                                                                                                                                cTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGaaa                                                           >  2:868614/1‑106 (MQ=255)
                                                                                                                                cTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGaaa                                                           <  1:868614/106‑1 (MQ=255)
                                                                                                                                 tCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAg                                                                              <  3:478461/86‑1 (MQ=255)
                                                                                                                                     aaaaCATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCa           <  4:244135/149‑1 (MQ=255)
                                                                                                                                     aaaaCATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCa           >  3:1276333/1‑149 (MQ=255)
                                                                                                                                     aaaaCATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCa           >  1:1843896/1‑149 (MQ=255)
                                                                                                                                     aaaaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCtctc            <  2:29676/148‑1 (MQ=255)
                                                                                                                                     aaaaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCtctc            >  1:29676/1‑148 (MQ=255)
                                                                                                                                      aaaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTta                                                                          <  1:570418/85‑1 (MQ=38)
                                                                                                                                       aaCATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCAt                                <  1:691162/126‑1 (MQ=255)
                                                                                                                                            ggCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTa    <  4:1111915/149‑1 (MQ=255)
                                                                                                                                              ccTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACt  >  3:362142/1‑149 (MQ=255)
                                                                                                                                               cTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATaa                                                                       <  4:608035/79‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
AAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACT  >  NC_000913/3276743‑3277033

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