Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,880	0	T	C	50.0%	58.6 / 23.2	58	K271E (AAA→GAA)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (22/7);  new base C (0/29);  total (22/36)
Fisher's exact test for biased strand distribution p-value = 5.55e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAG  >  NC_000913/357733‑358025                                                                                                                                                      |                                                                                                                                                      cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                                    <  2:1168753/151‑1 (MQ=255)    aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                                 >  2:673226/1‑151 (MQ=255)    aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                                  <  1:1095292/150‑1 (MQ=255)       ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaa                                                                                                                                              <  2:1230482/151‑1 (MQ=255)             aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAAc                                                                                                                                       <  1:192832/151‑1 (MQ=255)             aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAAc                                                                                                                                       >  2:803513/1‑151 (MQ=255)              aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAAcc                                                                                                                                      >  2:1002955/1‑151 (MQ=255)                cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAAccg                                                                                                                                     <  2:1088880/150‑1 (MQ=255)                 gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAAccgcc                                                                                                                                   >  1:907089/1‑151 (MQ=255)                 gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAAccgcc                                                                                                                                   <  2:905316/151‑1 (MQ=255)                 gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAAccgc                                                                                                                                    <  1:517911/150‑1 (MQ=255)                        aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTtctct                                                                                                                            >  1:1042438/1‑151 (MQ=255)                           aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTcc                                                                                                                          >  2:1115605/1‑150 (MQ=255)                             gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCagt                                                                                                                       >  1:492135/1‑151 (MQ=255)                                       ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                             >  2:1120553/1‑151 (MQ=255)                                       ccGAAGCCCTGTGTCAGGGCCTCCCGTGATCCTTTCCCGCCACCACCCGCGCTTTTACCCAACGCC‑TTGTGCAAAAATGCT‑TTCGCGGCGCCTGTCTGCCAGTTTTTTTCCCGCCCCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                <  1:346115/148‑1 (MQ=255)                                          cagccttgtgTCGGCGGTTACCGTGATTCATTCCCCCCAACCCCCGCCCTTTTACCCACCGCG‑TTTTAAAAAATTCCC‑TGCGCGGGGGCTGTCTGCCAGCTTTTTTCCCCCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                           <  1:780243/144‑1 (MQ=255)                                            gccttttgtCAGGGGCTCCCGTTTTCCTTTCCCCCCAACCCCCGCGCTTTTACCCACCGCA‑TTGTGAAAAATTCCC‑TCCGGGGGGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                          <  2:189059/143‑1 (MQ=255)                                             cccTGTGTCGGCGCCTCCCGTGATTCATTCCCCCCAACACCCGCCAATTTACCAACCCCC‑ATGTGAAAAAATGCT‑TTGGGGGCGGCTTTCTGCCGGCTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga‑g                                                                                                      <  1:970074/150‑1 (MQ=255)                                                   gtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGGGA‑G‑aaa                                                                                                  >  1:1163795/1‑148 (MQ=255)                                                    tCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCgggggcaaaggaaaa                                                                                                  >  2:70114/1‑138 (MQ=255)                                                     cagcgccccccGTGTTTCTTCCCCCCCAACACCCGCGCTTTTACCCACCGCC‑ATGTGAAAAAAGCCC‑TTCGGGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                          <  1:860479/131‑1 (MQ=255)                                                         gcctcccGTGATTCATTTCCCCCAACAACCGCGCTTTTACCCACCGCC‑ATGTGAAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGAAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTaa                                                                                         <  2:343693/146‑1 (MQ=255)                                                            tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑CTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGA‑G‑aaaaaaattttaaac                                                                                      >  2:756624/1‑139 (MQ=25)                                                             cccgtttttCATTCCCCCCAACACCCGCGCATTTACCCACCCCC‑TTGTGAAAAAATGCC‑TCCGCGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGc                                                                                              <  1:634222/135‑1 (MQ=255)                                                              ccGTGTTTCATTCCCCCAAACAACCGCGCATTTACCCAACGCC‑TTGT‑AAAAAACGCCGTTCGCGGGGGGTCTCTGCCGCCTATTTTCCCGCCACAAA‑AAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGcc                                                                                      <  2:486866/149‑1 (MQ=255)                                                              ccGTGATTCTTTCCCCCCAACACCCCCGCTTTTACCCACCCCC‑ATGTGAAAAAATCCC‑TTCGGGGGGGCTGTCTGCCAGTTTTTTTCCCCCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGcc                                                                                      <  2:1147132/149‑1 (MQ=255)                                                               cGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑AAAAAACTTTTAAccc                                                                                     >  2:1011073/1‑149 (MQ=25)                                                                 taaTTCTTTCCCCCCAACACCCGCGCTTTTATCCAACGCC‑ATGTAAAAAAATGCC‑TCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGt                                                                                   <  1:66639/147‑1 (MQ=255)                                                                    ttCATTCCCGCCAACACCCGCGTTTTTA‑CCAACCCCAATGTGAAAAAAGCCT‑TTCGCGGCGTCTTTGTGCAAGCTTTTTTCCCCCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtg                                                                               <  2:1145357/150‑1 (MQ=255)                                                                      cATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGG‑G‑AAAAAGAATTAAGacccccagttcg                                                                            >  1:155178/1‑142 (MQ=25)                                                                      cATTCCCCCCAACACCCGCGCATTTACCAACCGCC‑ATGTGCAAAAATCCT‑TTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtgt                                                                              <  2:895164/149‑1 (MQ=255)                                                                        ttcccccaaacaccCGCGCTTTTACCCAACGCC‑ATGTGCAAAAATCCC‑TCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtgttgt                                                                           <  2:283485/142‑1 (MQ=255)                                                                        ttcccGCCAACACCCGCGCTTTTACCCACCGCC‑ATGTGAAAAAATGCC‑TTCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtgttg                                                                            <  2:1194130/146‑1 (MQ=255)                                                                          cccGCCAACACCCGCGTATTTACCCACCGCC‑TTGTAAAAAATTGCT‑TCCGGGGGGGTGGTCTGCCAGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtgttg                                                                            <  1:1120003/146‑1 (MQ=255)                                                                           ccccaaacaccCGCGCTTTTACCCAACGCC‑ATGTAAAAAAATGCT‑TGCGCGGCGGCTGTCTGCAAGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTtgt‑gtg                                                                       <  1:931570/145‑1 (MQ=255)                                                                            cGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAA‑G‑AAAAAAATTTTAACCCCCCCTGTTTT‑GTggg                                                                     >  2:1127865/1‑150 (MQ=25)                                                                            cGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAAAAGCTTTAAACCCGCCCTGTTTT‑GTgg                                                                      >  1:200887/1‑150 (MQ=255)                                                                              ccAACACCCGCGCATTTCCCCCACGCC‑ATTGGAAAAAATGCC‑TTCGCGGCGGCTTTTTGCCCGTTTTTTTCCCGCCCCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGT‑GTGGCa                                                                    <  2:1009480/150‑1 (MQ=255)                                                                                   aaccgcaattttACCAAACCCC‑ATGTCAAAAAATCCT‑TGCGGGGCGGTTTTCTCCCAGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGT‑GTGGCAATCgc                                                               <  1:836903/141‑1 (MQ=255)                                                                                      cgcttttttACCCAACCCC‑ATGTACAAAAATCCC‑TTCGCGGCGCCTTTCTGCCAGCTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGT‑GTGGCAATc                                                                 <  1:879879/139‑1 (MQ=255)                                                                                         gtaTTTACCCAACGCC‑ATGTGAAAAAGTGCC‑TGCGCGGCGGCTGTCTGCAAATTATTTTCCCCCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCTGGt                                                        <  2:907730/149‑1 (MQ=255)                                                                                                   aaCGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGG‑AAAAAATTTTAACCCCG‑CATGTTGT‑TTGGGGATCGCTGCgggggaaaaagtgg                                              >  2:349295/1‑140 (MQ=25)                                                                                                    aCGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGG‑A‑AAAAAAATTTATACCCCTCCTTTTTT‑GGGGGGATCGCTGCTGGGAaaaaaatg                                               >  1:379811/1‑145 (MQ=19)                                                                                                     cGCC‑ATTTCAAAAATTGCT‑TCCGCGGGGGCTTCCTGCCAGTTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCTGGTAACAATGTGGaa                                            <  2:1168304/151‑1 (MQ=255)                                                                                                          aTGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑GAAAAAAATTAAACCCCTCCTGTTTT‑GTGGGAATCGCCGCTGGtgaaaaaatgg                                              >  2:1095333/1‑135 (MQ=25)                                                                                                           tgtaCAAAAATGCT‑TTGGCGCCGGTTGTCTCCCTTCTATTTTCCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                         <  2:183160/149‑1 (MQ=255)                                                                                                            gtgAAAAAATGCT‑TGCGGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCACAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                        <  2:719930/149‑1 (MQ=255)                                                                                                                aaaatTGCT‑TGCGCGCCGGTTTTCTGCCATTTTTTTTCCCGCCGAAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                                  <  2:699656/151‑1 (MQ=255)                                                                                                                aaaaaTGCC‑TTCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                    <  2:445070/149‑1 (MQ=255)                                                                                                                  aaaTGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAAAAAGTTTTAATCCCG‑CCTTTTGTGGGGGCAAT‑GCCGCTGGTGACAAAGTGGGAAGGGGAgggggg                                >  2:245170/1‑146 (MQ=255)                                                                                                                   aaTGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑A‑AAAAAAATTTTAACCCCCCCTGTTGT‑GTGGGGAACGCTGCTGGGGACAAAGTGGGAAGGGGAgggggg                                >  2:289343/1‑145 (MQ=255)                                                                                                                   aaTCCC‑TCCGGGGGGGTTGTCGCCCAGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggc                               <  1:246686/151‑1 (MQ=255)                                                                                                                        ctttCGGGGCGGTTGTCTGCCAGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggc                               <  2:1066737/145‑1 (MQ=255)                                                                                                                               ggcggcTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGCGCGAGA‑GAAAAAAGTTTTAAACCCGCCATTTGTT‑TTGGCAACTCCGCCTGGTAAAAAATGGGAAGGGGAAGGGGGGGGAAtaactcc                     >  2:421487/1‑145 (MQ=255)                                                                                                                                  gccTGTCCCCCAGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGa                             <  1:596731/137‑1 (MQ=255)                                                                                                                                            ccattATTTTCCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCgcgc          <  1:585288/145‑1 (MQ=255)                                                                                                                                                ttttttcCCCCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt     <  2:609120/147‑1 (MQ=255)                                                                                                                                                tATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGTGCGAGG‑G‑AAAAAAATTTTAAACCCTCCTTTTTT‑TTGGGGATCGCTGCTGGGGACAATATGGGAAGGGGAGGGGGGGGGATCAACCCCCCAACCGCGCCAAtt     >  1:429274/1‑149 (MQ=255)                                                                                                                                                  ttttcCCCCCGAAAC‑AAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg  <  1:1267810/150‑1 (MQ=255)                                                                                                                                                    tttCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGCGA‑G‑AAAAAAATTTAATCCCCCCATGTTGT‑GTGGGGAACGCTGCTGGGGAAAATGTGGGGAAGGGAGGGGGGGGGAGCAGCTCCACAAcccccccca       >  2:994692/1‑136 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAG  >  NC_000913/357733‑358025

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