Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,174,6050TC45.7% 114.4 / 20.4 81G187G (GGA→GGG) ycfTinner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (31/13);  new base C (0/37);  total (31/50)
Fisher's exact test for biased strand distribution p-value = 3.83e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCG  >  NC_000913/1174455‑1174749
                                                                                                                                                        |                                                                                                                                                 
aaaCAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCAtt                                                                                                                                                   >  1:691970/1‑151 (MQ=255)
   caacCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTccc                                                                                                                                                >  1:881304/1‑151 (MQ=255)
      ccAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCccac                                                                                                                                             <  2:779539/151‑1 (MQ=255)
        aGGCTCCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCccacca                                                                                                                                           <  1:960993/151‑1 (MQ=255)
         ggCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCcaccac                                                                                                                                          <  1:133848/151‑1 (MQ=255)
            accgaccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACgg                                                                                                                                        >  2:445835/1‑150 (MQ=255)
              cgaccgCCAGAACGGTAAGCAAAGAAGCCATCTGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGt                                                                                                                                     <  1:1042168/151‑1 (MQ=255)
               gaccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTg                                                                                                                                    >  2:995418/1‑151 (MQ=255)
                accgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTgg                                                                                                                                   >  1:1249841/1‑151 (MQ=255)
                  cgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGa                                                                                                                                 <  1:240070/151‑1 (MQ=255)
                   gCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGaa                                                                                                                                <  1:974613/151‑1 (MQ=255)
                    ccAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAAc                                                                                                                               <  1:229732/151‑1 (MQ=255)
                      aGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACg                                                                                                                              >  2:486957/1‑150 (MQ=255)
                        aaCGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGaaa                                                                                                                           <  2:56445/151‑1 (MQ=255)
                          cGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAAtt                                                                                                                         >  2:904458/1‑151 (MQ=255)
                          cGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAAtt                                                                                                                         >  1:984980/1‑151 (MQ=255)
                             tAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AAc                                                                                                                     >  1:10328/1‑151 (MQ=255)
                               aGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCg                                                                                                                   <  2:1140122/151‑1 (MQ=255)
                               aGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCg                                                                                                                   <  2:1158932/151‑1 (MQ=255)
                                    agaagCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCcac                                                                                                               >  2:881757/1‑150 (MQ=255)
                                      aagCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACt                                                                                                            >  2:364195/1‑151 (MQ=255)
                                      aagCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACt                                                                                                            >  1:958386/1‑151 (MQ=255)
                                       agCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTc                                                                                                           >  2:849400/1‑151 (MQ=255)
                                              agcagATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGt                                                                                                     >  2:863303/1‑150 (MQ=255)
                                               gaggagggcggcGCAACGCCACCACTTAAACAAAGGACATTGTTGTCGCGAACAACCAGGCCCCAAGCATTGAAAAGGGTAATTT‑GC‑GGAACCCATGCCTCATCCCCCCCCAGCGCTTGGACATGAAATT‑AACCGCCACACTCAGCAGTa                                                                                                    <  1:263527/144‑1 (MQ=255)
                                                  ggtgggggaccagcGGCCCCTCTTAAAAACAGGTATTTATTGTCGCCCAAAGCCAGCCGCAAGGGTTAAAAAAAGGAAAATT‑GC‑GGACCACCCTTTTCATCCCCCCCCACGGGGGGGAACCAAAATT‑AACCGCCACACTCAGCAGTacaa                                                                                                 <  1:252836/137‑1 (MQ=255)
                                                     ggcggcGCAACGGCACCTCTTTAACACAGGTAATTTTTGTGGCGCAAAACCAGGCGCCAAGGCTGAAAAAAGGAAAATT‑GC‑GGATCACACTGTTCATCCCCCCCCCCGGCGGGGAAACAAAATA‑AACCGCCACACTCAGCAGTa                                                                                                    <  1:1051508/144‑1 (MQ=255)
                                                      gcggcgCAACGGCCCCTTTTTAACACAGGTCATTATTGTTGCGCCAACCCAGGCGCCAAGGCTGTAAAAAGGAAAATT‑GC‑GGATCACACTGTTCTTCCCCCCCCACGGCGGGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAg                                                                                            <  2:766187/151‑1 (MQ=255)
                                                        ggcgCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGc                                                                                          >  2:672663/1‑151 (MQ=255)
                                                        ggcgCAACGGCACCTCTTTAAAACAGTTTTTTATTGTTGCGCCAACCCAGGCGCAAAGGCTGTAATAGGGCAAAGT‑GC‑GGACCCCGCTTTTCATCCCCCCCCCCGGGGGGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGc                                                                                          <  1:1296390/151‑1 (MQ=255)
                                                             accGGCACCTTTTAAACACAGGTCATTATTGTTGCGCAAACCCATGCGCCAGGGCTGAAAAAAGGCAAATG‑GC‑GGACCCCACTTTTCATCCCCCCCCGCGGGGGGGGAACGAAATT‑AGCCGCCACACTCAGCAGTACAAACAAGGCGAAt                                                                                      <  1:713374/148‑1 (MQ=255)
                                                             aacgccccttctTAAACAAAGTTATTTTTTGTGGCGCAAACCCAGGGGCAAGGGTTGTAA‑AAGGAAAATTGGC‑GAATCCCATTGTTCATCCCCCCCCCCGGCGGGGAAAGGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCg                                                                                         <  1:1208573/138‑1 (MQ=255)
                                                                  cACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGc                                                                                >  2:17280/1‑151 (MQ=255)
                                                                   aCCTCTTTAACCCAGATCATTATTGCTACGCCGAACCATGTGCCAAGGCTGTACTAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGCACGACATT‑AACCGCCACACCCAGAAGTCCAAACAAGGCGAATAGTgcc                                                                                >  1:684243/1‑148 (MQ=255)
                                                                   aCCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGc                                                                                >  1:845208/1‑150 (MQ=255)
                                                                    ccTTTTTACCCCGGGTAATTATTGTGGCGCAAAACCAGGCCCAAGGGTTGAAAAAAGGAAAATG‑GC‑GGATCCCCCTTTCCATCCCCCCCCCCGGCGTGGAAACAAAATA‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCa                                                                               <  2:548135/150‑1 (MQ=255)
                                                                       ctTTAACAAGGGTCATTATTGTGGCCCAAAACCATGCGCCAAGGCTGAAATAGGGCAAATG‑GC‑GGATCCCACTTTTCTTCCCCCCCCCCGGCTAGGAAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAgc                                                                             <  2:887455/149‑1 (MQ=255)
                                                                         ttAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGcc                                                                          >  2:547242/1‑150 (MQ=255)
                                                                           aaCACGGGTATTTATTGTGGCGCCAAACCAGGCGCCAAGGCTGTAAAAAGGAAAATT‑GC‑GGATCACACTTTTCATCCCCCCCCACGGGGGGGAAAGAAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGcc                                                                          <  1:107211/148‑1 (MQ=255)
                                                                              acaGGTCTTTTTTTTGGCGCAAACCAAGGCGCAAGGGCGGTAAAAAGGAAAATG‑GC‑GGATCACCTTGTTCTTCCCCCCCCCCGGGGGGGAAACAAAATA‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACg                                                                     <  1:773376/150‑1 (MQ=255)
                                                                              aaaGGTAATTATTGTTGCGCAAAACCAGGCGCAAGGGTTAAAATAA‑GAAAATT‑GCGGGATCACACTGTTCATCCCCCCCCACGGGGGGGAAACAAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACgg                                                                    <  2:909479/149‑1 (MQ=255)
                                                                                   ttattattGTGGGTCAAAACCATCGGCCAGGGTTGAAAAAAAAAAAATT‑GC‑GGACCCAGCTGTACATCCCCCCCCCCGGCATGAAAAAAAAATA‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTa                                                                 <  2:468656/147‑1 (MQ=255)
                                                                                   tCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaa                                                                >  1:829496/1‑150 (MQ=255)
                                                                                     attattGTTGCGCAAACCCATGCCCCAAGGCTGAAAAAGGGAAAATG‑GC‑GGATCACCCTTTTCATCCCCCCCCACGGCGTGGAAACAAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaaa                                                              <  2:57764/150‑1 (MQ=255)
                                                                                      ttattGTGGCGAAAACCCAGCCGCCAAGCGTAAAAAAAGGCAAAGG‑GC‑GGAGCACACTTTGCATCCCCCCCGGCGGCGTGGAACCAAAATA‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGAc                                                                      <  1:704320/141‑1 (MQ=255)
                                                                                         ttgttgCGCCAACCCATGCCCCAGGTCTGAAAAAAGGAAAATT‑GC‑GGATCACACTGTTCATTCCCCCCCCCGGCGTGGAAACGAAATA‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAtt                                                           <  1:649558/149‑1 (MQ=255)
                                                                                         ttgttgCGCCAAACCATGCGCCAAGGCTGTAATAAGGAAAATT‑GC‑GGATCACACTGTTCATCCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAtt                                                           <  2:501770/149‑1 (MQ=255)
                                                                                         ttgttgCGCAAAACCATGCCCCAAGCTTGTAATAAGCCAAATT‑GC‑GGATCACCTTTTTCATCCCCCCCCCCGGGGTGGAAACAAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGc                                                                   <  2:804868/141‑1 (MQ=255)
                                                                                         ttgtggCGCAAACCCATGCGCAAGGGCTGTAAAAAGGCAAATT‑GC‑GGATCACACTGTTCATCCCCCCCCACGGCGGGGAAACAAAATA‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAttt                                                          <  1:392328/150‑1 (MQ=255)
                                                                                              gcgcAAAACCTTGCGCCAAGGCTTTAATAAGGAAAATT‑GG‑GGATCACTCTTTTCATTCCCCACCACGGGGGGGAAACAAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATa                                                    <  2:40238/151‑1 (MQ=255)
                                                                                               cgcCAACCCATGCGCCAAGGCTGAAAAAAGGCAAATT‑GC‑GGATCACCCTTTTCATCCCCCCCCACGGCGGGGGAACAAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATa                                                    <  2:1192204/150‑1 (MQ=255)
                                                                                               cgcCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAc                                                   >  1:1186443/1‑151 (MQ=255)
                                                                                                gcCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAcc                                                  >  1:570830/1‑151 (MQ=255)
                                                                                                 ccAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAcc                                                  >  2:1197015/1‑150 (MQ=255)
                                                                                                   aaaCCAGGCGCCAAGTCTGAAAAAAGGAAAATT‑GC‑GGACCACTCTGTTCATCCCCCCCCACGGCGGGGGAACAAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACg                                               <  1:204930/151‑1 (MQ=255)
                                                                                                    aaCCATGCCCCAAGGCTGTAAAAAGGCAAATG‑GC‑GGATCACACTGTTCATCCCCCCCCGCGGCGTGGGAACAAAATA‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCATGa                                              <  1:199531/151‑1 (MQ=255)
                                                                                                      caaagcgcAAGGGCTGAAAGAAGGAAATGG‑GC‑GGTTCCCACTGTTCTTTCCCCACCACGGGGTGGAAAAGAAATT‑AGCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGaaa                                            <  2:97993/147‑1 (MQ=255)
                                                                                                       cATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGCAAt                                           >  2:373259/1‑151 (MQ=255)
                                                                                                        aTGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAAt                                           >  2:980449/1‑150 (MQ=255)
                                                                                                        aTGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAAt                                           >  2:1218531/1‑150 (MQ=255)
                                                                                                             caaaGGCTGTAAAAAGGAAAATT‑GC‑GGATCCCATTTTTCATCCCCCCCCACGGCGTGGGAACAAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGa                                        <  2:208999/146‑1 (MQ=255)
                                                                                                              cAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAAt                                    >  2:16919/1‑151 (MQ=255)
                                                                                                                       aataaGGAAAATT‑GC‑GGATCACCCTGTTCATCCCCCCCCACGGGGTGGAAACAAAATA‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATa                           <  1:464676/151‑1 (MQ=255)
                                                                                                                       aaaaaGGAAAATT‑GG‑GGATCCCCCTTTTCATCCCCCCCCCCGGGGGGGAAACAAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCAt                            <  1:959135/147‑1 (MQ=255)
                                                                                                                       aaaaaGCCAATTT‑GC‑GGACCACATTGTTCATCCCCCCCCCGGGGGTGGGAAAAAAATA‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACg                               <  1:1110929/144‑1 (MQ=255)
                                                                                                                          aaGGAAAATT‑GC‑GGATCACACTTTTCATCCCCCCCCCGGGGGGGGAAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACa                         <  1:1241584/150‑1 (MQ=255)
                                                                                                                               aattt‑GG‑GGATCACACTGTTCATCCCCCCCCACGGGGTGGAAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATa                     <  2:511562/146‑1 (MQ=255)
                                                                                                                               aaaTT‑GG‑GGACCACGCTTTTCTTCCCCCCCCACGGGGTGGGAA‑GAAATTAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGa                       <  1:203028/147‑1 (MQ=255)
                                                                                                                               aaaTT‑GC‑GGACCCCCCTTTTCATCCCCCCCCCCGGGGTGGAAACAAAATA‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGAt                      <  1:1260043/148‑1 (MQ=255)
                                                                                                                               aaaTT‑GC‑GGACCACACTGTTCATCCCCCCCCACGGCGGGGAAACAAAATA‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAcc                   <  1:121694/151‑1 (MQ=255)
                                                                                                                               aaaTT‑GC‑GGAACACACTGTTCATCCCCCCCCACGGCGGGGGAAAAAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAc                    <  2:829148/150‑1 (MQ=255)
                                                                                                                               aaaTG‑GC‑GGATCCCTCTTTTCTTCCCCCCCCTCGGCGTGGAACCGAAATA‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAcc                   <  2:308653/151‑1 (MQ=255)
                                                                                                                                  tt‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACcac                 >  1:634909/1‑150 (MQ=255)
                                                                                                                                   t‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACcaca                >  1:531602/1‑150 (MQ=255)
                                                                                                                                     tggggATCACTTTTTTCTTCCCCCCCCCCGGCGGGGGAGAAAATTT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACca                  <  1:1051421/145‑1 (MQ=255)
                                                                                                                                     gC‑GGATCACCCTTTTCATCCCCCCCCACGGGGTGGGAACAAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACaa               <  2:144402/150‑1 (MQ=255)
                                                                                                                                      gcgcatccctttGTTCATCCCCCCCCCCGGCGTGGGAACAAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACaa               <  2:1249921/139‑1 (MQ=255)
                                                                                                                                         gATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTg           >  2:478596/1‑151 (MQ=255)
                                                                                                                                           tCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGt         >  1:1272592/1‑151 (MQ=255)
                                                                                                                                              cacTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCt      >  1:790599/1‑151 (MQ=255)
                                                                                                                                                  gTTCTTCCCCCCCCACGGCGTGGAAACAAAATA‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTgg    <  1:556210/149‑1 (MQ=255)
                                                                                                                                                   ttCATCCCCCCCCACGGGGGGGAAACAAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCg  <  1:457572/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                 
AAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATT‑AACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCG  >  NC_000913/1174455‑1174749

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