Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 361,668 C→T 18.6% E65K (GAA→AAA)  lacA ← thiogalactoside acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913361,6680CT18.6% 145.2 / 6.6 71E65K (GAA→AAA) lacAthiogalactoside acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (10/46);  new base T (10/4);  total (21/50)
Fisher's exact test for biased strand distribution p-value = 2.38e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCG  >  NC_000913/361434‑361905
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                            
gACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTAcc                                                                                                                                                                                                                                                          <  2:16920/239‑1 (MQ=255)
      cgaatcaggttgtttttgcAATCGTTAGCGGAAAAGGATTCCTCCCGGCGTTTTTTCTAAATTCATGGTGTACAGGGTGTCCAGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCg                                                                                                                                                                                                                                                   <  1:170636/221‑1 (MQ=255)
      cGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCg                                                                                                                                                                                                                                                   <  2:77097/239‑1 (MQ=255)
         tCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGtt                                                                                                                                                                                                                                                <  2:156085/238‑1 (MQ=255)
                  tATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGGTTTTCCCAAGGCGTTTTCCCCtc                                                                                                                                                                                                                                        >  2:153072/1‑237 (MQ=255)
                          aTCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTTGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTACATGGTTGTTGGCACCGCAAGAAAAATAGACAGGCGGGTATACACCAGCGTTTTTCCCCCCCCGGGCaaa                                                                                                                                                                                                                           >  1:33304/1‑239 (MQ=255)
                            ctttttCGGAAAAGAGTACCTCTCGCCGTTTTTTCTCAATTCATGGTGTACCGGGTTGCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCCCCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACa                                                                                                                                                                                                                        <  1:104211/234‑1 (MQ=255)
                                                 ctcGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATCTGGATGTTGGAACCGGAAGAGAAATAGACAGGCGGTTCTTCCCAGGCGGTTTACCCTCCCCGGGGCAAACATTTTTTTAAtaggcgttt                                                                                                                                                                                                     >  1:131902/1‑229 (MQ=255)
                                                    gCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTAccc                                                                                                                                                                                                                                                         <  1:714974/188‑1 (MQ=255)
                                                    gCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTAccc                                                                                                                                                                                                                                                         >  2:714974/1‑188 (MQ=255)
                                                      cGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGAGCAATAAGTAAGTTATCACCGATTGTTACCGTGTAGTAACCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATCTGGATGTTTGAACCGTAAGAGAAATAAACAGGGGGTCTTTCCCCAGCGTTTTTTCCCACCCTGGGCAAACAttttttatacatggtttagtttttt                                                                                                                                                                                               >  1:583912/1‑220 (MQ=255)
                                                             tctcAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTTTACCCAGCGGGTTTTCCCTTCCGGGGGGAAAAATTTTTTATTTAGGCATTTTTCtttttttcact                                                                                                                                                                                     >  2:20891/1‑234 (MQ=255)
                                                                           gtACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGAAAGGCGGTTTTACCCCAGCCTTTTTCCCTACCCGCGCAAACAATTTTTTTAACTAGCttttttttttttttaccctttcttggtgggaa                                                                                                                                                                      >  1:147206/1‑219 (MQ=255)
                                                                               gggtggtCCGGTAAGGAAAAAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATTGACAATGGTTAAATTGAAATTTGCGTAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGTGGTTCGACACAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTATCAACTTCTGATGGATGCGAGTGAt                                                                                                                                                                 <  1:893868/233‑1 (MQ=255)
                                                                                    gtCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTTTACCCCGGCGGTTTTCCCTACCAGTGCCAACAATTTCTTAATTAGGCTTTTTTTTTTTTTAACCTCTtatggaaggcaaggtgtaaat                                                                                                                                                            >  1:590288/1‑224 (MQ=255)
                                                                                     tCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTTACCACGCGCTTTTCCCCTACCTGGGCAACAATTTTTTTCTATCGCTTTTCTTTTTTTTCAACtttttttgaatcgcagagataaaac                                                                                                                                                             >  1:353850/1‑216 (MQ=255)
                                                                                      cccGTAACGGAAATAGTATCGTTGGGGCCAATCATTTCCTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGTATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCATGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTc                                                                                                                                                         <  1:219832/239‑1 (MQ=255)
                                                                                      cccGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTTTTCCCCAGCGTTTTTCCCTCCCCGGGCCAACATTTTTTTTATATGGATTTTTCTTTTTTTAACTTTTgtttgcttcaggtgattcaaaca                                                                                                                                                           >  2:574882/1‑217 (MQ=255)
                                                                                      cccGTAACGGAAAGAATAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTc                                                                                                                                                         <  2:33304/239‑1 (MQ=255)
                                                                                      cccGTAACGGAAAAAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTc                                                                                                                                                         <  2:745885/239‑1 (MQ=255)
                                                                                       ccttAACGGAAAGGGGAACTTTGGGTTCAATCAGTACGTTAACACCGTTTTTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCa                                                                                                                                                        <  1:749242/236‑1 (MQ=255)
                                                                                       ccGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTTCGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCa                                                                                                                                                        <  2:88276/239‑1 (MQ=255)
                                                                                            aCGGAAAAAGTAACGTTGGGTTCAATTATTTCGTTAACACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACAt                                                                                                                                                  <  1:153072/239‑1 (MQ=255)
                                                                                                            tGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGa                                                                                                                                <  2:864563/239‑1 (MQ=255)
                                                                                                              ggTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTTTACCCCGCCGTTTTTCCCTACCCCGAGAAACACTTTTTTCATATCGATTTCTCTTTTTTTACATTTCGATGGGTGCAAGTGTTTAAAATCATaaaaaaacaatttcacccaaca                                                                                                                                 >  1:597693/1‑218 (MQ=255)
                                                                                                                          aCGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCgg                                                                                                                  <  2:47379/239‑1 (MQ=255)
                                                                                                                           cggtATCCCCGATTGTTACCGTGTAGTCATCGACAAGGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTa                                                                                                                <  2:353906/236‑1 (MQ=255)
                                                                                                                            gttacttgccttttttACCTTTTTGTCATCGACAATGGTTAAATTGAAATTTGCAAAATAATTGCGGCCTATATGGTTGTTTGAACCGGAAGAGAAATAGACAGGCGGTTCTTCCCAAGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCAAGAAGTCTTTTTTCCGGTa                                                                                                                <  1:439869/225‑1 (MQ=255)
                                                                                                                                 caacGATTGTTACCGTGTAGTCAGCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCt                                                                                                          <  2:578426/236‑1 (MQ=255)
                                                                                                                                  aCCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGt                                                                                                                 <  2:332606/232‑1 (MQ=255)
                                                                                                                                      ttttttcccGTTTTGTCTTCGACAATGTTTAAATTGAAATTTGCATAAAAATTGCGGCCTGTATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTTGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGca                                                                                                     <  1:336172/232‑1 (MQ=255)
                                                                                                                                      aTTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTTCCCTACCCGGGCAAACAATTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGGAGGATTAAAATCATAACTTAACGTTTTCCCCCGAAGTCTTTTTTTCGTTAAGCCTTctcc                                                                                                     >  2:784220/1‑236 (MQ=255)
                                                                                                                                       ttGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTTTACCCAGCCGGTTTTCCCACCCGTCGCAACAAATTTTTTAATATGGCGTTATCTTTTTTTAACTTCTTATGTGAGCGGTTGTTTATACTCATAAATTTAAATTTTCTCACGCTATTTTTTTTTCGGTAATCCTTCGcca                                                                                                    >  1:92975/1‑237 (MQ=255)
                                                                                                                                       ttGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGca                                                                                                     <  2:147206/238‑1 (MQ=255)
                                                                                                                                            aCCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCCCGCGTTTTTCCCCCCCGGGGGAAACAATTTTTTTTTCAGGCCTTCTTTTTTTTTAACTTTTTTTGGATTCGAGTGGTTTAACTACTCCCTTAAAGGTTTTCCACCAAAGCATTTTTTCCGGAAGccattcccccagt                                                                                                 >  1:332606/1‑227 (MQ=255)
                                                                                                                                              cGTGGAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCCACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCg                                                                                                                                                                       <  2:56744/170‑1 (MQ=255)
                                                                                                                                                tgtAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTTTACCCAAGCGGTTTCCCCTACCCTGGAAAACATTTCTTTTATTAGGCTTTCTTTTTTTTTCACCTCTTCTGGAATGCAGTGATTAAAATCATACATTAAAATTTTTCCCCAAAAACTTTTTTTTGTTAAACCTTTGccaaaaaatgta                                                                                           >  1:247777/1‑228 (MQ=255)
                                                                                                                                                tgtAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCa                                                                                                                                                        <  2:766561/182‑1 (MQ=255)
                                                                                                                                                     tCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACata                                                                                                 <  2:114240/228‑1 (MQ=255)
                                                                                                                                                      cATCGACAATGGTTAAATTTAAATTTGCATAAAAATTGCGGCCTATTTGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAg                                                                                      <  2:748504/238‑1 (MQ=255)
                                                                                                                                                          gACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTTCCCTACCCTGGCAAACATTTCTTTTATCAGGCTTTTTCTTTTTTCAACTTCTGATGGATGCGGGTGATTTAACTCCTACATTTACGTTTTTCCCCCAAGGCTTTTTTTCGGGAAGCCTTCCCACATATAGGTAAAAAGCTTg                                                                                  >  2:324576/1‑238 (MQ=255)
                                                                                                                                                           aCAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGcc                                                                                <  2:687329/239‑1 (MQ=255)
                                                                                                                                                            cAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGcc                                                                                <  2:278854/238‑1 (MQ=255)
                                                                                                                                                               tGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTTTACCCCGGCGGTTTTCCCCACCGTGGCAAACAATTTTTTTAATAGGCTTTTTCTTTTTTTAACTTCTGTTGGATGGGAGGGATTAAAACTCTACATTAACTGTTTTCCCCCAAAACTTTTTTTCGTTAAACCATTCCcacaaacgggaaaaaaagatctctta                                                                             >  1:114240/1‑217 (MQ=255)
                                                                                                                                                               tGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTTACCCGGCGGTTTTCCCCACCGTGGGAAAAATTTTTTTTATAAAGGCTTTTCTTTTTTTTAATTTCTTTTGGTGGTGGGGGTTAAAATCCAACAATTAACTTTTTCCCCCCAAGCTTTTTTTCTGGAAACCCTCCCCCAAATAGGTAAAAAACATGCCTg                                                                              >  1:720316/1‑237 (MQ=255)
                                                                                                                                                                        tGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCtt                                                                   <  2:247777/239‑1 (MQ=255)
                                                                                                                                                                               ttgaaaaaaaTTGGGGTCTATTGGGTTTTTGGAACCGTTAAGGAAAGAGACCGGCTGTTTTTCCCAGGCCTTTTCCCCTACCGGGGAAAACATTTCTTTAAGCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                                            <  1:125497/234‑1 (MQ=255)
                                                                                                                                                                               ttcaaaaaaaTTGCGGCCTTTATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                                            <  2:583912/234‑1 (MQ=255)
                                                                                                                                                                               ttcATAAAAATTGCGGCTTTTAGGGATTTTTGAACCGTAAGAGAAATAGACCGGCGGTTCTCCCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                                            <  1:784220/237‑1 (MQ=255)
                                                                                                                                                                               tGCATAAAAATTTCGGCCTATATGGATGTTGGAACCGTAAGAAAAATAGACAGGCGGTTCTTCCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                                            <  1:526733/239‑1 (MQ=255)
                                                                                                                                                                               tGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                                            <  2:250106/239‑1 (MQ=255)
                                                                                                                                                                               tGCAAATAAATTGTGGCCTTTATGGATGTTGGAACCGTAATAAAAATGGACAGGCGGTTCGACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                                            <  1:58808/239‑1 (MQ=255)
                                                                                                                                                                                   tAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTgg                                                        <  2:723700/239‑1 (MQ=255)
                                                                                                                                                                                             ggCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCaa                                               <  2:720316/238‑1 (MQ=255)
                                                                                                                                                                                               ccTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGTAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTa                                                                                           <  1:59803/192‑1 (MQ=255)
                                                                                                                                                                                                  tttttgATTTTGGAACCGTAAGAGAAATAGTCAGGCTGTTCTACCCATGCGTTTTCCCCTACCTGTGCAAACATTTCTTTAATCAGGTTTTCGCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGAt                                         <  1:308442/234‑1 (MQ=255)
                                                                                                                                                                                                      gggtGTTTGAACCCGTAGAGAAATAACCCGGGGGTTCTACCCATGCGTTTTCCCCTACCGTGGCAAACATTTCGTTAATCAGGCTTTCTCGTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACt                                     <  1:6285/236‑1 (MQ=255)
                                                                                                                                                                                                      ggATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTTCCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACt                                     <  2:353850/239‑1 (MQ=255)
                                                                                                                                                                                                      ggATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACt                                     <  2:92975/239‑1 (MQ=255)
                                                                                                                                                                                                      ggATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACt                                     <  2:590288/239‑1 (MQ=255)
                                                                                                                                                                                                         ttttGGAACCGTAAGAGAAATAGACAGGCGTTTCTTCCCAGGTGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGGGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                                  <  1:34791/237‑1 (MQ=255)
                                                                                                                                                                                                         ttttGAAACCGTAAGGAAAAATGAAAGGCGGTTTTCCCCAGGCGTTTTCCCCGACCGTGGCAAACATTTCTTTAAGCAGGCTTTCTCTTTTTTCAACTACTGATGGATGCGAGTGATTAAACTCATGCATTAACGTTTTCCCAAGAAGTCGTTTTTCCGGTAAGCCTTCGCACATAGCGGTAAATGGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTcc                                   <  1:20891/236‑1 (MQ=255)
                                                                                                                                                                                                         tGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCTTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                                  <  2:572170/239‑1 (MQ=255)
                                                                                                                                                                                                         tGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTcc                                   <  2:167977/238‑1 (MQ=255)
                                                                                                                                                                                                          gTTGGAACCGTAAGAGAAATAGACAGGCGGTTTTACCCAGGCTTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTTTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                                  <  2:346170/238‑1 (MQ=255)
                                                                                                                                                                                                           ttGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTTCCCCACCCTGGCAAAAATTTCTTTAATCAAGGTTTCTCTTTTTTTAACTTTTTATGGATGCGAGTGGTTAAACTCATACATTAAAGTTTTTCCCCGAAAGCTTTTTTTCGGTAAACCTTTGCCCATATCGGTAAAAAGCTTGCCTGCTCTTATTTTTTTGGTCATTGGCATGGTCAATGCTATCACTCCCtt                                >  2:119578/1‑239 (MQ=255)
                                                                                                                                                                                                             ggAACCGTAAGAGAAATAGACAGGCGGTTCTAACCAGGCGGTTTTCCCCACCCGGGGAAAAATTTTTTTTATTATGGTTTCTTTTTTTTTCACTTCTGATGGGAGCGGGGGATTAAACTCATTAATTAACGGTTTTCCCCCAAGTCTTTTTTTCCGGAACACTTCGCCCCCACTCTTAAAAAACTTGCCTTCTCTTATTCTTTTTTTCTTTGCGATGTTAAATgggctaaccccgcttt                              >  1:690866/1‑224 (MQ=255)
                                                                                                                                                                                                                    tAAGAGAAATAGACAGGCGGTTCTACCCCGGCGTTTTTCCCCACCCTGGGAAACATTTTTTTTATCAGGCTTTCTCTTTTTTTAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTTCCACGAAGTCTTTTTTTCCGTAAGCCTTCGCACCAATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGGTCAATGCGGTTACTCCCTTTTggttatt                       >  2:241763/1‑233 (MQ=255)
                                                                                                                                                                                                                          aaaTAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTc                                                 <  2:431272/207‑1 (MQ=255)
                                                                                                                                                                                                                          aaaTAAACAGGGTGTTTTCCCCAGGGGTTTTCCCCTACCGTGGGAAACATTTCTTTAAACAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGt                   <  1:574882/237‑1 (MQ=255)
                                                                                                                                                                                                                                          tCTAGCCAGGCGTTTTCCCCTTCCGTGGCAAACATTTCTTTAATCAGGTTTTCTCTTTTTTCAACTTCTGAGGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGGCTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATGCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTcgcg  <  2:805525/238‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                            
GACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCG  >  NC_000913/361434‑361905

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: