Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,890	0	C	A	40.7%	53.9 / 27.0	59	L267F (TTG→TTT)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (20/15);  new base A (0/24);  total (20/39)
Fisher's exact test for biased strand distribution p-value = 2.00e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATG  >  NC_000913/357741‑358037                                                                                                                                                       |                                                                                                                                                     tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                                      >  1:506113/1‑150 (MQ=255) tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                                      >  1:921270/1‑150 (MQ=255)     aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                                 <  1:1359615/151‑1 (MQ=255)          tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                                                                            >  2:935314/1‑151 (MQ=255)           aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                                                                           >  2:1326/1‑151 (MQ=255)                  cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                                    >  2:523498/1‑151 (MQ=255)                     gaattcccACCCGCAGACCTGTGTCAGCGGCTCCCGTGATTCATTCCCGCCAACACCCGCGCTTTTATCCAACGCCATGTGCAAAAATGCCTCC‑GGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                                 <  1:400853/145‑1 (MQ=18)                     caaaTCCCACCCGCAGACTTTTGTCAGCGGCTCCCGTGATTCTTTCCCCCCAACACCCGCGCATTTATCCACCGCCATGTCAAAAAATGCCTTC‑GCGGCGGCTGTCTGCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagt                                                                                                                                 <  1:845604/149‑1 (MQ=255)                        aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                              >  1:869322/1‑151 (MQ=255)                         ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCTGTAGTg                                                                                                                             >  1:215660/1‑151 (MQ=255)                          cccAGCCGCAGCCCTGTGCCAGCGGCTACCGTGATTAATTCCCGCCAACACCCGCGCATTTACCCAACGCCTTGTAAAAAAATGCCTCC‑GCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                                            <  1:205169/150‑1 (MQ=255)                           ccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                           >  2:1276439/1‑151 (MQ=255)                                cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                                      >  2:438042/1‑150 (MQ=255)                                 gCAGCCCTGTGTCAGCGGCTCCCGTGATTCATTCCCCCAAACACCCGCGCATTTACCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                                     <  2:497215/151‑1 (MQ=21)                                 gCAGACTTGTGTCAGCGCCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCACCGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                                     <  2:306067/151‑1 (MQ=21)                                   aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggac                                                                                                                    >  2:190467/1‑147 (MQ=255)                                    gcccTGTGTCAGCGGCTACCGTGATTCATTTCCGACAACAACCGCGCATTTCCCCAACGCCATGTACAAAAATGCCTTC‑GGGGCGGCTGTCTCCCGTTTATTTTTCCGCCGCAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                                  <  2:1096285/149‑1 (MQ=21)                                     cccTGTGTCAGCGGCTCCCGTGTTTCATTCCCGCCAACACCCGCGCTTTTACCCACCCCCATGTGAAAAAATGCCTTC‑GGGGGGGCTGTCTGCCAGTTATTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                                 <  1:721587/150‑1 (MQ=255)                                      ccTGTGTCAGCGACTCCCGTGATTCATTTCCGCCAACACCCGCCTATTTACCCACCGCCATGTGAAAAAATCCCTGC‑GCGGCGCCTGTCTGCCAGTTATTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgag‑a                                                                                                               <  2:193390/151‑1 (MQ=255)                                         gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGTGAGAG‑aaa                                                                                                             >  1:175550/1‑150 (MQ=255)                                         gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑aaa                                                                                                             >  2:1262651/1‑150 (MQ=255)                                           gtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGAa‑aaaaa                                                                                                           >  2:582414/1‑148 (MQ=255)                                            ccggCGGCTCCGGTGTTTCATTTCCCCCAACACCCGCGCATTTCCCCAACCCCATTTGAAAAAAGCCCCCG‑GCGGGGGCTCTCTCCCGCTTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGc                                                                                                         <  1:260247/148‑1 (MQ=255)                                                cGGCTACCGGGTTTATTTCCCCCCAACACCCGCGCATTTACCAACCGCCATGTCAAAAAATGCCTTC‑GCGGCGCTTGTCGGCAAGCTATTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCtttn                                                                                                     <  1:313007/151‑2 (MQ=255)                                                  gCTCCCGTGATTCTTTTCCCCCAACACCCGCGCATTTACCCAACCCCATGGGAAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAg                                                                                                   <  2:326654/151‑1 (MQ=255)                                                    tccggTGATCCATTCCCGCAAACACCCGCGCTTTTCCCCAACG‑CATGGTGAAAAAAGCCTTCGGGGGCGGTGGCTTGCCAGCTATTTTTCCGCGAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAg                                                                                                   <  1:1273520/145‑1 (MQ=255)                                                     aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑gaaaaaacattaaacc                                                                                                >  1:731275/1‑138 (MQ=25)                                                     aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑aaaaaaattttaaacc                                                                                                >  2:108562/1‑138 (MQ=25)                                                      cggTGATTCTTTCCCGCCAACACCCGCGCTTTTACCCACCGCCATGTAAAAAAATGCCTCC‑GCGGCGGCTGTCGCCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCg                                                                                               <  1:131205/149‑1 (MQ=255)                                                      ccGTGATTCATTCCCCCCAACACCCGCGCATTTACCCAACGCCATGTCCAAAAATGCTTTC‑GCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCg                                                                                               <  2:755731/151‑1 (MQ=255)                                                        gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑GAAAAAATTTaaacgcggc                                                                                             >  1:280143/1‑144 (MQ=25)                                                          gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAACTTTTAAccccccat                                                                                           >  2:538333/1‑146 (MQ=25)                                                           attcattcCCGCCAACACCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtg                                                                                          <  2:501814/151‑1 (MQ=255)                                                           attcattcCCCCCAACACCCGCGCTTTTATCCAACCCCATGTGAAAAAATGCCTCC‑GGGGCGGCTGTCGGCCAGCTATTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtg                                                                                          <  2:1322341/151‑1 (MQ=255)                                                            ttCATTCCCCCCAACACCCGCGCTTTTACCCACCCCCTTGTAAAAAAATGCCTTC‑GCGGCGGCTGTCTGCAAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtg                                                                                          <  1:762902/150‑1 (MQ=255)                                                                    cGCCAACACCCGCGCTTTTACCCAACGCCATGTGAAAAAAGGCCTTC‑GGGGCGGCTGTCTCCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGc                                                                                 <  2:758031/151‑1 (MQ=255)                                                                         ccacccGCGCATTTCCCCACCCCGATGCGAAAAAATCCCTCC‑GCGGCGCTTCTCTGCCAGCTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                                                            <  1:396006/147‑1 (MQ=255)                                                                         acaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGAG‑AAAAAGAGTTATACCCGCCATGTTTGGTGGGAATcc                                                                            >  1:186325/1‑150 (MQ=25)                                                                             ccGCTCTTTTATCCAACGCCAGGTGAAAAAATGCCTCC‑GCGGCGGCTGTCTGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgc                                                                        <  1:779763/151‑1 (MQ=255)                                                                                     ttATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTTAAACCCG‑CCTGTTTTGTGGGAAACGCCGCTgggtac                                                                 >  2:705248/1‑146 (MQ=25)                                                                                       aTCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAACTTTAAACCCCTCCTGTTTTGTGGGGATCGCTGCTGGTAaaaa                                                               >  2:1225082/1‑150 (MQ=37)                                                                                        acccacccccgtatgaaaaaaGCCTTCGGCGGCGGCTGTCGGCCGGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                              <  1:1195386/136‑1 (MQ=255)                                                                                           aaCGCCATGTGAAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                                          <  2:128052/151‑1 (MQ=255)                                                                                            ccGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                                         <  2:516676/150‑1 (MQ=255)                                                                                                 aTGTGAAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGTTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                                    <  2:512904/151‑1 (MQ=255)                                                                                                  tgtgCAAAAATGCCTGC‑GCGGCGGCTGTCTGCCGTCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                                     <  2:255303/149‑1 (MQ=255)                                                                                                       aaaaaTGCCTTC‑GCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                                              <  2:134841/151‑1 (MQ=255)                                                                                                        aaaaTGCCTTC‑GCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg                                             <  2:506793/151‑1 (MQ=255)                                                                                                         aacggCCTTC‑GGGGCGGCTGTTGGCCAGTTATTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                                            <  2:281469/147‑1 (MQ=255)                                                                                                             gCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑GAAAAAATTTTAACCCGCCCTTTTTTGTTGGGAACGCTGCGGGGGAAAAAGTGGGAAGGGGAgggggggggg                                        >  1:1272346/1‑142 (MQ=255)                                                                                                                ttC‑GGGGCGGTTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                                       <  1:864728/149‑1 (MQ=255)                                                                                                                  ccggggcggcTGTCTGCCGGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc                                      <  1:119666/145‑1 (MQ=255)                                                                                                                        cggcTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg                              <  2:599407/151‑1 (MQ=255)                                                                                                                        cggcTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg                              <  2:295787/151‑1 (MQ=255)                                                                                                                          gcTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa                            <  1:183293/151‑1 (MQ=255)                                                                                                                            tGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAG‑AAAAAGATTTAAGCCCCCCCTGTTTTGTTGGGATCGCTGCTGGGGAAAATGTGGGAAGGGGAGGGGGGGGGAGCAGATCCAGaaac                          >  2:861969/1‑149 (MQ=255)                                                                                                                             gtcggcaggttATTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg                         <  1:1276440/141‑1 (MQ=255)                                                                                                                             gTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                          <  2:588428/150‑1 (MQ=255)                                                                                                                                    agttATTTTCCCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt                  <  1:1262652/148‑1 (MQ=255)                                                                                                                                        tttttcCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa               <  2:921274/149‑1 (MQ=255)                                                                                                                                        tttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa               <  1:935310/149‑1 (MQ=255)                                                                                                                                                    caacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATg  <  2:175544/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                     TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATG  >  NC_000913/357741‑358037

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