Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,556,2000GA45.6% 95.1 / 52.4 93S217F (TCC→TTC) rtcARNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base G (40/10);  new base A (0/42);  total (40/53)
Fisher's exact test for biased strand distribution p-value = 5.51e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCC  >  NC_000913/3556051‑3556343
                                                                                                                                                            |                                                                                                                                               
ccAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGa                                                                                                                                                <  1:521908/151‑1 (MQ=255)
     tgtgCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAAc                                                                                                                                           >  1:316055/1‑151 (MQ=255)
      gtgCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACt                                                                                                                                          <  2:110324/151‑1 (MQ=255)
      gtgCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACt                                                                                                                                          <  1:217201/151‑1 (MQ=255)
      gtgCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACt                                                                                                                                          >  1:1312740/1‑151 (MQ=255)
        gCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTcc                                                                                                                                        >  2:1262653/1‑151 (MQ=255)
        gCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTc                                                                                                                                         <  1:711049/150‑1 (MQ=255)
         ccGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTccc                                                                                                                                       <  1:220912/151‑1 (MQ=255)
            cAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGc                                                                                                                                    >  1:1082582/1‑151 (MQ=255)
              acccccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGcc                                                                                                                                   <  2:314239/150‑1 (MQ=255)
                 accTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGCGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAg                                                                                                                                 >  2:1226314/1‑149 (MQ=255)
                    tCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAgtgt                                                                                                                              >  2:246688/1‑149 (MQ=255)
                                cttttttCCCCGCCAAAAAAAAAGCGTTG‑GGTGATATTTTCTCTTTCAACTCAAAGCGAGGCGGTATCCCCC‑GGCCC‑CTGG‑CGGCGCGGAGGG‑TTAAGAATTTT‑TTGTT‑CATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACg                                                                                                                  <  1:1236951/146‑1 (MQ=255)
                                cgcTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCATGc                                                                                                                 >  1:1300234/1‑150 (MQ=255)
                                 gcTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                               >  2:194475/1‑151 (MQ=255)
                                   tttcccccGACAAAAAAAAGGGTTTC‑GTTAATTTTTTCCTTTTCGACTTAAGGCGAGACGTTATTCCCC‑GGCCC‑CTGT‑TCGCGCGGCAGT‑TTATGAATTTT‑CTGTT‑CATGCGGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                               <  1:736380/143‑1 (MQ=255)
                                    ttttccccaagaaaaaaaaGTGTTCGGGTGATATTTCCCCTTTCGCCTTCAACCGAGCCGTTATTCCCG‑GCCCC‑CTGG‑CCGCGCGGCAGG‑TTAAAAAAATT‑TTTTT‑CAGGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCa                                                                                                              <  1:146980/137‑1 (MQ=255)
                                    tttCCCCGACGAAAAAAAAGCGTTG‑GGTAATTTTTTCCTTTTCTACTCAAAGCGAGAGGGTTTTCCCC‑GGCCC‑CTGG‑TCGCGGGGCAGT‑TTATGAATATT‑TTTTT‑CATGAAGGGAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAg                                                                                                             <  1:59341/150‑1 (MQ=255)
                                    tttCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGc                                                                                                            >  2:868161/1‑151 (MQ=255)
                                       cACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAAc                                                                                                         >  2:854480/1‑151 (MQ=255)
                                         ccGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAt                                                                                                       >  1:191761/1‑151 (MQ=255)
                                         ccGACAAAAAAAAGGCGTGC‑GGTGATTTTTCCACTTCCGACTTCAAGCGAGCCGTTATTCCCC‑GCCCC‑TTGG‑TCGCGGGGCAGT‑TTATAAATATT‑CTGTT‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACa                                                                                                        <  1:869360/150‑1 (MQ=255)
                                          cgacaaAAAAAAAGGGTTG‑GTTGATTTTTTCGCTTCCGCCTAAAACAGAGAGGTTATTCCCC‑GCCCG‑CGGG‑GCGCGCGGAAGT‑TAATAAATATT‑TTGTT‑CAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg                                                                                                      <  2:434108/151‑1 (MQ=255)
                                              gacAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                                  >  1:648531/1‑151 (MQ=255)
                                               acAAAAAGGCGTTG‑GGTAATTTTTCCCCTTCCGCCTTCAAGCGAGCCGGTATCCCCC‑GCCCC‑CGGG‑TCGCGGGGCAGT‑TTATGAATTTT‑CTGTT‑CAGGAAGAGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgg                                                                                                     <  2:151643/147‑1 (MQ=255)
                                                cAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCgg                                                                                                >  1:17437/1‑151 (MQ=255)
                                                 aaaaaaGTGTCG‑GGTAATATTTTCCCTTCCGACTCAAAGGGAGCGGGTATTCCCC‑GGCCC‑CGGG‑TCGCCCGGCAGG‑TTATGAATTTT‑CTGTT‑CAGGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                               <  1:255222/151‑1 (MQ=255)
                                                 aaaaaaGCTTTC‑GGTATTTTTTTCACTTGCGGCTGAAGAGGGGACGGTATTCCCC‑GCCCC‑CGGG‑CCGCG‑GGCGGGTTTTTGATTTTT‑CTGTT‑CACGGGGAAAAAAACACCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                 <  1:539026/149‑1 (MQ=255)
                                                 aaaaaaGCGTTC‑GTTGTTTTTTCCCCTTTCGACTAAAAGGGAGACGGTATT‑CCC‑CGCCC‑CTGGTTCCCGCGGCAGG‑TTAGAAATATT‑CGGTT‑CAGGCAGGAAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                               <  2:1147321/151‑1 (MQ=255)
                                                 aaaaaaGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                               >  2:140756/1‑151 (MQ=255)
                                                 aaaaaaGCGTTC‑GGGAATATTTTCACTTTCGCCTCCAAGCGAGACGGTATTCCCC‑GCCCC‑CGGG‑TCGCGCGGCAGT‑TTATAAATATT‑CTGTT‑CAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                               <  1:1064369/151‑1 (MQ=255)
                                                 aaaaaaGCGTCG‑GGTGATTTTTTCCCTTCCGTCTACAAGGGGGACGGTATTCCCC‑GCCCC‑CTGC‑CCGCGGGGCAGT‑TTATGAATTTT‑CTGTC‑CAGGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                               <  1:426650/151‑1 (MQ=255)
                                                  aaaaaGCTTTC‑GGTGATATTTTCTTTTCCGACTTCAAGCGAGCCGGTATCCCCC‑GCCCC‑CGGG‑TCGCGCGGCAGT‑TTATAAATTTT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                              <  1:271383/151‑1 (MQ=255)
                                                   aaaaGCGTTG‑GGTGATTTTTCCCCTTTCGCCTTAAAGCGAGACGGTATTCCCC‑GCCCC‑GTGG‑TCGCGCGGAAGT‑TTATAAATTTT‑CTGTT‑AATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcac                                                                                             <  1:1494724/151‑1 (MQ=255)
                                                    aaaGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcac                                                                                             >  1:684368/1‑150 (MQ=255)
                                                     aaGCGTTC‑GGTGATTTTTTCCCTTTCGACTTAAAGCGAGACGGTTTTCCCC‑GGCCC‑CTGG‑TCGCGCGGCAGT‑TTATAAATATT‑CTGTT‑AATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCacac                                                                                           <  1:480563/151‑1 (MQ=255)
                                                          ttC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCAAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAg                                                                                        >  2:84309/1‑149 (MQ=255)
                                                              ggTGATTTTTTCATTTTCGACTTCAAGCGAGACGTTATTCCCC‑GGCCC‑CGGG‑TCGCGCGGCAGT‑TTATAAATTTT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                      <  1:1501601/148‑1 (MQ=255)
                                                               gTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAAt                                                                                   >  1:1278163/1‑150 (MQ=255)
                                                                 gATTTTTTCCCTTTCGCCTTCAAGCGAGACGGTATTCCCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATa                                                                                  <  1:1159522/149‑1 (MQ=255)
                                                                  atatTTTCTCTTTCGACTTCAAGCGAGACGGTATTCCCC‑GGCCC‑CTGG‑TCGCGCGGCGGT‑TTATGAATTTT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa                                                                                <  1:1183612/150‑1 (MQ=255)
                                                                   ttttttCTCTTCCGCCTACAAGCGAGACGTTATTCCCC‑GCCCC‑CTGG‑TCGCGCGGCAGT‑TTATGAATTTT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa                                                                                <  1:70533/147‑1 (MQ=255)
                                                                   ttttttCCCTTTCGACTTAAAGGGAGACGTTATTCCCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTTTT‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg                                                                                 <  1:889041/146‑1 (MQ=255)
                                                                   ttttttCACTTCCGACTTCAAGCGAGACGGTTTTCCCC‑GGCCC‑CTGG‑CCGCGGGGCAGT‑TTAGGAATTTT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTa                                                                                     <  2:698479/142‑1 (MQ=255)
                                                                   tatTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                              >  2:1380111/1‑151 (MQ=255)
                                                                     ttttCCTTTTCGCTTTCAAGCGGGACGTTTTCCCCC‑GGCCC‑CTGG‑CCGCGCGGCGGT‑TTATAAATTTT‑CTGTT‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg                                                                                 <  2:1005951/146‑1 (MQ=255)
                                                                     ttttCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACtt                                                                            >  1:214451/1‑151 (MQ=255)
                                                                      tttCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtc                                                                           >  2:964792/1‑151 (MQ=255)
                                                                         cACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTcc                                                                        >  1:1157525/1‑151 (MQ=255)
                                                                          aCTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTcc                                                                        >  1:463437/1‑150 (MQ=255)
                                                                           cTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCAc                                                                      >  1:867618/1‑151 (MQ=255)
                                                                            tttCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                     >  2:558435/1‑151 (MQ=255)
                                                                             ttCGCCTTCAAGCGAGGCGGTATTCCCC‑GCCCC‑CGGG‑TCGCGGGGCAGG‑TTATGAATTTT‑CTTTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                     <  1:457597/150‑1 (MQ=255)
                                                                             ttCGACTTCAAGCGAGACGGTATTCCCC‑GGCCC‑CTGG‑TCGCGCGGCAGT‑TTATGAATTTT‑CTGTT‑AATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                     <  1:1471712/150‑1 (MQ=255)
                                                                               cGACTCCAAGCGAGACGGTTTTCCCC‑GCCCC‑CTGG‑CCGCGCGGCAGT‑TTATAAATTTT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa                                                                   <  1:295877/150‑1 (MQ=255)
                                                                                 aCTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                >  2:844646/1‑151 (MQ=255)
                                                                                 aCTTAAAGCGAGGCGGTATTCCCC‑GGCCC‑CTGG‑CCGCGCGGCGGT‑TTATGAATTTT‑CTGTT‑CAGGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCAc                                                                      <  1:374326/145‑1 (MQ=255)
                                                                                  cTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                               >  1:598991/1‑151 (MQ=255)
                                                                                   ttaaaGCGAGCCGGTATTCCCC‑GGCCC‑CTGG‑TCGCGCGGCAGT‑TTATAAATTTT‑CTGTT‑CATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                <  1:194664/146‑1 (MQ=255)
                                                                                   tcaaaGCGAG‑CGGTTTTTCCCGGGCCC‑CTGG‑TCGCGCGGCGGT‑TTATAAATTTT‑CTTTC‑CAGGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa                                                                   <  1:798265/143‑1 (MQ=255)
                                                                                   taaaggCGAGACGGTATTCCCG‑GGCCC‑CTGG‑CCGCGGGCAAGT‑TTATAAATTTT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                <  2:1356657/144‑1 (MQ=255)
                                                                                     cAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                            >  2:1424216/1‑151 (MQ=255)
                                                                                      aggCGAGACTGTATTCCCC‑GCCCC‑CTGG‑TCGCGGGGCAGT‑TTATAAATTTT‑CTGTT‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                            <  1:618276/148‑1 (MQ=255)
                                                                                      aaGCGAGACGTTATTCCCC‑GGCCC‑CTGG‑TCGCGCGGCAGT‑TTATGAATATT‑CTTTT‑CATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcaca                                                           <  1:1000512/151‑1 (MQ=255)
                                                                                      aaGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                            >  2:1258614/1‑150 (MQ=255)
                                                                                          gagaCGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                       >  1:225779/1‑151 (MQ=255)
                                                                                           agaCGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt                                                      >  1:943573/1‑151 (MQ=255)
                                                                                             gccgatttcccc‑GGCCCTCTGC‑GCGCGCGGCGGT‑TGTTAAATTTT‑GTGTT‑CATGCGGCAAAACCCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATg                                                        <  2:531806/139‑1 (MQ=255)
                                                                                               ggTATCCCCC‑GGCCC‑CGGG‑TCGCGCGGCAGT‑TTATAAATATT‑CTGTT‑CAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                  <  2:744205/151‑1 (MQ=255)
                                                                                                  tttgcccgggCCC‑TTGG‑CCGCGCGTCGTG‑TTA‑GAATTTTCCTGTT‑CATGCAGGCAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccc                                                   <  2:48184/140‑1 (MQ=255)
                                                                                                   tccccc‑GGCCC‑CGGG‑TCGCGCGGGAGT‑TTATAAATTTT‑CTTTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                  <  2:1494489/144‑1 (MQ=255)
                                                                                                     cccc‑GCCCC‑CGGG‑CCGCGGGGAGGT‑TTATAAATTTT‑CTTTT‑CAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                  <  2:1089043/144‑1 (MQ=255)
                                                                                                      ccc‑GGCCC‑CTGG‑CCGCGCGGCAGG‑TTATAAATATT‑‑TGTTCCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                             <  2:889506/149‑1 (MQ=255)
                                                                                                      ccc‑GGCCC‑CGGG‑TCGCGCGGAAGT‑TTATAAATATT‑TTGTT‑CAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                             <  2:273097/149‑1 (MQ=255)
                                                                                                      ccc‑GCCCC‑CTGG‑TCGCGCGGCAGG‑TTATAAATATT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                             <  1:1415862/149‑1 (MQ=255)
                                                                                                             cc‑cTGG‑CCGCGCGGCAGT‑TTATAAATATT‑CTGTT‑CATGCAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGtt                                     <  1:1262671/151‑1 (MQ=255)
                                                                                                                       gcgcgGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgtg                             >  2:1224937/1‑151 (MQ=255)
                                                                                                                            ggcggtttATAAATATT‑CTGTT‑CATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaa                         <  2:1503162/145‑1 (MQ=255)
                                                                                                                            gCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaa                         >  1:311862/1‑150 (MQ=255)
                                                                                                                            gCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaa                         >  2:1482981/1‑150 (MQ=255)
                                                                                                                               gg‑TTATGAATATT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                       <  1:868177/149‑1 (MQ=255)
                                                                                                                                  ttATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATg                   >  2:340545/1‑151 (MQ=255)
                                                                                                                                  ttATGAATATT‑CTGTT‑CATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                     <  1:209819/149‑1 (MQ=255)
                                                                                                                                   tATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                  >  2:574375/1‑151 (MQ=255)
                                                                                                                                   tATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                  >  2:595155/1‑151 (MQ=255)
                                                                                                                                   tATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                  >  1:515550/1‑151 (MQ=255)
                                                                                                                                    aTGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc                 >  1:609764/1‑151 (MQ=255)
                                                                                                                                       aaTTTT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                <  2:463431/149‑1 (MQ=255)
                                                                                                                                       aaTATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg              >  1:647125/1‑151 (MQ=255)
                                                                                                                                          atT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc           >  2:1061257/1‑151 (MQ=255)
                                                                                                                                            t‑CTGTT‑CATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCg          <  1:1136124/150‑1 (MQ=255)
                                                                                                                                               tGTT‑CAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgag        <  2:1275822/150‑1 (MQ=255)
                                                                                                                                                       gCAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc  <  1:194084/149‑1 (MQ=255)
                                                                                                                                                            |                                                                                                                                               
CCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCC  >  NC_000913/3556051‑3556343

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: ‑