Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,880	0	T	C	50.9%	57.7 / 84.5	55	K271E (AAA→GAA)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (21/6);  new base C (0/28);  total (21/34)
Fisher's exact test for biased strand distribution p-value = 3.52e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-10

TAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAT‑GTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357735‑358026                                                                                                                                                   |                                                                                                                                                       tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                                   >  2:46980/1‑151 (MQ=255) tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                                   >  1:776044/1‑151 (MQ=255)  aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                                   <  1:500047/150‑1 (MQ=255)   ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                                                                                                                 >  2:299063/1‑151 (MQ=255)    gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                                                                                                                                <  1:391361/151‑1 (MQ=255)      cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                              >  2:919363/1‑151 (MQ=255)          tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                           >  1:951993/1‑150 (MQ=255)           aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                          >  1:132722/1‑150 (MQ=255)            aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGAGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                         >  2:884333/1‑150 (MQ=255)             gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                       <  1:1008604/151‑1 (MQ=255)                    gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct                                                                                                                                >  2:360138/1‑151 (MQ=255)                                 ccAGCCGCAGCCCTGTGTCAGCGGCTCCCGTGATTCATTTCCGCCAACACCCGCGCTTTTACCCACCGCCATGTGCAAAAATGCTTC‑CGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                     <  2:1006094/149‑1 (MQ=255)                                   aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                 >  1:758245/1‑151 (MQ=255)                                     ccccgcaCCTTTGACAGCGGTTCCGGTTATTCATTCCCCCAAACACCCCCACTTTTCTCCACCCCCATGTCAAAAAAGGCCTT‑CGCGGCGGCTGTCGGCCATTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                 <  1:115071/143‑1 (MQ=255)                                     ccccAGCCCTGTTTCCGCGGCTCCCGTGTTTAATTCCCGCAACCACCCGCGATTTTCTCCACCGCCATGGGAAAAAATCCCTT‑CGGGGGGGTTGTCTGCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                               <  1:57281/148‑1 (MQ=255)                                     ccGCAGCCCTGTGTCAGCGGCTCCGGTGATTCATTCCCGCCAACACCCGCGCTTTTACCCAACGCCATGTGAAAAAATGCCTC‑CGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagta                                                                                                                        <  2:40550/142‑1 (MQ=255)                                     ccGCAGCCCTGTGTCAGCGGCTCCCGTGTTTCATTCCCCCCAACACCCGCGCTTTTACCCAACCCCTTGTGAAAAAATGCCTT‑CGCGGGGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                  <  1:919383/148‑1 (MQ=255)                                     ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                               >  2:258910/1‑151 (MQ=255)                                        cAGCCCTGTGTCGGGGGCTCCCGTGATTCATTCCCGCAAACACCCCCGCTTTTACCCACCGCCTTGTGAAAAAATGCCTC‑CGGGGGGGCTGTCTCCCGTCTATTTTTCCGCGGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                 <  1:950404/146‑1 (MQ=255)                                         aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggc                                                                                                             >  2:8990/1‑147 (MQ=255)                                                  tCAGCGGCTCCCTTGATTCTTTCCCCCCAACACCCGCCCTTTTACCCAACCCCATGTGAAAAAATGCCTC‑CGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAg                                                                                                 <  2:33752/150‑1 (MQ=255)                                                  tCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑aaaaa                                                                                                  >  2:302734/1‑147 (MQ=255)                                                      cGGCTCCCGTGATTCATTTCCCCCAACACCCGCGCTTTTACCCACCGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCGTTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑ttt                                                                                            <  2:503977/150‑1 (MQ=25)                                                       ggCTCCCGTGATTCATTCCCGCCAACACCCGCGCATTTACCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTa                                                                                           <  2:237164/150‑1 (MQ=255)                                                         ctcccGTGATTCATTCCCCCCAACACCCGCGCTTTTACCCAACGCCATGTGCAAAAATGCCTT‑CGCGGGGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTaa                                                                                          <  1:823704/146‑1 (MQ=25)                                                         cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑AAAAAAC‑ttt                                                                                            >  2:628029/1‑147 (MQ=255)                                                          tcccGTGATTCATTCCCGCAAACACCCGCGCATTTACCCACCGCCATGTGCAAAAAGGCCTC‑CGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAg                                                                                         <  1:938377/147‑1 (MQ=255)                                                            ccGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑aaaaaaaatttaaac                                                                                        >  1:61807/1‑137 (MQ=25)                                                              gTGATTCTTTCCCGCCAACACCCGCGCATTTACCCAACGCCATGTGAAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGt                                                                                    <  2:637840/150‑1 (MQ=255)                                                                       ttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCATGACG‑AAAAAGC‑TTTaaaacatcatctt                                                                               >  1:161670/1‑136 (MQ=25)                                                                        cccGCCAACAACCGCGCATTTATCCAACGCCATGTGAAAAAATGCCTC‑CGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTtgtgtg                                                                         <  2:53662/150‑1 (MQ=255)                                                                         cccccacaacccGCGCATTAACCCAACGCATGGAGAAAAAATGCTTC‑CGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTtgtg                                                                           <  1:461832/138‑1 (MQ=255)                                                                             caacaaCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑G‑AAAAAAA‑TTTAAACCCGGCATGTTGTGTGGGaa                                                                     >  2:360363/1‑150 (MQ=37)                                                                               acaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGG‑G‑AAAAAAA‑TTTAAAACCCCCCTGTTGTGTgggga                                                                     >  2:92488/1‑145 (MQ=25)                                                                                  acCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCATAA‑GAAAAAAGC‑TTTAAACCC‑CCTTGTTGTGTGGcgaacc                                                                  >  2:47901/1‑143 (MQ=25)                                                                                   ccGCGCTTTTACCCACCGCCTTGTGAAAAAATCCTTC‑CGCGGCGGCTTTCCGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATc                                                                   <  1:779138/146‑1 (MQ=255)                                                                                   ccGCGCTTTTACCCACCCCCTTGTGAAAAAATGCCTT‑CGGGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgctgc                                                              <  1:280207/151‑1 (MQ=255)                                                                                       gcttttACCCAACCCCATGTCAAAAAATCCTTT‑CGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt                                                          <  2:164231/148‑1 (MQ=255)                                                                                             aTCCACCGCCATGTGCAAAATTCCCTT‑CGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa                                                      <  2:405279/149‑1 (MQ=255)                                                                                                caccGCCATGTG‑AAAAATTCCTTCCGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt‑gtg                                                <  2:698493/148‑1 (MQ=255)                                                                                                 aaCGCCATGTAAAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAT‑GTgg                                               <  2:892476/151‑1 (MQ=255)                                                                                                          tgaaaaaaTGCCTT‑CGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAT‑GTGGAAAGGGa                                        <  2:776130/146‑1 (MQ=255)                                                                                                             aaaaaTGCCTT‑CGCGGCGGCTGTCTCCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAT‑GTGGAAAGGGa                                        <  1:755450/146‑1 (MQ=255)                                                                                                               aaaTGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGAAAGA‑G‑AAAAAGA‑TTTTAACCCCCCCTGTTTTGTGGGGAACGCTGCATGTAACAATGGGGGAAAGGGagggggg                                  >  2:910811/1‑145 (MQ=255)                                                                                                                   gCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑G‑AAAAAAC‑TTTTAACCCCCCCTGTTGTGTGGGGATAGCTGCTGGGAACAAA‑GTGGAAAGGGGAgggggggg                               >  2:873545/1‑142 (MQ=255)                                                                                                                    cttt‑CGGGGCGGTTGTTTGCTAGATTTTTTCCCCCCGAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAT‑GTGGAAAGGGAAGTGCGGCGaa                             <  2:629261/148‑1 (MQ=255)                                                                                                                     cttccgcgGCGGCTGTCTGCCGGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAT‑GTGGAAAGGGAAGTGCGGCGa                              <  2:374370/145‑1 (MQ=255)                                                                                                                          gcggcgGCTTTCTGCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAT‑GTGGAAAGGGAAGTGCGGCGAATCa                          <  2:80058/148‑1 (MQ=255)                                                                                                                          gcggcgGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGCTCTCTCCAGTAGTGGCGGGGGCAAA‑G‑AAAAAACATTTAAACCC‑CCATGTTGTGTGGGGAACTCTGGTGGTAAAAAA‑GTGGGAAGGGGAgggggg                                 >  2:436789/1‑136 (MQ=255)                                                                                                                            ggcggcTTTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAT‑GTGGAAAGGGAAGTGCGGCGAATCAGCt                       <  2:296178/149‑1 (MQ=255)                                                                                                                            ggcggcTTCCTGCCGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAT‑GTGGAAAGGGAAGTGCGGCGAATCAGCTc                      <  1:629295/150‑1 (MQ=255)                                                                                                                                    tcggCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAT‑GTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc               <  2:823446/146‑1 (MQ=255)                                                                                                                                        ccAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAT‑GTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg            <  2:686178/148‑1 (MQ=255)                                                                                                                                        ccAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAT‑GTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc             <  2:1020616/147‑1 (MQ=255)                                                                                                                                        ccAGCTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAT‑GTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc             <  2:647921/147‑1 (MQ=255)                                                                                                                                                tttcCCGCCAAAAGAAAACCGCCGTTCTCGCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAT‑GTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt  <  1:884232/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                       TAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAT‑GTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357735‑358026

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