Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 361,670 A→T 18.9% V64E (GTA→GAA)  lacA ← thiogalactoside acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913361,6700AT18.9% 115.5 / 5.9 70V64E (GTA→GAA) lacAthiogalactoside acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (12/41);  new base T (10/4);  total (22/48)
Fisher's exact test for biased strand distribution p-value = 1.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTC  >  NC_000913/361439‑361889
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                   
ccGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTTccccaga                                                                                                                                                                                                                              >  1:118795/1‑235 (MQ=255)
ccGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGCTGGGCGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAAACGAAAAAGAAACAAACAGGCGGTTTTATCCacac                                                                                                                                                                                                                              >  1:379530/1‑236 (MQ=255)
          tgTTAATGCCAATCGTTATCGGAAAAGAGTTCATCTCGGCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTccc                                                                                                                                                                                                                     <  1:975179/236‑1 (MQ=255)
              aTTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATTGATGTTGGAACCGTAAGAGAAATAGAAAGGCGGATCTCACCAAGGGGTTTTtccacacc                                                                                                                                                                                                                >  2:894396/1‑234 (MQ=255)
                     aTCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCTTCCTATATGGCTTTTGGAACCTTAAGTGAAATAGCCATGCGGGGTTTCCCCGGCTTTTTccccccacccgggaaa                                                                                                                                                                                                        >  2:500311/1‑227 (MQ=255)
                       tgTTATCGGAAAAGGGTACATCTCGCCGTTTTTGCTCAATTCATGGTGTACAGTGTGTCCCGTAACGGAAAGAGTAAAGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACa                                                                                                                                                                                                   <  2:853359/238‑1 (MQ=255)
                                       tACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGCCGGGGCTTCCCCGCGGTGTTTCCCTCCCCTGGGAAAACATTTTTTAATAAggta                                                                                                                                                                                   >  2:691349/1‑237 (MQ=255)
                                                                     tgtACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATTAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGa                                                                                                                                                   <  1:507610/239‑1 (MQ=255)
                                                                          aGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCATTACGTTATCGCCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGAt                                                                                                                                              <  1:919086/239‑1 (MQ=255)
                                                                          aGGGTGTCCCGTAAAGGAAAGGAGAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGGTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGAt                                                                                                                                              <  1:408838/239‑1 (MQ=255)
                                                                            ggTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTTCCCCCGCGGTTTCCCCCCCCGCGGCAAAAAATTCTTTTATCCAGCATTTTCTTTTTTTAACTTCTGTTGGATGGGCGTGGTTa                                                                                                                                            >  2:931228/1‑239 (MQ=255)
                                                                               gtCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTATTTTCCCAGCGCTTTTCCCCTCCGTGGCCCACATTTCTTTTATATGACGTTCTCTTTTTTTAACtttttttcgtagcggggagtataac                                                                                                                                           >  2:367870/1‑216 (MQ=255)
                                                                               gtCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTA‑‑AGAAATAGACAGGCGGGTCTACCCAGGCGGTTTTCCCCTCCGTGGCAAAACATTTTTTTAATAGGATTTCTCTATTTTCAACTTTCGTTGGATTGGAGttgatatact                                                                                                                                        >  1:922358/1‑229 (MQ=255)
                                                                                 cccGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTTTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCCACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTc                                                                                                                                       <  1:847511/239‑1 (MQ=255)
                                                                                 cccGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGCGAAATAGACAGGCGGTTTTACCCACACGCTTTCCCCCCCCTCGGCAAACAATTCTTTTTATAGGATTTTTCTTTTTTTTACTTTTgtgtgttaggggtggtaaaaaca                                                                                                                                         >  2:630355/1‑217 (MQ=255)
                                                                                 cccGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTc                                                                                                                                       <  2:894630/239‑1 (MQ=255)
                                                                                               agTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTTCCCCAGCGGTTTTCCCCACCCCGGCAAACATTTCTTTTTTTAGGCTTTTTCTTTTTTTTCATTCTCTTGGATGAGAGTGATTAAACTCCTAAATTaaatttt                                                                                                                        >  2:746293/1‑234 (MQ=255)
                                                                                                 tAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATATACAGGCGGGTTTACCCCCGCGGTTTCCCCTTCCCTGGGCAAAAATTTTTTTATTAGGCGTTTTCTTTTTTTTACTTCTGTTGGTTTGTAGTGTTTTAAATCATCAAATAAAGtttt                                                                                                                      >  2:782579/1‑238 (MQ=255)
                                                                                                       ttggTGCAATCTGTTCTTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGa                                                                                                               <  2:385819/237‑1 (MQ=255)
                                                                                                         ggTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGGTTTCCCCTTCCGTGGCAAACATATCTTTAATCAGGCTTTCTCTTTTTTTAAATTCTGGTGGCTGCGAGGTGTTAAACACATACCTTTAAGTTTTTCCCCgcaa                                                                                                             >  1:842942/1‑236 (MQ=255)
                                                                                                                tCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTTTAACCAAGCGTTTTTCCCCTCCGGGGCAAACATTTTTTTAATTAGGCTTTTCTTTTTTTTAACTTCTGATGGTATGCGGGGGTTTAAACCCATAATTAACCATTTTCCCCAAAGGtttttt                                                                                                      >  1:299954/1‑239 (MQ=255)
                                                                                                                     aCGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTTTAACCAGGCGGTTTTCCCCAACCGGGGAAACAATTCTTTTATCAGGCTTTCTTTTTTTTTACCTTCTTATTGGTGCGCGGCGTTAAAATCAACAATTAAAGTTTTTCCACCAAGACTTTTTTTCgg                                                                                                 >  1:550415/1‑239 (MQ=255)
                                                                                                                       gTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTa                                                                                               <  1:88016/239‑1 (MQ=255)
                                                                                                                         tATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTTCCCCAGCGTTTTTCCCCCCCCTGGCCAAAATTTCTTTATTACAGCTTTTTCTTTTTTTTACTTCTTTTGGATTGGCGTGGTTTAAATCAAACATTTAATTTTTTCCCCCAAATCttttttttcttaa                                                                                              >  2:88016/1‑231 (MQ=255)
                                                                                                                              ccGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTTTCCAAGCGTTTTTCCCCACCGCGGCAAACATTTCTTTTATTCAGCTTTCCTTTTTTTTCACTTTTTATTGATTCGCGTGATTAAAATCCTACATTTAACGTTTCCCCCCCAGTCCTTTTTTCGGTAAACCtt                                                                                         >  1:358087/1‑238 (MQ=255)
                                                                                                                              ccGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTTCCCCGGCGTGTTTCCCTACCGTGGCAAAAATTTCTTTTATCACGCTTTCTCTTTTTTTCACTTCTTCTGTAGTCGTGTGTTTAAAATCATACACTTAACTTTTTCCACCAACGCTTTTTTTCGGTTAACCTTc                                                                                        >  2:666559/1‑239 (MQ=255)
                                                                                                                              ccGATTGTTACCGTGTAGTCATAGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGATCTACCCTCGCGGTTTCCCCCTCCCGCGCAAAACTTTATTTTTATAGGCCTTTTCTTTTTTTCACTATTGATGGACGCAGGGGGTTATACTCAAACATTTACTTTTTTTCCCCACAGTCTTTTTTcgttaaacatt                                                                                          >  2:504585/1‑228 (MQ=255)
                                                                                                                                 tttttttccgtgttggtaacGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGca                                                                                     <  2:550415/220‑1 (MQ=255)
                                                                                                                                 atggTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCCGGAGGTTTCCCCCACCCGTGCCAAAAATTCTTTTATCAGGCTTTCTCTTTTTTTCACTTCTGAATGGATGAGGGGATTTAAATCCTACATTTACCGTTTTCCCCCCAAGTCTTTTTTCGGGAAACCTTTcc                                                                                       >  1:579098/4‑237 (MQ=255)
                                                                                                                                 aTTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAAGCGTTTTTCCCTACCCGGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTTAACTTCTGTTTGTGGCCACGGATTAAACTCCTAAGTTAACATTTTTCCCCCAAGGCTTTTTTTCCGTTAGCCTTCgac                                                                                     >  1:428813/1‑237 (MQ=255)
                                                                                                                                   tttttccGTGTAGTCAATGACAATTGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTTCCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTGCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGcaca                                                                                   <  2:505559/234‑1 (MQ=255)
                                                                                                                                   tGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGcaca                                                                                   <  1:746293/239‑1 (MQ=255)
                                                                                                                                   tGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGcaca                                                                                   <  1:329366/239‑1 (MQ=255)
                                                                                                                                      tACCGTGGAGTCATCGACAATGGTTAAATTGAAATTGGCATAAAAATGGCGGCCTATATGGATGGTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACatat                                                                                <  1:500311/239‑1 (MQ=255)
                                                                                                                                        ccGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGcaca                                                                                   <  1:782575/234‑1 (MQ=255)
                                                                                                                                                      aaaaTGGTTAAATTGAAATTTGCATAAAAATTGGGGCCTTTATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACCCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGcc                                                                <  2:413753/237‑1 (MQ=255)
                                                                                                                                                       gattGGTTAAATTGAAATTTGCAAAAAAATTGCGGCGTATAAGGATGTTGGAACCGTAAGGGAAATAGTCAGGCGGGTGTACCCTGGCGTTTGCGTCTTCCTGGGCAAAAATTTCGTTAATCAGGCTTTCTCTTTTTTCAACGTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTACCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCTCATATCGGTAAATAGCTTGcc                                                                <  2:507874/235‑1 (MQ=255)
                                                                                                                                                       cAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACtt                                                                                                                                                                  <  1:210597/144‑1 (MQ=255)
                                                                                                                                                       cAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGcc                                                                <  1:894396/238‑1 (MQ=255)
                                                                                                                                                       cAATGGTTAAATTGAAATTGGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGcc                                                                <  1:1083305/238‑1 (MQ=255)
                                                                                                                                                           tgTTAAATTGAAATTTGCATAAAAATTGCGGCCTTTATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTTCCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGctc                                                           <  1:596445/238‑1 (MQ=255)
                                                                                                                                                                    gAAATTTGCATAAAAATTGCGGCCTATAGCGATGTTGGAACCGTAAGAGAAAAAGTATGGCGGTTCTACCCATGCGTTTTCCCCTACCGTTGCAAACATGTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCttt                                                  <  1:504585/239‑1 (MQ=255)
                                                                                                                                                                        attGAATAAAAATTTGGGACTATAGGGATGTTGGAACCGTA‑‑AGAAATGTACAGGCGGTTCTCCCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATTGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                            <  2:922358/238‑1 (MQ=255)
                                                                                                                                                                          ttcaaagaaaTTGGGGCCTATATGGATGTTTGACCCGTAAGAGAAATAGACAGGCGGTGCTACCCAGGCGTTTTCCCCTTCCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGGCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                            <  2:118795/232‑1 (MQ=255)
                                                                                                                                                                          tgactAAAAATTGCGGCCTATTTGGTTGTTGGAACCGTAAGAGAAATAGACCGGCGGTTCTCCCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                            <  2:75292/235‑1 (MQ=255)
                                                                                                                                                                          tGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                            <  1:533688/239‑1 (MQ=255)
                                                                                                                                                                               aaaaaTTGGGGGCTTTATTGATGTTGGACCCCGAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGc                                       <  2:569451/239‑1 (MQ=255)
                                                                                                                                                                               aaaaaTTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGc                                       <  1:931228/239‑1 (MQ=255)
                                                                                                                                                                                  aaTTGCGGCCTATTTGGATGTTTGAACCGTAAGAGAAATAGACAGGCTGTTTTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATg                                    <  2:229719/239‑1 (MQ=255)
                                                                                                                                                                                         gCCTATATGGATGTTGGAACCGTAAGAGAAAAAGTCAGGCGGTTTTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAAt                              <  2:381648/238‑1 (MQ=255)
                                                                                                                                                                                             atgtGGATTTTGGAACCGTAAGGGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGAt                         <  1:560783/236‑1 (MQ=255)
                                                                                                                                                                                             atatGGATGTTTGAACCGTAAGAGAAATAGACAGGCGGTTCTCACCAGGCGTGTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATTGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTGGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGAt                         <  1:409652/239‑1 (MQ=255)
                                                                                                                                                                                              tatGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTAcc                                                                                                                                                                                                               <  2:415234/62‑1 (MQ=255)
                                                                                                                                                                                              tatGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTc                                                 <  1:861215/214‑1 (MQ=255)
                                                                                                                                                                                               tgggttttttggaCCGTAAGGGAAATAGCCAGGCTGTTTTTACCATGCGTTTTCCCCGACCGTGGCAAACATTTCTTTAATCAGGCTTTCTTTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCa                       <  2:379530/227‑1 (MQ=255)
                                                                                                                                                                                               tggggtttttgAACCGTAAGAGAAATAGACAGGCGGTTCTTCCCTGGCTTTTTCCCCTACCGTGGCAAACATTTCTTTAATCCGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCa                       <  2:985772/229‑1 (MQ=255)
                                                                                                                                                                                               atGGATGTTGGAACCGTAAGAGAAATAGACCGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTTATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCa                       <  1:813402/239‑1 (MQ=255)
                                                                                                                                                                                                 ggATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACt                     <  1:889057/239‑1 (MQ=255)
                                                                                                                                                                                                    ttttGGAACCTTAAGAGAAATAGTAAGGCGGTTCTACCCAGGCGTTTTCACCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACGTCTGATTGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCTGTAAGCCTTCGCTCATATCGGTAAATAGCTTGTCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTcc                   <  2:470274/236‑1 (MQ=255)
                                                                                                                                                                                                    ttttGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                  <  2:579098/237‑1 (MQ=255)
                                                                                                                                                                                                     ttttgAACCGTTAGGGAAATATACCGGTGGTTTTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                  <  2:474976/234‑1 (MQ=255)
                                                                                                                                                                                                     gTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                  <  1:674025/238‑1 (MQ=255)
                                                                                                                                                                                                       tGGAACCGTAAGGAAAATAGACAGGCGGTTCCACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTa               <  2:435365/239‑1 (MQ=255)
                                                                                                                                                                                                       tGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGtt                <  1:933058/238‑1 (MQ=255)
                                                                                                                                                                                                        ggAACCGTAAGAGAAATAGACAGGCGGTTTTTCCCAGGCGGTTTTCCCTACCCTGGCcaaaa                                                                                                                                                                                                   >  1:415234/1‑60 (MQ=255)
                                                                                                                                                                                                        ggAACCGTAAGAGAAATAGACAGGCGGTTTTACCCAGGCGTTTTCCCCTACCGGGGCAAACATTTCTTTTATCAGGCTTTCTCTTTTTTTAACTTCTGATGGATGCGAGGGATTAAACTCATACATTAACGTTTTTCCCCGAAGTACTTTTTTCGGTAAGCCTTTCCACATATCGGTAAATAGCTTTCCTGCTCTTTTTTTTTCTGTCCTTTGCATTTTCCATGCGAGCACTCCGtttt              >  1:591705/1‑237 (MQ=255)
                                                                                                                                                                                                            ccGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGata          <  1:64759/239‑1 (MQ=255)
                                                                                                                                                                                                               tAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGt       <  1:25054/239‑1 (MQ=255)
                                                                                                                                                                                                                     aaaTAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACTTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGca                                                                                     <  1:792572/155‑1 (MQ=255)
                                                                                                                                                                                                                      aaTAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTTGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTc  <  1:371709/237‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                   
CCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTC  >  NC_000913/361439‑361889

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