Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8800TC42.1% 92.4 / 14.3 57K271E (AAA→GAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (26/7);  new base C (0/24);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 3.60e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357731‑358026
                                                                                                                                                     |                                                                                                                                                       
cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc                                                                                                                                                        >  2:749319/1‑151 (MQ=255)
cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc                                                                                                                                                        >  2:423311/1‑151 (MQ=255)
  cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAGAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                                      >  1:737974/1‑151 (MQ=255)
  cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                                      <  1:424421/151‑1 (MQ=255)
  cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                                       <  2:631388/150‑1 (MQ=255)
      ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                                   <  2:233697/150‑1 (MQ=255)
      ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                                                                                                                  <  1:916575/151‑1 (MQ=255)
           gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                              <  2:354731/149‑1 (MQ=255)
              aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                         >  2:531989/1‑151 (MQ=255)
                gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                       >  1:103575/1‑151 (MQ=255)
                gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                       <  1:127086/151‑1 (MQ=255)
                               cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                         >  1:206278/1‑150 (MQ=255)
                                   tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                                    >  1:1030446/1‑151 (MQ=255)
                                    ccggCCGCAGACTTGTGTCAGCGCCTCCCGTGATTATTTCCCGCCAACACCCGCGCTTTTACCCACCGCCATGTGAAAAATGGCCTTCGCGGCGGCTGTCGCCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagta                                                                                                                        <  2:1135744/143‑1 (MQ=255)
                                    ccgccccaatccttgtgTCGGCGGTTACGGTGTTTCTTTCCCCCCAACACCCGCGCATTTATCCACCCCCTTGTGAAAAAAGGCCTTCGCGGCGGCTTTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                                    <  2:629497/137‑1 (MQ=255)
                                    ccAGCCGCAGCCGTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACACCCGCGCTTTTATCCACCGCCATGTGCAAAAAGGCCTGCGCGGCGGCTTTCTCCCAGTTTTTTTCCCGCCGCAACAAACCCGCCGTTCTCTCCagtagt                                                                                                                      <  2:601833/148‑1 (MQ=18)
                                    ccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                   >  1:352860/1‑151 (MQ=255)
                                     cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                 >  1:469389/1‑151 (MQ=255)
                                          gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggg                                                                                                            >  1:1053028/1‑149 (MQ=255)
                                          gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                               >  2:907298/1‑148 (MQ=255)
                                             gACCTGTGTCAGCGGCTCCCGTGATTCATTCCCGCCAACACCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTTTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                          <  2:259372/150‑1 (MQ=255)
                                                 tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGaa                                                                                                           >  2:1089006/1‑145 (MQ=255)
                                                 tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaa                                                                                                     >  2:274660/1‑150 (MQ=255)
                                                 tgtgTCAGCGGCTACCGTAATTCATTCCCCCCACCAACCGCGCATTTACCCAACCCTATGTAAAAAAATCCCTCCGCGGCGCCTGTCTGCCGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                       <  1:712839/148‑1 (MQ=255)
                                                   tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggcaagg                                                                                                        >  1:401349/1‑140 (MQ=255)
                                                     tCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggcggggaaaa                                                                                                   >  2:581736/1‑138 (MQ=255)
                                                      cAGCGGCTCCCGTGATTCTTTCCCCCCAACACCCGCGCTTTTCCCCACCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                                                                                                  <  1:893008/148‑1 (MQ=21)
                                                         cGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggaaaggaa                                                                                                     >  1:456186/1‑134 (MQ=255)
                                                          ggCTCCCGTGATTCTTTCCCCCCAACACCCGCGCATTTACCCACCGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                                                                                           <  2:257026/151‑1 (MQ=255)
                                                          ggCTCCCGTATTTCATTCCCCCCAACACCCGCGCATTTACCCACCGCCATGTGCAAAAATGCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                                                   <  2:352865/143‑1 (MQ=255)
                                                          ggCTACCGTAATCCATTCCCGCAAACACCCCCGCTTTTACCCACCGCCTTGTACAAAAAGGCCTTCGGGGCGGTTGCCTGCCAGTTATTTTCCCGCGACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTa                                                                                            <  1:560066/150‑1 (MQ=255)
                                                            ctccggtgtTTCTTTCCCCCAAACACCCGCGCTTTTCCCCACCGCCATGTCAAAAAATCCCTCCGCGGCGGCTTTCCGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                                                                                           <  1:650957/144‑1 (MQ=255)
                                                               cggtgtTTCTTTCCCGCAACCACCCGCGCATTTCCCCACCGCCATGTCCAAAAATCCCTCCGGGGCGGCTGCCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGccc                                                                                       <  2:504087/148‑1 (MQ=255)
                                                               ccgtttttCTTTCCCGCCAACACCCGCTTATTTATCCACCGCCTTGTGCAAAAATGCCTTCGCGGGGGCTGTCGCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg                                                                                          <  1:706254/141‑1 (MQ=255)
                                                               ccGTGTTTCATTCCCCCAAACACCCGCCAATTTACCCACCCCCATGTAAAAAAAGGCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCg                                                                                      <  2:648613/151‑1 (MQ=255)
                                                                cGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaaaaaattttaaccc                                                                                       >  2:617142/1‑136 (MQ=25)
                                                                          ttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGAAAAAAATTTTAACCCCTCAtgtt                                                                               >  2:96361/1‑147 (MQ=25)
                                                                                  acaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAGGAAAAACCTTTAGCCCGGCAtggttt                                                                             >  2:2500/1‑137 (MQ=25)
                                                                                    acccGCGCATTTACCCACCGCCATGTAAAAAAATGCTTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                                                   <  2:691192/147‑1 (MQ=255)
                                                                                           cATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGGAAAAACTTTAAACCCCCCCTGTTTTGTGGGGATCGCCGCTgg                                                           >  2:1106340/1‑149 (MQ=37)
                                                                                             tttccccACCGCCATGTAAAAAATGGCCTCCGCGGCGCCTGTCTGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                                <  1:42919/138‑1 (MQ=255)
                                                                                                    accccctgtgaaaaaaaTGCCTCCGCGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                                  <  1:821399/136‑1 (MQ=255)
                                                                                                        cagtgtcaaaaaaTGCCTCGGCGGCGGCTTCCTCCCGGTTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAg                                           <  1:252244/143‑1 (MQ=255)
                                                                                                            gtgAAAAAATGCTTCCGGGGCGGTTGTCTGCCACTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                         <  2:796812/149‑1 (MQ=255)
                                                                                                            gtgAAAAAATGCGTGCGCGGCGCTTGCCGCCGTGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                                 <  1:774316/141‑1 (MQ=255)
                                                                                                              gCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAGCTTTTAAACCCTCCTGTTTTGTGGGGAACGCCGCTTGGTAAAAAGTGGGAAggg                                         >  2:156126/1‑147 (MQ=255)
                                                                                                                aaaaaTGCCTGCGGGGGGGTTTTCTGCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                     <  1:738208/149‑1 (MQ=255)
                                                                                                                 aaaaTGCTTGCGCGGCGGCTGTCGGCAAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg                                  <  2:830004/151‑1 (MQ=255)
                                                                                                                 aaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGGAAAAAATTTTAAACCCCCCTGTTGTGTGGGGAACGCCGCTGGGGAAAAAATGGAAAGGGGAggggg                                  >  1:351484/1‑147 (MQ=255)
                                                                                                                   aaggCCTCCGGGGGGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggc                                <  1:815605/148‑1 (MQ=255)
                                                                                                                       ccTTTGCGCCGCCTTTCTGCCACCTTTTTTCCCGCAGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                            <  2:175171/151‑1 (MQ=255)
                                                                                                                           cgcgGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAAACTTTAAACCGTCCTGGTGTGTGGGTATAGCCGCTGGTGACAAATTTGGAAggggggggggggg                                >  1:912286/1‑135 (MQ=255)
                                                                                                                                    tGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAGGAGAAATAAATTTTAACCCCTTACGTTTTGGGGGAAACGCTGCCGGTGAGAACGGGGGAAGGGAGGAGGGGGGGGTCAGCTctccaa                 >  2:837820/1‑144 (MQ=255)
                                                                                                                                      tCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc               <  2:382457/149‑1 (MQ=255)
                                                                                                                                        tGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGCCCGGGGAAAAAATTTAAATCCCCCCTGTTTTGTGGGGAACTCTGGTGGGGACAAAATGGGAAGGGGAGGGGGGGGGGACAGCTCCAGAAAcgc             >  1:497982/1‑149 (MQ=255)
                                                                                                                                         gCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAGAAAAAGATTTTAACCCCCCCTGTTTTGTGGGAAACGCTGCTGGTGAAAAAATGGGAAGGGGAGGGGGGCGGGTCAACTCCCCAAcccc             >  2:581556/1‑146 (MQ=255)
                                                                                                                                             gttATTTTTCGCCCCAAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc             <  2:480100/142‑1 (MQ=255)
                                                                                                                                                  ttttCCCCCCCAAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa    <  1:617129/149‑1 (MQ=255)
                                                                                                                                                  tttcCCGCCGCAACAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt  <  2:832696/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357731‑358026

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