Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,890	0	C	A	56.6%	14.7 / 67.4	53	L267F (TTG→TTT)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (17/6);  new base A (0/30);  total (17/36)
Fisher's exact test for biased strand distribution p-value = 3.12e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-07

GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGAC  >  NC_000913/357742‑358039                                                                                                                                                         |                                                                                                                                                        gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaac                                                                                                                                                         <  1:1232236/149‑1 (MQ=255)  aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑aaa                                                                                                                                                     <  1:1301364/151‑1 (MQ=255)    aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑aaaa                                                                                                                                                    <  1:1205811/150‑1 (MQ=255)    aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAAc                                                                                                                                                   >  1:1303713/1‑151 (MQ=255)       cGTAGCGCATCAGGCAATTCCACCCGCTGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAAcc                                                                                                                                                  <  2:1219159/149‑1 (MQ=255)             gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTtc                                                                                                                                           >  2:554059/1‑150 (MQ=255)                tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑TCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTtctctc                                                                                                                                       >  1:1057746/1‑151 (MQ=255)                  aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTcc                                                                                                                                      >  1:1185853/1‑150 (MQ=255)                        ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCCCTCCAGGAGTg                                                                                                                               >  2:674692/1‑151 (MQ=255)                          ccACCCGCAGCCTTGTGTCGGCGGCTCCCGTGTTTCTTTCCCGCCAACAC‑CCGCGCTTTTACCAACCGCCTGGCGAAAAAATCCCTC‑C‑GGGGGGGTTGTCTGCCAGTTTTTTTC‑CCGCCCCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                             <  1:1324225/151‑1 (MQ=255)                           cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑g                                                                                                                           >  2:715539/1‑151 (MQ=255)                              ccGCAGCCTTGTCTCGGCGCCTCCTGTTATTTATTCCCGCAAACAC‑CCGCGCATTAATCCAACGCCTTGTG‑AAAAATGTCTTCC‑GCGGCGCCTGTCTGCAATTTTTTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCagtag                                                                                                                                 <  2:171586/143‑1 (MQ=255)                               cccagccatgtgCCGGCGGCTGCCGTGATCCATCCCCCCCAACAC‑CCGCGCTTTTACCCACCCCGATGTAAAAAAGTTCGTC‑G‑GCGGCGCCTGCCGCCCCGCTATTTTC‑CCCGCCAAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑ggg                                                                                                                         <  2:1456042/141‑1 (MQ=255)                               cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                                                        >  1:547115/1‑150 (MQ=255)                                   gcccTGTGTCAGCGCCTACCGTTTTTCTTTCCCGCCAACAC‑CCGCGCTTTTACCCACCGCCATGTAAAAAAATGCCTT‑C‑GCGGCGGCTTTCTGCCAGCTTTTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAag                                                                                                                    <  1:657477/148‑1 (MQ=255)                                      cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                                                        >  1:585558/1‑143 (MQ=255)                                         tttgggCGTCTACCGTGATTCTTCCCCCCCAACAC‑CCGCGCTTTTATCCAGCCCCATGGGAAAAAGTCCCTC‑C‑GCGGCGGCTGTGTCCCAGTTATTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGaaa                                                                                                               <  1:664317/144‑1 (MQ=255)                                         tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑ggggggagggggaa                                                                                                              >  2:692297/1‑140 (MQ=255)                                              gcgcCTCCCGTGTTTCATTCCCCCCACCAC‑CCGCGCTTTTCCCCACCCCTTTGTAAAAAAAT‑CCTT‑CGGCGGCGGCTGTCTCCCAGTTTTTTTC‑CCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAg                                                                                                            <  2:98300/147‑1 (MQ=255)                                                ggCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGaaaaaaattt                                                                                                        >  1:858129/1‑142 (MQ=255)                                                ggCTACCGTAATTCTTTCCCGCCAACAC‑CCGGGCTTTTACCCACCCCCCTGTAAAAAAACGCCTG‑C‑GCGGGGGCTGCCTCCCAGTTTTTTTC‑CCGCCCAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAAt                                                                                                              <  1:15446/143‑1 (MQ=255)                                                gcctcccGTGTTTCTTTCCCGCCACCAC‑CCCCGCTTTTACCAACCGCCAGGT‑CAAAAAAGCCTTCC‑GCGGCGGTTGCCTGCCAGTTTTTTTC‑CCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGaa                                                                                                                <  2:230456/136‑1 (MQ=255)                                                gcctcccGTGATTAATTCCCGCAACCAC‑CCGCGCTTTTACCCACCCCCATGTGAAAAAAGCCCTT‑C‑GCGGCGGCTTTCTGCCAGTTTTTTTC‑CCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCt                                                                                                          <  1:949156/142‑1 (MQ=255)                                                 gCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGAGgaaaaaattt                                                                                                        >  1:181921/1‑141 (MQ=255)                                                  cccccGTGTTTCTTTCCCGCCAACAC‑CCCCGATTTTACCCACCGCTATGTCAAAAAATCCTTC‑C‑GCGGCGGCTTTCTGCCAGTTTTTTTC‑CCGCCACAAA‑AAATCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAg                                                                                                     <  2:1217260/147‑1 (MQ=255)                                                  cTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAGAAGaaaaaaactt                                                                                                        >  2:1323704/1‑140 (MQ=25)                                                    ccggttaTTAATTCCCGCCAACAC‑CCCCCCATTTACCCACCCCATTGTGAAAAAATCCCTC‑G‑GCGGCGGCTTTCGGCCAGTTTTTTTC‑CCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGccc                                                                                                  <  2:1346563/145‑1 (MQ=255)                                                    aCCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAGAGAAAAAACTTTAAAccc                                                                                                  >  1:598663/1‑151 (MQ=255)                                                      cggtttttaTTTCCCCCAACAACCCGCGCATTTACCCACCCCCTTGTAAAAAATTCCCTC‑C‑GGGGCGGTTGTCGGCCAGTTTTTTTC‑CCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCtt                                                                                                         <  2:1254825/135‑1 (MQ=255)                                                               ttccccccAAAAC‑CCGCCCTTTTACCCACCCCAATGTAAAAAAATCCCTG‑C‑GCGGGGGCTGTCTGCCAGTTTTTTTC‑CCCCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                                                                                           <  1:592878/141‑1 (MQ=255)                                                               tccccGCCACCAC‑CCGCTTTTTTACCCCACGCCATGTAAAAAAATCCTTC‑C‑GCGGGGGCTGTCTGCCAGTTTTTTTC‑CCGCCACACA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGACCGGCAtgttg                                                                                         <  2:608206/146‑1 (MQ=255)                                                                  cccccAACAC‑CCGCGCTTTTACCCACCGCCTTGTGAAAAATTGCTTT‑G‑GGGGCGGCTTTCTGCCAGTTTTTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgt                                                                                      <  1:252997/146‑1 (MQ=255)                                                                   ccccaaccacccccgcTTTTACCCACCCCCATGGGAAAAATTCCCTG‑C‑GGGGCGGTTGTC‑GCCAATTTTTTTTCCCGCCACAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                                                                                           <  2:954952/131‑1 (MQ=255)                                                                       aacac‑cCGCGCATTTACCCACCCCCTTGTGAAAAAATTCCTT‑C‑GGGGCGGCTTTCTGCCATTTTTTTTC‑CCCCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt                                                                                <  1:1368108/150‑1 (MQ=255)                                                                             ccgcacttttccccACCGCCATGTAAAAAAGTGCATG‑G‑GCGGGGGCGGTCGGCGGTTTTTTTGC‑CCGCCACAAA‑AAAGCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                               <  2:688549/134‑1 (MQ=255)                                                                               gcgcATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAAAAAAAAACTTTAA‑CCCCTCCTGTTTTGTTGGAATccc                                                                             >  2:1309730/1‑144 (MQ=25)                                                                                       aTCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGAGAAAAAAATTTTAACCCCTCCATTTGTGTTGCGAACGCTGCTggggaaaa                                                                 >  2:774083/1‑144 (MQ=25)                                                                                         ccaccgccagggaaaaaaaTGCCTG‑C‑GCGGCGCCTGTCTCCCAGTTTTTTTC‑CCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                                   <  2:1276166/132‑1 (MQ=255)                                                                                             ccccttgtaaaaaaaTGCTTC‑C‑GCGGCGGTTGCCGGCCAGTTTTTTTC‑CCCCCGAAAA‑AAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                                               <  2:852803/136‑1 (MQ=255)                                                                                               ccagggaaaaaaaTGCCTC‑C‑GGGGCGGTTTTCGGCCAGTTTTTTTC‑CCCCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                                           <  2:751548/140‑1 (MQ=255)                                                                                                   gtcaaaaaaTCCCTG‑C‑GGGGCGGTTGTCTCCCATTTTTTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                                       <  1:62683/144‑1 (MQ=255)                                                                                                   gtaaaaaaaGGCCTC‑C‑GCGGGGCTTTTCTGCCAGTTTTTTTC‑CCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                                       <  2:75230/144‑1 (MQ=255)                                                                                                    taaaaaaaTCCTTC‑C‑GCGGCGCTTGTCTGCCAGTTTTTTTC‑CCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                                      <  2:476208/145‑1 (MQ=255)                                                                                                    taaaaaaaTCCCTC‑C‑GCGGCGGCTTTCTCCCATTTTTTTTT‑CCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATAGCTGCTGGTAACAATGTGGAAAGGGa                                                     <  1:97284/146‑1 (MQ=255)                                                                                                    aacAAAAATGCCTC‑C‑GCGGCGGTTGTCTGCCGGTTATTTTC‑CCGCCACAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                                   <  2:1143154/149‑1 (MQ=255)                                                                                                         aaaTGCCTG‑C‑GGGGCGGCTGTCTGCCAGCTATTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                                              <  2:588244/151‑1 (MQ=255)                                                                                                          aaTGCCTT‑CGGGGGCGCCTTTCTCCCGTCT‑TTTTT‑CCCCCCAAACAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                                              <  2:1513976/151‑1 (MQ=255)                                                                                                               cTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGAAAAAAATTTTAAACCCCCCTGTTTTGTGGGGATTGCTGCTGGGAACAAAGGGGGAAGGGGAgggggggg                                            >  2:1343207/1‑140 (MQ=255)                                                                                                                       ggcggcTGTCTGCCAGTTTTTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc                                   <  1:437978/150‑1 (MQ=255)                                                                                                                       ggcggcTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGACGAGGAAAAGAATATTACCCGACATATTGTGTGGGGCTCGGCGCTGGGGACAATATTGAAAAGGGAGGGCGGGGAGTCAGc                                     >  2:166503/1‑148 (MQ=255)                                                                                                                                     atcttttttc‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa                     <  1:190109/145‑1 (MQ=255)                                                                                                                                        ttttttc‑CCGCCACAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg                  <  2:165898/148‑1 (MQ=255)                                                                                                                                             ttcccgccgAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt              <  1:1205206/147‑1 (MQ=255)                                                                                                                                                 cgccgaAAC‑AAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCt         <  1:1103636/151‑1 (MQ=255)                                                                                                                                                        ac‑aaAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGac  <  1:585739/151‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                        GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGAC  >  NC_000913/357742‑358039

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