Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,556,207	0	T	C	48.0%	130.0 / 32.1	100	S215G (AGT→GGT)	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (45/7);  new base C (0/48);  total (45/55)
Fisher's exact test for biased strand distribution p-value = 2.34e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCAC  >  NC_000913/3556059‑3556351                                                                                                                                                      |                                                                                                                                                 gCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTcc                                                                                                                                                <  1:951223/151‑1 (MQ=255) gCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTcc                                                                                                                                                >  1:113463/1‑151 (MQ=255)  ccGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTccc                                                                                                                                               <  1:1345276/151‑1 (MQ=255)    gCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCg                                                                                                                                             <  1:950764/151‑1 (MQ=255)    gCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCg                                                                                                                                             >  1:488690/1‑151 (MQ=255)    gCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCg                                                                                                                                             >  1:1426510/1‑151 (MQ=255)      aaCCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGcc                                                                                                                                           >  2:187450/1‑151 (MQ=255)      aaCCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGcc                                                                                                                                           >  2:1152591/1‑151 (MQ=255)      aaCCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGcc                                                                                                                                           <  1:1271022/151‑1 (MQ=255)         caccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAgt                                                                                                                                        >  2:1405264/1‑151 (MQ=255)          accTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAgtg                                                                                                                                       >  1:1001527/1‑151 (MQ=255)           ccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAgtgt                                                                                                                                      >  1:1236238/1‑151 (MQ=255)            cTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAgtgt                                                                                                                                      >  2:847814/1‑150 (MQ=255)            cTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAgtgt                                                                                                                                      >  2:1450769/1‑150 (MQ=255)            cTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAgtgt                                                                                                                                      >  1:1444152/1‑150 (MQ=255)              cGGCACTGCCCCCCTTTCCACCGACGACAAAAAAGCGTTCGGTGATTTTTTCTCTTCCGACTCCAAGCGAGACGTTATCCCCCGCCCCCTGGCCGCGCGGCAGTTTAGAAATATTCTTTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGc                                                                                                                                   <  1:1006602/151‑1 (MQ=255)                    ggACGCGCTTTCCACCGACAAAAAAAAAGCGTGGGGTGATTTTTTCTTTTTCGACTTCAAGCGAGACGTTATCCCCCGCCCCCGGGCCGCGGGGCAGTTTATAAATATTCTGTTCAGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTc                                                                                                                             <  1:1118953/150‑1 (MQ=255)                     gACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTc                                                                                                                             >  2:1112260/1‑150 (MQ=255)                       cgcgCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACg                                                                                                                          >  2:466296/1‑151 (MQ=255)                         cgcTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCt                                                                                                                        >  2:803339/1‑151 (MQ=255)                          gcTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                                       >  2:1154894/1‑151 (MQ=255)                           cTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                                       >  1:164920/1‑150 (MQ=255)                             ttcccccGACAAAAAAAAGGCGTTGGGTGATTTTTTCCCTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCGGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCa                                                                                                                      <  2:374571/144‑1 (MQ=255)                              tcccccGACAAAAAAAAAGCTTTCGGTGATTTTTCCCCTTCCGCCTACAAGCGAGCCGGTATCCCCCGCCCCCTGGTCGCGCGGCAGTTTATGAATATTCTGTTCATGCAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGc                                                                                                                    <  1:691981/146‑1 (MQ=255)                                ccccgccaaaaaaaagGCGTTGGGTAATATTTTCACTTCCGCCTCCAAGCGAGACGTTATCCCCCGGCCCCGGGCCGCGCGGCGGTTTATAAATTTTCTGTTCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAAc                                                                                                                 <  1:899704/141‑1 (MQ=255)                                cACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAAc                                                                                                                 >  1:567036/1‑151 (MQ=255)                                 aCCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAAc                                                                                                                 >  1:1049009/1‑150 (MQ=255)                                  ccGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAt                                                                                                               >  2:79206/1‑151 (MQ=255)                                  ccGACGAAAAAAAAGCGTTCGGTGATTTTTTCCCTTCCGCCTTCAAGCGAGACGTTATTCCCCGGCCCCTGGCCGCGCGGCAGTTTATGAATATTCTGTTCATGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAt                                                                                                               <  2:997417/151‑1 (MQ=255)                                   cgacaaAAAAAAAGCGTTGGGTGATATTTTCACTTTCGACTCCAAGCGAGACGTTATTCCCCGGCCCCTGGTCGCGCGGCGGGTTATGAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg                                                                                                              <  1:262256/151‑1 (MQ=255)                                    gacgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATAAATATTCTGTTCATGCGGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgg                                                                                                             <  2:544616/151‑1 (MQ=255)                                     acgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc                                                                                                            >  2:573178/1‑151 (MQ=255)                                     acgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc                                                                                                            >  2:892492/1‑151 (MQ=255)                                       gacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                                          >  2:1251709/1‑151 (MQ=255)                                           aaaaaGCGTTCGGTGATTTTTTCACTTTCGACTTAAAGCGAGACGGTATTCCCCGCCCCCTGGTCGCGCGGCGGTTTATAAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                                      <  2:1463112/151‑1 (MQ=255)                                             aaaGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcac                                                                                                     >  1:462035/1‑150 (MQ=255)                                              aaGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCacac                                                                                                   >  2:1102534/1‑151 (MQ=255)                                                 cGTTCGGTGATTTTTTCACTTTCGACTCCAAGCGAGACGGTATTCCCCGCCCCCTGGTCGCGCGGAAGTTTATGAATATTCTTTTCATGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAg                                                                                                <  2:1338777/151‑1 (MQ=255)                                                  gTTGGGTAATTTTTTCACTTTCGACTTCAAGCGAGAGGGTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATAAATATTTTGTTCATGAAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGc                                                                                               <  2:1083748/151‑1 (MQ=255)                                                      ggTGATTTTTTCCCTTTCGCCTTCAAGGGAGACGTTTTCCCCCGGCCCCTGGTCGCGCGGCAGTTTATAAATTTTCTGTTCAGGCGGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAAt                                                                                           <  2:741299/151‑1 (MQ=255)                                                      ggTGATTTTTTCCCTTTCGCCTTCAAGCGAGCCGGTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATAAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAAt                                                                                           <  2:1235694/151‑1 (MQ=255)                                                       gTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAAt                                                                                           >  1:134223/1‑150 (MQ=255)                                                        tGATTTTTTCCCTTTCGACTAAAAGCGAGGCGTTTTCCCCCGCCCCCGGGTCGCGGGGCGGTTTATAAATTTTCTGTTCATG‑AGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATa                                                                                          <  1:320781/150‑1 (MQ=255)                                                         aatttttccccTTTCGACTACAAGCGAGACGTTATTCCCCGGCCCCTGGCCGCGCGGCAGGTTATAAATTTTCTGTTCAGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa                                                                                        <  1:733053/141‑1 (MQ=255)                                                          attttttCTCTTTCGACTTAAGGGGAGACGTTATTCCCCGCCCCGTGGTGGCGCGGCGGGTTATAAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGaa                                                                                       <  1:1461652/148‑1 (MQ=255)                                                            tttttCTCTTTCGCTTTAAAGCGAGAGGGTATTCCCCGCCCCCTGGCCGCGCGGCGGTTTATGAATTTTCTTTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                                      <  2:530315/149‑1 (MQ=255)                                                            tttttCACTTCCGCCTTAAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATAAATATTCTGTGCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                                      <  2:371846/149‑1 (MQ=255)                                                              tttCTCTTTCGCCTTCAAGCGAGACGGTATTCCCCGGCCCCTGGCCGCGCGGCGGTTTATGAATTTTCTGTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACtt                                                                                    <  1:50879/150‑1 (MQ=255)                                                                  aCTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTcc                                                                                >  1:946670/1‑150 (MQ=255)                                                                     ttCGACTTCAGGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATAAATTTTCTGTTCATGCAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                               <  2:761216/148‑1 (MQ=255)                                                                       cGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                                          >  2:727022/1‑151 (MQ=255)                                                                         aCTTCAAGCGAGACGGTATTCCCCGGCCCCTGGCCGCGCGGCAGGTTATGAATATTCTGTTCATGAAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                        <  2:1464784/151‑1 (MQ=255)                                                                            acAAGCGAGCCGGTATTCCCCGGCCCCTGGCCGCGGGGCGGTTTATAAATTTTCTGTTCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                                       <  1:1428876/148‑1 (MQ=255)                                                                              aaGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATAAATTTTCTGTTCATGCAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcaca                                                                   <  2:1013969/151‑1 (MQ=255)                                                                                  gagcggttttcccccGGCCCCGGGCGGCGCGGCAGTTTATAAATATTTTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                               <  1:554316/139‑1 (MQ=255)                                                                                   agaCGTTATTCCCCGTCCCCTGGGCGCGCGGCTTTTAATGATTATTCTGTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt                                                              <  1:276549/151‑1 (MQ=255)                                                                                   agaCGGTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATGAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt                                                              <  2:512177/151‑1 (MQ=255)                                                                                    gaCGTTATTCCCCGCCCCCTGGTCGCGGGGCAGGTTATAAATTTTTTTTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTc                                                             <  1:1242381/151‑1 (MQ=255)                                                                                      cGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccc                                                           >  1:623758/1‑151 (MQ=255)                                                                                       ggTATTCCCCGGCCCCGGGTCGCGCGGCAGTTTATGAATATTCTGTTCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                          <  2:156654/151‑1 (MQ=255)                                                                                       ggTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                          >  1:397467/1‑151 (MQ=255)                                                                                          aTTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCgcgc                                                       >  2:717180/1‑151 (MQ=255)                                                                                                cGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                                  >  1:75840/1‑150 (MQ=255)                                                                                                cGGCCCCTGGCCGCGCGGCAGTTTATGATTATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCa                                                 <  1:139178/151‑1 (MQ=255)                                                                                                  cccccTGGTCGCCCGGCAGTTTATAATTTTTCTTTCAATGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg                                               <  1:1530305/150‑1 (MQ=255)                                                                                                   ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGt                                              >  2:6143/1‑151 (MQ=255)                                                                                                    cccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGtt                                             >  1:14261/1‑151 (MQ=255)                                                                                                    cccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGtt                                             >  2:158173/1‑151 (MQ=255)                                                                                                       tgccgccgcgGAAGTTAATGAATTTTCTTTTCATGGAGAAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTg                                            <  1:558668/143‑1 (MQ=255)                                                                                                        ggTCGCGCGGCAGGTTATGAATATTCTTTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                                         <  2:650589/151‑1 (MQ=255)                                                                                                         gccgcgcgGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAg                                        <  2:1480162/149‑1 (MQ=255)                                                                                                              gcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt                                   >  1:751596/1‑151 (MQ=255)                                                                                                               gggCAGTTTATGAATATTCTGTTCATGAAGAAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTa                                  <  1:709844/150‑1 (MQ=255)                                                                                                               cgGCAGGTTATGATTATTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTa                                  <  1:385707/151‑1 (MQ=255)                                                                                                                ggCAGGTTATAATTTTTTTTTTCAGGAAGAAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaa                                 <  1:1071331/151‑1 (MQ=255)                                                                                                                 gCAGTTTATAAATTTTCTTTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaaa                                <  2:690679/151‑1 (MQ=255)                                                                                                                 gCAGGTTATGAATTTTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaaa                                <  2:366517/151‑1 (MQ=255)                                                                                                                 gCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaa                                 >  1:71563/1‑150 (MQ=255)                                                                                                                  cAGTTTATGAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                               <  2:786191/151‑1 (MQ=255)                                                                                                                  cAGGTTATAAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                               <  1:917407/151‑1 (MQ=255)                                                                                                                   ggtttATAAATTTTTTTTCCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACg                              <  1:912939/148‑1 (MQ=255)                                                                                                                   gggTTATAAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaaaag                              <  2:1132053/150‑3 (MQ=255)                                                                                                                   aGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACg                              >  2:1478973/1‑151 (MQ=255)                                                                                                                   aGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACg                              >  1:1033799/1‑151 (MQ=255)                                                                                                                   aGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACg                              >  2:1051695/1‑151 (MQ=255)                                                                                                                      ttATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATg                           >  2:1037460/1‑151 (MQ=255)                                                                                                                        ataaaTTTTCTGTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                          <  2:1425187/147‑1 (MQ=255)                                                                                                                             tatTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg                    >  2:1196989/1‑151 (MQ=255)                                                                                                                              atTTTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc                   <  2:1495341/151‑1 (MQ=255)                                                                                                                              atTCTTTTCAGGCAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg                    <  1:749297/150‑1 (MQ=255)                                                                                                                                 ttGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCga                 <  1:1299742/149‑1 (MQ=255)                                                                                                                                   gTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGAc              >  2:1116457/1‑151 (MQ=255)                                                                                                                                     tCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACt             <  1:1152419/150‑1 (MQ=255)                                                                                                                                      cATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTc           >  1:1494508/1‑151 (MQ=255)                                                                                                                                       aTGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc          <  2:1426766/151‑1 (MQ=255)                                                                                                                                        ggaaGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCg         <  1:487534/148‑1 (MQ=255)                                                                                                                                          cAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGtt       <  2:978483/151‑1 (MQ=255)                                                                                                                                           aGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTg      >  1:450325/1‑151 (MQ=255)                                                                                                                                             aaaaaac‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGcc    <  1:178553/149‑1 (MQ=255)                                                                                                                                               aaaaa‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGccac  >  2:1391338/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                 GCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCAC  >  NC_000913/3556059‑3556351

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