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

CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  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)
           gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGc‑aac                                                                                                                                              <  2:354731/149‑1 (MQ=255)
              aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAAc                                                                                                                                         >  2:531989/1‑151 (MQ=255)
                gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAAccg                                                                                                                                       >  1:103575/1‑151 (MQ=255)
                gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAAccg                                                                                                                                       <  1:127086/151‑1 (MQ=255)
                               cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCagt                                                                                                                         >  1:206278/1‑150 (MQ=255)
                                   tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                                    >  1:1030446/1‑151 (MQ=255)
                                    ccggCCGCAGACTTGTGTCAGCGCCTCCCGTGATTATTTCCCGCCAACACCCGCGCTTTTACCCACCGCCATGTGAAAAATGGCCTTCGCGGCGGCTGTCGCCCAGCTATTTTCCCGCCGC‑AAAAAAACCGCCGTTCTCTCCagta                                                                                                                        <  2:1135744/143‑1 (MQ=255)
                                    ccgccccaatccttgtgTCGGCGGTTACGGTGTTTCTTTCCCCCCAACACCCGCGCATTTATCCACCCCCTTGTGAAAAAAGGCCTTCGCGGCGGCTTTCTGCCAGCTTTTTTCCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                                    <  2:629497/137‑1 (MQ=255)
                                    ccAGCCGCAGCCGTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACACCCGCGCTTTTATCCACCGCCATGTGCAAAAAGGCCTGCGCGGCGGCTTTCTCCCAGTTTTTTTCCCGCCGC‑AACAAACCCGCCGTTCTCTCCagtagt                                                                                                                      <  2:601833/148‑1 (MQ=18)
                                    ccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                   >  1:352860/1‑151 (MQ=255)
                                     cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑g                                                                                                                 >  1:469389/1‑151 (MQ=255)
                                          gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggggg                                                                                                            >  1:1053028/1‑149 (MQ=255)
                                          gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggg                                                                                                               >  2:907298/1‑148 (MQ=255)
                                             gACCTGTGTCAGCGGCTCCCGTGATTCATTCCCGCCAACACCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTTTCTGCCAGCTATTTTCCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAag                                                                                                          <  2:259372/150‑1 (MQ=255)
                                                 tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGaa                                                                                                           >  2:1089006/1‑145 (MQ=255)
                                                 tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGA‑Gaa                                                                                                     >  2:274660/1‑150 (MQ=255)
                                                 tgtgTCAGCGGCTACCGTAATTCATTCCCCCCACCAACCGCGCATTTACCCAACCCTATGTAAAAAAATCCCTCCGCGGCGCCTGTCTGCCGGTTTTTTTCCCGCCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAaga‑g                                                                                                       <  1:712839/148‑1 (MQ=255)
                                                   tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggggcaagg                                                                                                        >  1:401349/1‑140 (MQ=255)
                                                     tCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggggggcggggaaaa                                                                                                   >  2:581736/1‑138 (MQ=255)
                                                      cAGCGGCTCCCGTGATTCTTTCCCCCCAACACCCGCGCTTTTCCCCACCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATa                                                                                                  <  1:893008/148‑1 (MQ=21)
                                                         cGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggggggaaaggaa                                                                                                     >  1:456186/1‑134 (MQ=255)
                                                          ggCTCCCGTGATTCTTTCCCCCCAACACCCGCGCATTTACCCACCGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTaa                                                                                           <  2:257026/151‑1 (MQ=255)
                                                          ggCTCCCGTATTTCATTCCCCCCAACACCCGCGCATTTACCCACCGCCATGTGCAAAAATGCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAAt                                                                                                   <  2:352865/143‑1 (MQ=255)
                                                          ggCTACCGTAATCCATTCCCGCAAACACCCCCGCTTTTACCCACCGCCTTGTACAAAAAGGCCTTCGGGGCGGTTGCCTGCCAGTTATTTTCCCGCGAC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTa                                                                                            <  1:560066/150‑1 (MQ=255)
                                                            ctccggtgtTTCTTTCCCCCAAACACCCGCGCTTTTCCCCACCGCCATGTCAAAAAATCCCTCCGCGGCGGCTTTCCGCCAGTTTTTTTCCCGCCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTaa                                                                                           <  1:650957/144‑1 (MQ=255)
                                                               cggtgtTTCTTTCCCGCAACCACCCGCGCATTTCCCCACCGCCATGTCCAAAAATCCCTCCGGGGCGGCTGCCTGCCAGTTTTTTTCCCGCCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGccc                                                                                       <  2:504087/148‑1 (MQ=255)
                                                               ccgtttttCTTTCCCGCCAACACCCGCTTATTTATCCACCGCCTTGTGCAAAAATGCCTTCGCGGGGGCTGTCGCCCAGTTTTTTTCCCGCCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAg                                                                                          <  1:706254/141‑1 (MQ=255)
                                                               ccGTGTTTCATTCCCCCAAACACCCGCCAATTTACCCACCCCCATGTAAAAAAAGGCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCg                                                                                      <  2:648613/151‑1 (MQ=255)
                                                                cGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGA‑Gaaaaaaattttaaccc                                                                                       >  2:617142/1‑136 (MQ=25)
                                                                          ttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCGAGA‑GAAAAAAATTTTAACCCCTCAtgtt                                                                               >  2:96361/1‑147 (MQ=25)
                                                                                  acaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAG‑GAAAAACCTTT‑AGCCCGGCAtggttt                                                                             >  2:2500/1‑137 (MQ=25)
                                                                                    acccGCGCATTTACCCACCGCCATGTAAAAAAATGCTTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                                                   <  2:691192/147‑1 (MQ=255)
                                                                                           cATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGA‑GGAAAAACTTTAAACCCCCCCTGTTTTGTGGGGATCGCCGC‑Tgg                                                           >  2:1106340/1‑149 (MQ=37)
                                                                                             tttccccACCGCCATGTAAAAAATGGCCTCCGCGGCGCCTGTCTGCCAGTTTTTTTCCCGCCAA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                                <  1:42919/138‑1 (MQ=255)
                                                                                                    accccctgtgaaaaaaaTGCCTCCGCGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAAC‑AAtg                                                  <  1:821399/136‑1 (MQ=255)
                                                                                                        cagtgtcaaaaaaTGCCTCGGCGGCGGCTTCCTCCCGGTTTTTTTCCCGCCGA‑AAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAAC‑AATGTGGAAAg                                           <  1:252244/143‑1 (MQ=255)
                                                                                                            gtgAAAAAATGCTTCCGGGGCGGTTGTCTGCCACTTTTTTTCCCGCCGA‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAAC‑AATGTGGAAAggg                                         <  2:796812/149‑1 (MQ=255)
                                                                                                            gtgAAAAAATGCGTGCGCGGCGCTTGCCGCCGTGCTTTTTTCCCGCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAAC‑AAtgt                                                 <  1:774316/141‑1 (MQ=255)
                                                                                                              gCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAGA‑GAAAAAGCTTTTAAACCCTCCTGTTTTGTGGGGAACGCCGCTTGGTAA‑‑AAAGTGGGAAggg                                         >  2:156126/1‑147 (MQ=255)
                                                                                                                aaaaaTGCCTGCGGGGGGGTTTTCTGCCATTTTTTTTCCCGCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAAC‑AATGTGGAAAGGGAAGt                                     <  1:738208/149‑1 (MQ=255)
                                                                                                                 aaaaTGCTTGCGCGGCGGCTGTCGGCAAGTTTTTTTCCCCCCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAAC‑AATGTGGAAAGGGAAGTgcg                                  <  2:830004/151‑1 (MQ=255)
                                                                                                                 aaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGG‑GAAAAAATTTTAAACCC‑CCCTGTTGTGTGGGGAACGCCGC‑TGGGGAA‑AAAATGGAAAGGGGAggggg                                  >  1:351484/1‑147 (MQ=255)
                                                                                                                   aaggCCTCCGGGGGGGCTGTCTGCCAGCTTTTTTCCCGCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAAC‑AATGTGGAAAGGGAAGTgcggc                                <  1:815605/148‑1 (MQ=255)
                                                                                                                       ccTTTGCGCCGCCTTTCTGCCACCTTTTTTCCCGCAGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAAt                            <  2:175171/151‑1 (MQ=255)
                                                                                                                           cgcgGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGGAAAAAACTTTAA‑ACCGTCCTGGTGTGTGGGTATAGCCGC‑TGGTGACAAATTTGGAAggggggggggggg                                >  1:912286/1‑135 (MQ=255)
                                                                                                                                    tGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAGGA‑GAAATAAATTTTAACCCCTTACGTTTTGGGGGAAACGCTGC‑CGGTGAG‑AACGGGGGAAGGGAGGAGGGGGGGGTCAGCTctccaa                 >  2:837820/1‑144 (MQ=255)
                                                                                                                                      tCTGCCAGTTTTTTTCCCGCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc               <  2:382457/149‑1 (MQ=255)
                                                                                                                                        tGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGTGCCCGG‑GGAAAAAATTTAAATCCCCCCTGTTTTGTGGGGAACTCTGG‑TGGGGAC‑AAAATGGGAAGGGGAGGGGGGGGGGACAGCTCCAGAAAcgc             >  1:497982/1‑149 (MQ=255)
                                                                                                                                         gCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAAA‑GAAAAAGATTTTAACCCCCCCTGTTTTGTGGGAAACGCTGC‑TGGTGAA‑AAAATGGGAAGGGGAGGGGGGCGGGTCAACTCCCCAAcccc             >  2:581556/1‑146 (MQ=255)
                                                                                                                                             gttATTTTTCGCCCCA‑AACAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc             <  2:480100/142‑1 (MQ=255)
                                                                                                                                                  ttttCCCCCCCAAACAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa    <  1:617129/149‑1 (MQ=255)
                                                                                                                                                  tttcCCGCCGC‑AACAAATCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt  <  2:832696/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  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: ‑