Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,890	0	C	A	45.6%	45.1 / 41.1	69	L267F (TTG→TTT)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (24/13);  new base A (0/31);  total (24/45)
Fisher's exact test for biased strand distribution p-value = 2.57e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC  >  NC_000913/357741‑358033                                                                                                                                                            |                                                                                                                                                   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)                                                                                                                                      ccaagttttttgcccgccgCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa                        <  1:1558385/131‑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)                                                                                                                                          gttttttgcccgccaCAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt                 <  2:1706919/138‑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)                                                                                                                                                ttcccgccgCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt        <  3:110444/147‑1 (MQ=255)                                                                                                                                                   cccccgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg      <  3:440625/146‑1 (MQ=255)                                                                                                                                                     gccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGCGAG‑AAAGAAC‑ATTAAGAACCTCACGCTGTGGGGGAATTCCTGCT‑TGGAACACTTTGGTAAGGGGAGGGGCGCGAATCAACTCCCCAACCCCGCCAAATTAGtttccccc  >  3:34247/1‑145 (MQ=255)                                                                                                                                                        acaaaaaACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGTCGAATCAGCTCCAGAACCGCGCTAATTGAGtt        <  1:1349779/141‑1 (MQ=255)                                                                                                                                                          aacaaACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGcc    <  1:51732/144‑1 (MQ=255)                                                                                                                                                            |                                                                                                                                                   TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTT‑ATCCAACGCCATGTG‑CAAAAATG‑CCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC  >  NC_000913/357741‑358033

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