Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,174,6100AC48.0% 111.0 / 19.0 77W186G (TGG→GGG) ycfTinner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (30/10);  new base C (0/37);  total (30/47)
Fisher's exact test for biased strand distribution p-value = 4.05e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCA  >  NC_000913/1174461‑1174756
                                                                                                                                                       |                                                                                                                                                   
ccAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCccac                                                                                                                                                    <  2:779539/151‑1 (MQ=255)
  aGGCTCCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCccacca                                                                                                                                                  <  1:960993/151‑1 (MQ=255)
   ggCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCcaccac                                                                                                                                                 <  1:133848/151‑1 (MQ=255)
      accgaccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACgg                                                                                                                                               >  2:445835/1‑150 (MQ=255)
        cgaccgCCAGAACGGTAAGCAAAGAAGCCATCTGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGt                                                                                                                                            <  1:1042168/151‑1 (MQ=255)
         gaccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTg                                                                                                                                           >  2:995418/1‑151 (MQ=255)
          accgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTgg                                                                                                                                          >  1:1249841/1‑151 (MQ=255)
            cgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGa                                                                                                                                        <  1:240070/151‑1 (MQ=255)
             gCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGaa                                                                                                                                       <  1:974613/151‑1 (MQ=255)
              ccAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAAc                                                                                                                                      <  1:229732/151‑1 (MQ=255)
                aGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACg                                                                                                                                     >  2:486957/1‑150 (MQ=255)
                  aaCGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGaaa                                                                                                                                  <  2:56445/151‑1 (MQ=255)
                    cGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAAtt                                                                                                                                >  2:904458/1‑151 (MQ=255)
                    cGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAAtt                                                                                                                                >  1:984980/1‑151 (MQ=255)
                       tAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAAc                                                                                                                            >  1:10328/1‑151 (MQ=255)
                         aGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCg                                                                                                                          <  2:1158932/151‑1 (MQ=255)
                         aGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCg                                                                                                                          <  2:1140122/151‑1 (MQ=255)
                              agaagCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCcac                                                                                                                      >  2:881757/1‑150 (MQ=255)
                                aagCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACt                                                                                                                   >  2:364195/1‑151 (MQ=255)
                                aagCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACt                                                                                                                   >  1:958386/1‑151 (MQ=255)
                                 agCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTc                                                                                                                  >  2:849400/1‑151 (MQ=255)
                                        agcagATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGt                                                                                                            >  2:863303/1‑150 (MQ=255)
                                         gaggagggcggcGCAACGCCACCACTTAAACAAAGGACATTGTTGTCGCGAACAACCAGGCCCCAAGCATTGAAAAGGGTAATTTGCGGAACCCATGCCTCATCCCCCCCCAGCGCTTGGACATGAAATTAACCGCCACACTCAGCAGTa                                                                                                           <  1:263527/144‑1 (MQ=255)
                                            ggtgggggaccagcGGCCCCTCTTAAAAACAGGTATTTATTGTCGCCCAAAGCCAGCCGCAAGGGTTAAAAAAAGGAAAATTGCGGACCACCCTTTTCATCCCCCCCCACGGGGGGGAACCAAAATTAACCGCCACACTCAGCAGTacaa                                                                                                        <  1:252836/137‑1 (MQ=255)
                                               ggcggcGCAACGGCACCTCTTTAACACAGGTAATTTTTGTGGCGCAAAACCAGGCGCCAAGGCTGAAAAAAGGAAAATTGCGGATCACACTGTTCATCCCCCCCCCCGGCGGGGAAACAAAATAAACCGCCACACTCAGCAGTa                                                                                                           <  1:1051508/144‑1 (MQ=255)
                                                gcggcgCAACGGCCCCTTTTTAACACAGGTCATTATTGTTGCGCCAACCCAGGCGCCAAGGCTGTAAAAAGGAAAATTGCGGATCACACTGTTCTTCCCCCCCCACGGCGGGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAg                                                                                                   <  2:766187/151‑1 (MQ=255)
                                                  ggcgCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGc                                                                                                 >  2:672663/1‑151 (MQ=255)
                                                  ggcgCAACGGCACCTCTTTAAAACAGTTTTTTATTGTTGCGCCAACCCAGGCGCAAAGGCTGTAATAGGGCAAAGTGCGGACCCCGCTTTTCATCCCCCCCCCCGGGGGGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGc                                                                                                 <  1:1296390/151‑1 (MQ=255)
                                                       accGGCACCTTTTAAACACAGGTCATTATTGTTGCGCAAACCCATGCGCCAGGGCTGAAAAAAGGCAAATGGCGGACCCCACTTTTCATCCCCCCCCGCGGGGGGGGAACGAAATTAGCCGCCACACTCAGCAGTACAAACAAGGCGAAt                                                                                             <  1:713374/148‑1 (MQ=255)
                                                       aacgccccttctTAAACAAAGTTATTTTTTGTGGCGCAAACCCAGGGGCAAGGGTTGTAAAAGGAAAATTGGCGAATCCCATTGTTCATCCCCCCCCCCGGCGGGGAAAGGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCg                                                                                                <  1:1208573/138‑1 (MQ=255)
                                                            cACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGc                                                                                       >  2:17280/1‑151 (MQ=255)
                                                             aCCTCTTTAACCCAGATCATTATTGCTACGCCGAACCATGTGCCAAGGCTGTACTAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGCACGACATTAACCGCCACACCCAGAAGTCCAAACAAGGCGAATAGTgcc                                                                                       >  1:684243/1‑148 (MQ=255)
                                                             aCCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGc                                                                                       >  1:845208/1‑150 (MQ=255)
                                                              ccTTTTTACCCCGGGTAATTATTGTGGCGCAAAACCAGGCCCAAGGGTTGAAAAAAGGAAAATGGCGGATCCCCCTTTCCATCCCCCCCCCCGGCGTGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCa                                                                                      <  2:548135/150‑1 (MQ=255)
                                                                 ctTTAACAAGGGTCATTATTGTGGCCCAAAACCATGCGCCAAGGCTGAAATAGGGCAAATGGCGGATCCCACTTTTCTTCCCCCCCCCCGGCTAGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAgc                                                                                    <  2:887455/149‑1 (MQ=255)
                                                                   ttAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGcc                                                                                 >  2:547242/1‑150 (MQ=255)
                                                                     aaCACGGGTATTTATTGTGGCGCCAAACCAGGCGCCAAGGCTGTAAAAAGGAAAATTGCGGATCACACTTTTCATCCCCCCCCACGGGGGGGAAAGAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGcc                                                                                 <  1:107211/148‑1 (MQ=255)
                                                                        acaGGTCTTTTTTTTGGCGCAAACCAAGGCGCAAGGGCGGTAAAAAGGAAAATGGCGGATCACCTTGTTCTTCCCCCCCCCCGGGGGGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACg                                                                            <  1:773376/150‑1 (MQ=255)
                                                                        aaaGGTAATTATTGTTGCGCAAAACCAGGCGCAAGGGTTAAAATAAGAAAATTGCGGGATCACACTGTTCATCCCCCCCCACGGGGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACgg                                                                           <  2:909479/149‑1 (MQ=255)
                                                                             ttattattGTGGGTCAAAACCATCGGCCAGGGTTGAAAAAAAAAAAATTGCGGACCCAGCTGTACATCCCCCCCCCCGGCATGAAAAAAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTa                                                                        <  2:468656/147‑1 (MQ=255)
                                                                             tCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaa                                                                       >  1:829496/1‑150 (MQ=255)
                                                                               attattGTTGCGCAAACCCATGCCCCAAGGCTGAAAAAGGGAAAATGGCGGATCACCCTTTTCATCCCCCCCCACGGCGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaaa                                                                     <  2:57764/150‑1 (MQ=255)
                                                                                ttattGTGGCGAAAACCCAGCCGCCAAGCGTAAAAAAAGGCAAAGGGCGGAGCACACTTTGCATCCCCCCCGGCGGCGTGGAACCAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGAc                                                                             <  1:704320/141‑1 (MQ=255)
                                                                                   ttgttgCGCCAACCCATGCCCCAGGTCTGAAAAAAGGAAAATTGCGGATCACACTGTTCATTCCCCCCCCCGGCGTGGAAACGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAtt                                                                  <  1:649558/149‑1 (MQ=255)
                                                                                   ttgttgCGCCAAACCATGCGCCAAGGCTGTAATAAGGAAAATTGCGGATCACACTGTTCATCCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAtt                                                                  <  2:501770/149‑1 (MQ=255)
                                                                                   ttgttgCGCAAAACCATGCCCCAAGCTTGTAATAAGCCAAATTGCGGATCACCTTTTTCATCCCCCCCCCCGGGGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGc                                                                          <  2:804868/141‑1 (MQ=255)
                                                                                   ttgtggCGCAAACCCATGCGCAAGGGCTGTAAAAAGGCAAATTGCGGATCACACTGTTCATCCCCCCCCACGGCGGGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAttt                                                                 <  1:392328/150‑1 (MQ=255)
                                                                                        gcgcAAAACCTTGCGCCAAGGCTTTAATAAGGAAAATTGGGGATCACTCTTTTCATTCCCCACCACGGGGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATa                                                           <  2:40238/151‑1 (MQ=255)
                                                                                         cgcCAACCCATGCGCCAAGGCTGAAAAAAGGCAAATTGCGGATCACCCTTTTCATCCCCCCCCACGGCGGGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATa                                                           <  2:1192204/150‑1 (MQ=255)
                                                                                         cgcCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAc                                                          >  1:1186443/1‑151 (MQ=255)
                                                                                          gcCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAcc                                                         >  1:570830/1‑151 (MQ=255)
                                                                                           ccAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAcc                                                         >  2:1197015/1‑150 (MQ=255)
                                                                                             aaaCCAGGCGCCAAGTCTGAAAAAAGGAAAATTGCGGACCACTCTGTTCATCCCCCCCCACGGCGGGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACg                                                      <  1:204930/151‑1 (MQ=255)
                                                                                              aaCCATGCCCCAAGGCTGTAAAAAGGCAAATGGCGGATCACACTGTTCATCCCCCCCCGCGGCGTGGGAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCATGa                                                     <  1:199531/151‑1 (MQ=255)
                                                                                                caaagcgcAAGGGCTGAAAGAAGGAAATGGGCGGTTCCCACTGTTCTTTCCCCACCACGGGGTGGAAAAGAAATTAGCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGaaa                                                   <  2:97993/147‑1 (MQ=255)
                                                                                                 cATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGCAAt                                                  >  2:373259/1‑151 (MQ=255)
                                                                                                  aTGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAAt                                                  >  2:980449/1‑150 (MQ=255)
                                                                                                  aTGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAAt                                                  >  2:1218531/1‑150 (MQ=255)
                                                                                                       caaaGGCTGTAAAAAGGAAAATTGCGGATCCCATTTTTCATCCCCCCCCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGa                                               <  2:208999/146‑1 (MQ=255)
                                                                                                        cAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAAt                                           >  2:16919/1‑151 (MQ=255)
                                                                                                                 aataaGGAAAATTGCGGATCACCCTGTTCATCCCCCCCCACGGGGTGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATa                                  <  1:464676/151‑1 (MQ=255)
                                                                                                                 aaaaaGGAAAATTGGGGATCCCCCTTTTCATCCCCCCCCCCGGGGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCAt                                   <  1:959135/147‑1 (MQ=255)
                                                                                                                 aaaaaGCCAATTTGCGGACCACATTGTTCATCCCCCCCCCGGGGGTGGGAAAAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACg                                      <  1:1110929/144‑1 (MQ=255)
                                                                                                                    aaGGAAAATTGCGGATCACACTTTTCATCCCCCCCCCGGGGGGGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACa                                <  1:1241584/150‑1 (MQ=255)
                                                                                                                         aatttGGGGATCACACTGTTCATCCCCCCCCACGGGGTGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATa                            <  2:511562/146‑1 (MQ=255)
                                                                                                                         aaaTTGGGGACCACGCTTTTCTTCCCCCCCCACGGGGTGGGAAGAAATTAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGa                              <  1:203028/147‑1 (MQ=255)
                                                                                                                         aaaTTGCGGACCCCCCTTTTCATCCCCCCCCCCGGGGTGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGAt                             <  1:1260043/148‑1 (MQ=255)
                                                                                                                         aaaTTGCGGACCACACTGTTCATCCCCCCCCACGGCGGGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAcc                          <  1:121694/151‑1 (MQ=255)
                                                                                                                         aaaTTGCGGAACACACTGTTCATCCCCCCCCACGGCGGGGGAAAAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAc                           <  2:829148/150‑1 (MQ=255)
                                                                                                                         aaaTGGCGGATCCCTCTTTTCTTCCCCCCCCTCGGCGTGGAACCGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAcc                          <  2:308653/151‑1 (MQ=255)
                                                                                                                            ttGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACcac                        >  1:634909/1‑150 (MQ=255)
                                                                                                                             tGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACcaca                       >  1:531602/1‑150 (MQ=255)
                                                                                                                               tggggATCACTTTTTTCTTCCCCCCCCCCGGCGGGGGAGAAAATTTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACca                         <  1:1051421/145‑1 (MQ=255)
                                                                                                                               gCGGATCACCCTTTTCATCCCCCCCCACGGGGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACaa                      <  2:144402/150‑1 (MQ=255)
                                                                                                                                gcgcatccctttGTTCATCCCCCCCCCCGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACaa                      <  2:1249921/139‑1 (MQ=255)
                                                                                                                                   gATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTg                  >  2:478596/1‑151 (MQ=255)
                                                                                                                                     tCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGt                >  1:1272592/1‑151 (MQ=255)
                                                                                                                                        cacTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCt             >  1:790599/1‑151 (MQ=255)
                                                                                                                                            gTTCTTCCCCCCCCACGGCGTGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTgg           <  1:556210/149‑1 (MQ=255)
                                                                                                                                             ttCATCCCCCCCCACGGGGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCg         <  1:457572/150‑1 (MQ=255)
                                                                                                                                                   ccccccccCGGGGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGatca  <  2:255663/151‑1 (MQ=255)
                                                                                                                                                   ccccccccCGGCGTGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGatc   <  2:340365/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
CCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCA  >  NC_000913/1174461‑1174756

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: