Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 361,668 C→T 19.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,6680CT19.6% 52.0 / 4.6 44E65K (GAA→AAA) lacAthiogalactoside acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (13/22);  new base T (7/2);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 5.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

ATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑CTACCCAGGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTC‑ATACATT‑AACG‑TTTTCCCACGAA‑GTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGG  >  NC_000913/361453‑361887
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                  
aTTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGAGGTTGGAACCGTAAGAGAAATAGAAAGGCGGG‑TT‑CTACCAAgcggt                                                                                                                                                                                                                         >  2:427954/1‑224 (MQ=255)
     cAATCGTTAACGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑CTACCCAGGC‑GTTTTccc                                                                                                                                                                                                                  <  2:528851/228‑1 (MQ=255)
                   aaaGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGCTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGG‑TT‑TTCCCCAGGG‑TGTTTccccacccgggaaaaa                                                                                                                                                                                                     >  1:93202/1‑216 (MQ=255)
                                 gCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TC‑CTACCCACAC‑GGTTTACCCCACCGTGGCAAAA‑AATTTTTTATTCAg                                                                                                                                                                                     >  1:195589/1‑228 (MQ=255)
                                 gCCGTTTTTTCTCAATTCATGGAGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGGGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAACATTGCGGCCTATCTGGATGTTGGAACCGTAAGACAAATAAACATAC‑GG‑GC‑CTACACTGGC‑GATTTTCCCGCCAGTTGGAAAC‑AAATCTTTTAttcg                                                                                                                                                                                     >  1:451238/1‑225 (MQ=255)
                                                 tCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑CTACCCAGGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTTc                                                                                                                                                                     <  2:359479/227‑1 (MQ=255)
                                                          aCAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑AT‑CTAACCCGGG‑GTTTTCCCCTACCGTGGGAAAC‑AATTCTTTAATCAGGC‑TTTTTCTTTTTTCAACTTCTgat                                                                                                                                                           >  2:476/1‑228 (MQ=255)
                                                                 tgcccGTAAAGGAAAGAGGAACGTTGGGTGCAATTAGTACGTTATCACCGATTGTTACCGTGGAGTCATTGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑CTACCCAGGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTTCAACTTCTGATGGATGCg                                                                                                                                                    <  1:427954/226‑1 (MQ=255)
                                                                 gtCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑TTACCCAAGC‑GTTTTTCCCTTCACTGGCAAACAATTTTTTTATACAGG‑‑TTTCtttttttttcattttttttggtggc                                                                                                                                                     >  2:50131/1‑210 (MQ=255)
                                                                 gtCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑CTACCCCGGC‑GTTTTTCCCTACCGTGGCAAAC‑AATTCTTTAATTAGGG‑ATTCTCTTTTTTTAACTTCTtatgggggtc                                                                                                                                                    >  2:96077/1‑223 (MQ=255)
                                                                   cccTTAACGGAAAGAGTAACGTTTGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑CTACCCAGGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTTCAACTTCTGATGGATGCGa                                                                                                                                                   <  1:606118/227‑1 (MQ=255)
                                                                   cccGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TTCCTACCAAGCC‑GTTTT‑TCCTACCCTGGGAAAC‑AATTTTTTTTTAAGACGTTTCTTTTTTTTCACCtgtccttattcacca                                                                                                                                                   >  1:378357/1‑214 (MQ=255)
                                                                    ccGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑CTACCCAGGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTtc                                                                                                                                                                              <  2:378357/199‑1 (MQ=255)
                                                                       tAACGGAAAGAGTAAAGTTGGGTGCAATCAGTTAGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑CTACCCAGGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGAt                                                                                                                                              <  1:198194/228‑1 (MQ=255)
                                                                          cGGAAAGAGTAACGGTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑CTACCCAGGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGAt                                                                                                                                              <  2:534052/225‑1 (MQ=255)
                                                                              aaGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGATAAATAGACAGGC‑GG‑TT‑CTACCCAGGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTc                                                                                                                                       <  2:174545/228‑1 (MQ=255)
                                                                                    agggTTGGGTTGAATTAGTACGTTAACCTCGGTTGTTACCGTGGAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑CTACCCAGGC‑GTTTTCCCCAACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTC‑ATACAt                                                                                                                                <  1:80172/225‑1 (MQ=255)
                                                                                      cGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑CTACCCCGGC‑GTTTTTCCCTACCGTGGCAAAC‑ATTTTTTTTATCTGAC‑TTTCTCTTTTTTCAACTTCTGATGGGTGTGAGTGATTATAATC‑TTACAtt                                                                                                                               >  2:19252/1‑227 (MQ=255)
                                                                                                 aTCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGAGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GGCTT‑TTCCCCGGGG‑GTTTTCCCC‑ACCCTGGGAAAC‑ATTTCTTTTAGTAGGC‑TTTCTCTTTTTTCCACTTCTGATGGGTTGGAGTGATTTAACTC‑ATACATTAAAAG‑TTTTcct                                                                                                                  >  2:80172/1‑227 (MQ=255)
                                                                                                     gTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑CTACCCAAGC‑GGTTTTCCCTACCGGGGCAAAC‑ATTTTCTTAAACAGGC‑ATTCTCTTTTTTTTACTTCTGATGGAGGGCAATGATTTAACTC‑ATACCTTAAACA‑TTTTcccccaa                                                                                                              >  2:570792/1‑224 (MQ=255)
                                                                                                           tATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑TTATCCAGGC‑GTTTTCCCCTACCACGGGAAAA‑AATTCTTTAATCAGGA‑CTTCTCTTTTTTTAACATCCGGTGGCTACGGTTGGTTAAACTCAAAAAATT‑AACG‑‑TTTCCCAccgcatttt                                                                                                        >  1:512835/1‑219 (MQ=255)
                                                                                                           tATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑CTACCCAGGC‑GTTTTTCCCTACCGTGGCCAAC‑ATTTCTTTAATTAAGC‑TTTCTCTTTTTTCAACTTCTGATGGTTGCGAGTGATTAAAATC‑CTACATT‑TACG‑TTTTTCCCCCAA‑GTCtt                                                                                                       >  2:545830/1‑227 (MQ=255)
                                                                                                            aTCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAACAATTGCGGCCTATATGGATGTTGGAACCGGAAGAGAAATAGACAGGC‑GG‑TT‑TTTACCAGAG‑GTTTTTCCCTACCACGGTGAAG‑ATTTTTTTAATACGGC‑ACTTTCTTTTTTCTACTTTTATTTGATTCAAGGTATTAAAACC‑ATACATT‑TACA‑CTTTTCCACgtaatgttt                                                                                                       >  2:574704/1‑219 (MQ=255)
                                                                                                                ccGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑GG‑CTTACCAGACGGTTTTCCCCAACC‑TGGGAAAA‑AATTTTTTTATCAAGA‑CTTTTCTTTTTTTACCTTTTCGTGGGAGTGTGTGCTTAAACTC‑ATACATT‑AACGTTTTTCCCACAAG‑Attttttttt                                                                                                  >  1:637974/1‑227 (MQ=255)
                                                                                                                               tAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑GT‑CTTCCCAGGA‑GTTTTTTCCTACCCTGGGAAAC‑ATTTCTTTTATCAAGC‑TTTTTCTTTTTTTAACTTCTTTTGTAGTCGAGCGATTTAAACC‑ATCAACT‑TACT‑ATTTTCCACCAA‑GACtttttttcccttacacattcaac                                                                                  >  1:587427/1‑211 (MQ=255)
                                                                                                                                            tGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑TTAACCAAGC‑GTTTTTCCCTACCCTGGGAAAC‑ATTTTTTTTATTAGGC‑TTTTCCTTTTTTTAACTTTTCATGGGGGCGACTGTGTTAACTC‑ATACATT‑TACC‑TTTTTCCCCTAAGGTCTTTTTT‑CTGAAAACACTCGCACATAACTGTAAAt                                                                      >  2:302557/1‑227 (MQ=255)
                                                                                                                                            tGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACAGTAAGAGAAATAGACAGGC‑GG‑TT‑TTACCCAGGC‑GGTTTCCCCTACCCAGGGAAAC‑ATTTTTTTAATCAAGG‑ATTTTCTTTTTTCACCTTCTGATGGGTGCGGGTGATTAAAACT‑ATACAAT‑AACA‑GTTTTCCACTAA‑GGCTTTTTTTCTGGAAACATTTGTACATATCTGTGAATa                                                                     >  1:91959/1‑228 (MQ=255)
                                                                                                                                                          tttGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑CTACCCAGGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTC‑ATACATT‑AACG‑TTTTCCCACGAA‑GTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGctct                                                        <  2:426883/227‑1 (MQ=255)
                                                                                                                                                                   aaaTTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑CTTACCAGGC‑GGTTTTCCCTATCCTGGCAAAA‑AATTCTTTAATACGGC‑TTTCTCTTTTTTTTACTTTTGTTGTATGCCAGTGGTTAAACTC‑ATACATT‑TAAG‑TTTTTCCCCGAA‑GTGTTTTTTCTGGTATACCTTCTCACATATCTGTAAATAGCTTTCCTGCTCCTTTTATTTct                                              >  2:651506/1‑227 (MQ=255)
                                                                                                                                                                       tgtgtcttattTGGAGTTTGGAACCCGAAAAGAAATATACAGGC‑TG‑TT‑CTACCCAGGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTC‑ATACATT‑AACG‑TTTTCCCACGAA‑GTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                          <  1:476/221‑1 (MQ=255)
                                                                                                                                                                            tctattggggtgtttgagacGTTAGGGAAAAAGACAGGC‑GT‑TT‑TTACCCTGTC‑GTTTTCCCCTTCCGTGGCAAAC‑CTTTCTTTAATCAGGC‑TTTCTCTTTGTTCAACTTCTGATGGATGCGAGTGATTAAACTC‑ATCCATT‑AACG‑TTTTCCCACGAA‑GTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGc                                     <  1:31183/209‑1 (MQ=255)
                                                                                                                                                                                          tgAACCTTAAGAGAAATAGACAGGC‑TG‑TT‑CTACCCAGGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTC‑ATACATT‑AACG‑TTTTCCCACGAA‑GTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGAt                       <  1:89881/227‑1 (MQ=255)
                                                                                                                                                                                          ggAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑CTACCCAGGCGGTTTTCCCCT‑CCGTGGGAAAC‑ATTTCTTTAATCAGGC‑GTTCTCTTTTTTCAACTTCTGATGGATGCGCGGGAGTAAAGTA‑CTAACTT‑AACG‑GTTTTTCACGAA‑AGTTTTTTTCTGTGAACACTTTCCCCCTAACTGTAAATAAATTGCCTGGTCTTGTACTTTCTGTCAtttgaaagtgtaatacgat                       >  1:178034/1‑211 (MQ=255)
                                                                                                                                                                                              ccGTAAGAGAAATAGACCGGC‑GG‑TT‑CTACCCAGGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTC‑ATACATT‑AACG‑TTTTCCCACGAA‑GTCTTTTTTCCGGTAAGCCTTCTCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACt                   <  2:639334/228‑1 (MQ=255)
                                                                                                                                                                                              ccGTAAGAGAAATAGACAGGC‑GG‑TT‑CTACCCAGGC‑GTTTTCCCCTACCGTTGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTC‑ATACATT‑AACG‑TTTTCCCACGAA‑GTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACt                   <  2:195589/228‑1 (MQ=255)
                                                                                                                                                                                              ccGTAAGAGAAATAGACAGGC‑GG‑TT‑CTACCCAGGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTC‑ATACATT‑AACG‑TTTTCCCACGAA‑GTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACt                   <  2:93202/228‑1 (MQ=255)
                                                                                                                                                                                              ccGTAAGAAAAATAGAAAGGC‑TG‑TT‑CTACCCAGGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTC‑ATACATT‑AACG‑TTTTCCCACGAA‑GTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGTTCATTGGCATGTTCAATGCGACCACt                   <  2:637974/228‑1 (MQ=255)
                                                                                                                                                                                               ccgtaagagaaaattgtgggggg‑TT‑CTACCCAGGC‑TTTTTCCCCTACCGTGGCAAAC‑ATTTCCTTAATCAGGC‑TTTCTATTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTC‑ATACATT‑AACG‑TTTTCCCACGAA‑GTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACt                   <  2:316380/207‑1 (MQ=255)
                                                                                                                                                                                                gTAAGAGAAATAGACCGGC‑GG‑TT‑CTACCCAGGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTC‑ATACATT‑AACG‑TTTTCCCACGAA‑GTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACt                   <  2:248454/226‑1 (MQ=255)
                                                                                                                                                                                                 taaggaaaaTAGAACGGC‑GG‑TT‑CTACCCAAGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTC‑ATACATT‑AACG‑TTTTCCCACGAA‑GTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                <  2:91959/222‑1 (MQ=255)
                                                                                                                                                                                                 tAAGAGAAATAGACCGGC‑GG‑TT‑TGACCCATGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTC‑ATACATT‑AACG‑TTTTCCCACGAA‑GTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                <  1:387081/228‑1 (MQ=255)
                                                                                                                                                                                                  aaggaaattAGACAGGC‑GG‑TT‑CTACCCCTGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTGCAACTTCTGATGGATGCGAGTGATTAAACTC‑ATACATT‑AACG‑TTTTCCCACGAA‑GTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                <  1:96077/219‑1 (MQ=255)
                                                                                                                                                                                                  aaGGGAAATAGACAGGC‑GG‑TT‑CTTCCCAGGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCATTTTTCAACTTCTGATGGATGCGAGTGATTAAACTC‑ATACATT‑AACG‑TTTTCCCACGAA‑GTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                <  2:512835/227‑1 (MQ=255)
                                                                                                                                                                                                                 gC‑GG‑TT‑CTAGCCAGGC‑TTTTTCCCCTACCGTGGCAAAA‑ATTTCTTTAAACAGGC‑TTTCTCTTTTTTCAACTTCTGATGGATGCTAGTGATTAAACTC‑ATTCATT‑AACG‑TTTTCCCACGAA‑GTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTgg  <  1:50131/226‑1 (MQ=255)
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                  
ATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TT‑CTACCCAGGC‑GTTTTCCCCTACCGTGGCAAAC‑ATTTCTTTAATCAGGC‑TTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTC‑ATACATT‑AACG‑TTTTCCCACGAA‑GTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGG  >  NC_000913/361453‑361887

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