Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8800TC45.3% 74.2 / 22.1 64K271E (AAA→GAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (27/8);  new base C (0/29);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 4.46e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357731‑358026
                                                                                                                                                          |                                                                                                                                                      
cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTAttttt                                                                                                                                                        >  1:951325/1‑150 (MQ=255)
     aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑ccgccg                                                                                                                                                 <  2:228031/151‑1 (MQ=255)
      tgCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑Ccgccgc                                                                                                                                                >  3:229641/2‑151 (MQ=255)
       gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑Ccgccgc                                                                                                                                                >  1:514694/1‑150 (MQ=255)
          tgctAAGCGTAGCGCATCCGGCAATTCCAGACGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaa                                                                                                                                              >  1:89881/4‑149 (MQ=255)
          tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaac                                                                                                                                             >  4:296096/1‑150 (MQ=255)
            aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaacaa                                                                                                                                           >  1:782945/1‑150 (MQ=255)
            aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACaaa                                                                                                                                          <  3:358069/151‑1 (MQ=255)
              aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACaaaa                                                                                                                                         <  1:1443762/150‑1 (MQ=255)
              aaGCGTAGCGCATAAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACCACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGC‑CGGCAGTCCTCCAGCTATTTTT‑CCGCCGCAACaaaa                                                                                                                                         >  1:431767/1‑150 (MQ=18)
                gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAcc                                                                                                                                       >  1:287807/1‑150 (MQ=255)
                 cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTT‑ATCCACCGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgc                                                                                                                                     <  1:622152/151‑1 (MQ=255)
                    aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgccg                                                                                                                                   >  2:1414879/1‑150 (MQ=255)
                       gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTtct                                                                                                                               >  2:1487317/1‑151 (MQ=255)
                       gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTtc                                                                                                                                >  1:764947/1‑150 (MQ=255)
                           cagccattcccAGCCGGTGACTTGTGCCAGCGGCTACCGTGATTCATTCCCGCCACCACCGCCGCATTT‑ACCCACCCCCATGGG‑CAAAAAGG‑CCTT‑CGCGG‑CGGTTTTCTGCCAGCTTTTTTT‑CCGCCGCAACAAAACCGCCGTtc                                                                                                                                <  1:1625269/137‑1 (MQ=12)
                               cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                        >  2:51741/1‑150 (MQ=255)
                               acaTTCCACCCGCACACCTGTTCCAGTGTCTCCAGTTTTTAATCTCCGCAAACACCGGCGCATTT‑ATCC‑ACCCCATGTGCCAAAAGTG‑CCTT‑GCGGG‑CGCCTCTCTGCCAATATTTTTT‑CCCCCAAAAGAAAACCGCCGTTCTCTCCagt                                                                                                                        <  2:1740576/148‑1 (MQ=255)
                                attcccaccggCAGTCGTGTGTGGGCGGCCCCCTTTATTCTTTCCCCCAAACACCCGC‑CATTTCACCCACCGCAATGTG‑‑AAAAATGCCCTC‑CGGGG‑GGGTTGTCGCCCAGTTTTTTTC‑CCGCCCCAACAAAACCGCCGTTCTCTCCag                                                                                                                         <  1:1591554/138‑1 (MQ=255)
                                aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                        >  1:1354919/1‑149 (MQ=255)
                                    ccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                  >  4:391851/1‑151 (MQ=255)
                                    ccACCCGCAGCCTTGTGTCGGGGGCTCCGGTGTTTCTTTCCCCCCAACACCCGCGATTTT‑ACCCACCCCCTTGTG‑AAAAAATG‑CCTC‑CGGGG‑CGGCTTTCTCCCAGTTATTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                  <  2:964332/151‑1 (MQ=255)
                                      aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                >  2:1468973/1‑151 (MQ=255)
                                           cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggcg                                                                                                              >  2:1322673/1‑146 (MQ=255)
                                            agagttgggtCAGGGCCTCCCGTGATTTATCTCCGCAAACACCCGGCCATTT‑ACCCACCCCTATGTG‑CAAAAATG‑CTTT‑CGGGG‑GGGTTGTCCGCCAGCTTTTTTC‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                              <  4:76039/139‑1 (MQ=255)
                                            accccTGTGTCAGCGGCTCCCGTGTTTCTTCCCCGCAAACCCCCGCGCATTT‑ACCCACCGCCATGTA‑AAAAAACG‑CCTC‑CGCGG‑CGCCTGTCGCCCGGTTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                            <  1:1313644/146‑1 (MQ=255)
                                            aGACCTGTGCCGGCGGCTCCGGTGATTCATTCCCCCCAACACCCGCCTTTTT‑ACCCAACCCCAGGTG‑AAAAAAGC‑CTTT‑GGGGG‑GGGTTGTCGGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                           <  2:726091/150‑1 (MQ=255)
                                             gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                              >  3:141947/1‑146 (MQ=255)
                                                  gtgtCAGCGGTTACCGTGTTCCATTCCCCCCAACCGCCGCGCTTTT‑ACCCAACCCAATGGC‑CAAAATTG‑CCTC‑CGCGG‑CGCCTGCCTGCCAGTTTTTTTT‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑aa                                                                                                     <  2:1229653/149‑1 (MQ=255)
                                                   tgtCGGCGCCTACCGTATTTCTTTCCCCCCAACACCCCCCCTTTT‑ACCCACCCCCTTGTG‑‑AAAAATG‑CCTT‑CGGGGCCGGTTTTCTCCCGTTTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAAt                                                                                                   <  4:592157/150‑1 (MQ=255)
                                                     tCAGCGGCTCCCGTGATTCATTCCCCCCAACACCCGCGCTTTT‑ACCCACCGCCATGTG‑AAAAAATG‑CCTT‑CGGGG‑CGGTTTTCGGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAg                                                                                                 <  3:442843/150‑1 (MQ=255)
                                                      cAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggcggggaaaa                                                                                                   >  1:719173/1‑137 (MQ=255)
                                                       aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggtagtgagacaac                                                                                                 >  1:1375636/1‑136 (MQ=255)
                                                        gCGGCTCCGGTTATTCTTTCCCCCCAACACCCGCGCTTTT‑ACCCACCGCCATGTG‑AAAAAATG‑CCTT‑CGGGG‑CGG‑TTTTTTCCAGCTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAg                                                                                                 <  2:518525/147‑1 (MQ=255)
                                                          ggCTCCGGTGATTCTTTCCCCCAAACACCCGCGCTTTT‑ACCCACCGCAATGTG‑CAAAAATG‑CTTC‑CGGGG‑CGGCTGTCGGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑tt                                                                                             <  2:1729899/148‑1 (MQ=255)
                                                          ggCTCCCGTGATTCATTCCCCCAAACACCCGCGCATTT‑ACCCAACGCCATGTA‑AAAAAATC‑CCTT‑CGGGG‑CGGCTGTCTCCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAAt                                                                                                   <  2:1597961/143‑1 (MQ=255)
                                                           gCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGACGTTCTCTCCAGTAGTGGCggtg                                                                                                              >  1:952332/1‑130 (MQ=255)
                                                             accggTGATTCTTTCCCGCCAACACCCGCGCATTT‑ACCCACCGCCATGTG‑AAAAATGG‑CCTT‑CGGGG‑GGGTTGTCGCCCAGCTTTTTTC‑CCCCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGc                                                                                        <  2:1388773/146‑1 (MQ=255)
                                                               ccgtatttCATCCCCCCCAACACCCGCGCATTT‑ACCCACCGCCATGTG‑AAAAAATC‑CCTC‑CGCGG‑CGGCTGTCTGCCGGTTATTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCg                                                                                     <  4:50371/145‑1 (MQ=255)
                                                                cGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑GAATAAC‑TTTTAAccccc                                                                                    >  1:1201560/1‑149 (MQ=25)
                                                                     ttCTTTTCCGCAAACACCCGCGCTTTT‑ACCCAACGCCTTGTG‑AAAAAATC‑CCTT‑GGCGG‑CGGTTGTTTGCCTTCTATTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑ttt                                                                                            <  3:506470/138‑1 (MQ=255)
                                                                     ttCTTTCCCCCCAACACCCGCCTATTT‑ATCCACCGCCATGTG‑AAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCAt                                                                                 <  3:226223/149‑1 (MQ=255)
                                                                       catccccccaaacaccCGCACTTTT‑ACCCACCGCC‑TGTT‑AAAAAAAT‑CCTTGCGCGG‑CGGCTGTCGGCCATTTATTTTC‑CCGCCACAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCAtg                                                                                <  1:137982/138‑1 (MQ=255)
                                                                           tCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAgaatttaaaccgccc                                                                                   >  1:202671/1‑128 (MQ=25)
                                                                            ccccaaacaccCGCTCTTTT‑ACCCACCCCAATGTG‑AAAAAAGG‑CCTT‑CGCGG‑GGGCTGTCGGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTgg                                                                        <  2:644180/146‑1 (MQ=255)
                                                                             ccccAACACCCGCGATTTT‑ACCCACCGCCATGTG‑AAAAAATG‑CTTT‑CGCGG‑CGGTTGTCTGCCAGTTATTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGc                                                                       <  3:174882/149‑1 (MQ=255)
                                                                                      ccgctatttt‑ACCCACCCCCAGGTA‑AAAAAGTC‑CCTC‑CGGGG‑CGGCTTTCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgctgc                                                              <  2:1539586/144‑1 (MQ=255)
                                                                                          gctttt‑ACCAAACCCCTTGTA‑AAAAATTC‑CCTG‑CGCGG‑GGGTGGTCTGCAGGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                             <  2:580097/145‑1 (MQ=255)
                                                                                             ttt‑ATCCACCGCCATGTG‑CAAAAATC‑CCTC‑CGGGG‑CGCCTTTCTCCCAGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTaa                                                       <  2:456331/150‑1 (MQ=255)
                                                                                             ttt‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAAA‑TTTTAAACCCTCCTCTTGTGTgg                                                                        >  2:1309664/1‑134 (MQ=25)
                                                                                             ttt‑ACCCAACGCCATGTG‑AAAAATTG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTaa                                                       <  2:1041295/150‑1 (MQ=255)
                                                                                                 aTCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAG‑AAAAAAA‑TTTAAAACCCGCCTTTTGTGTGGGGAACGCTGCT‑ggggaaaa                                                    >  2:711177/1‑144 (MQ=21)
                                                                                                       cGCCTTGTA‑AAAAAGTG‑CCTC‑CGGGG‑CGGCTGTCTGCCGTCTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGa                                              <  4:382090/150‑1 (MQ=255)
                                                                                                         ccATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAG‑AAAAAAA‑TTTAAAACCCCCCTGTTTTGTGGGGATCGCCGGT‑GGGGACaa                                                    >  2:1225633/1‑142 (MQ=25)
                                                                                                             gtg‑AAAAAATG‑CCTC‑GGGGG‑CGGCTGTCTGCCAGTTATTTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGa                                        <  2:772815/150‑1 (MQ=255)
                                                                                                                  aaaaaTG‑CTTT‑CGCGG‑GGGGTCGCTCCCGGCTTTTTTC‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAg                                           <  1:1079752/143‑1 (MQ=255)
                                                                                                                  aaaaaGG‑CTTT‑CGGGG‑CGGCTGTCTGCCAGCTATTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTg                                    <  3:72326/150‑1 (MQ=255)
                                                                                                                      aTG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAG‑AAAAAACATTTAAACCC‑TCCTGGTGTGTGGGGAGCGCTCCTGGGTACAAATGGGGAAAGGGagggggg                                  >  3:413335/1‑143 (MQ=255)
                                                                                                                                  gg‑cggcTTTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTcc                     <  3:576118/151‑1 (MQ=255)
                                                                                                                                      ggcTGTCTGCCAGCTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                    <  2:1495079/149‑1 (MQ=255)
                                                                                                                                           tCTGCGAGTTTTTTCC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa                 <  1:1642425/147‑1 (MQ=255)
                                                                                                                                             cccctgtttttttc‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCgcgc           <  1:1426048/142‑1 (MQ=255)
                                                                                                                                              gCCAGGTATTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc               <  4:413561/146‑1 (MQ=255)
                                                                                                                                               ccattttttttc‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt          <  2:460283/143‑1 (MQ=255)
                                                                                                                                               ccAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAG‑GAAAAAC‑TTTTAAACCCCCCTGTTGTGTGGGAATCGCCACT‑GGGAAAAAAGTGGGAAGGGGAGGGGGGGGGAACAACCCCAGaa                 >  2:1593069/1‑143 (MQ=255)
                                                                                                                                                cggtttttttc‑CCGCAGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCgcgc           <  1:616437/142‑1 (MQ=255)
                                                                                                                                                 aGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAAC‑TTTAA‑CCCCTCCTTTTGTTTTGGGATCCCCGCT‑GGGTACAAATTGGGAAGGGGAGAGCGGGGGAACAGcccccgca                 >  2:1691015/1‑137 (MQ=255)
                                                                                                                                                  gttATTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt      <  2:362924/149‑1 (MQ=255)
                                                                                                                                                      ttttt‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTTGCAGAG‑AAAAAAA‑ATTAAAACCCTCCTTTTTTTTTGGAATAGCTGCT‑GGTAAAAATATGGAAAAGGGAGGGGGGGGGAGCAGCCCCCCAACCGCGCCAAATGAGt  >  4:396701/1‑151 (MQ=255)
                                                                                                                                                      ttttc‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg   <  2:202702/150‑1 (MQ=255)
                                                                                                                                                       tttc‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt  <  2:1538440/150‑1 (MQ=255)
                                                                                                                                                          |                                                                                                                                                      
CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  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 < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑