Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,556,207	0	T	C	42.1%	79.7 / 20.1	57	S215G (AGT→GGT)	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (28/5);  new base C (0/24);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.85e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACC  >  NC_000913/3556057‑3556352                                                                                                                                                        |                                                                                                                                                  gtgCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑Ct                                                                                                                                                   >  1:903191/1‑151 (MQ=255)       cAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGc                                                                                                                                             >  2:744031/1‑151 (MQ=255)        aaCCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAATAA‑CTCCCCGcc                                                                                                                                            >  1:541696/1‑151 (MQ=255)          ccaccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGcc                                                                                                                                            >  2:653083/1‑149 (MQ=255)               tCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAgtgt                                                                                                                                       >  1:270923/1‑149 (MQ=255)                  gCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCg                                                                                                                                   >  2:23857/1‑150 (MQ=255)                   cACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGAt                                                                                                                                 >  2:325413/1‑151 (MQ=255)                    aCTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGAtt                                                                                                                                >  1:963108/1‑151 (MQ=255)                      tgccgcgTTTTCCACCGACGAAAAAAAAGCGTTCGGTGATTTTTTCTCTTTCGACTTCAAGCGAGACGTTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAAA‑CTCCCCGCCAGTGTAGCGAtt                                                                                                                                <  2:635829/146‑1 (MQ=255)                      tgccgcgCTTTTCCCCGACGAAAAAAAAGCGTTGGGTGATTTTTTCCCTTTCGACTTCAAGCGAGACGGTATTCCCCGCCCCCTGGTCGCGCGGCAGTTTATAAATATTCTGTTCATG‑AGGGAAAAACCCCCCCGCCAGTGTAGCGAtt                                                                                                                                <  1:871699/146‑1 (MQ=255)                                    ccGACAAAAAAAAAGCGTTGGGTGATTTTTTCCCTTTCGACTCAAAGCGAGACGGTATTCCCCGCCCCCGGGTCGCGCGGCAGTTTATAAATATTCTGTTCATGAAGGGAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCaa                                                                                                                   <  1:455922/148‑1 (MQ=255)                                     cgaaaacAAAAAGGCGTTTGGTATTATTTCCACTTTGGACTAAAAGAGGCACGGTATTGCCCGTCCCTTGGTCGCGGGGACGGTTATTATTTTTCTCTACACGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg                                                                                                               <  1:141459/146‑1 (MQ=255)                                       acgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc                                                                                                             >  1:980027/1‑151 (MQ=255)                                        cgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc                                                                                                             >  2:390471/1‑150 (MQ=255)                                        cgacAAAAAACGGTTGGTTAATTTTTTCCTTTCCGACTTAAAGGGAGCCGTTATTCCCCGCCCCCTGGTCGCGCGGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACa                                                                                                                 <  1:916575/146‑1 (MQ=255)                                         gacAAAAATGCGTTCGCTGATTTTTCCCCTTTCGCCTCCAATAATACCGTTATGCCCCGCCCCCGGGTGGCCGGGCAGGTTATAAATTTTCTTTTCAAGGAGGAAAAAC‑CCCTCCGCCAGTGTAGCGATTTCACGCTCAGCAACa                                                                                                                 <  1:29154/145‑1 (MQ=255)                                          acAAAAAAGCGTTGGGTGATTTTTTCCCTTTCGCCTTCAAGCGAGACGGTATTCCCCGCCCCCTGGTCGCGGGGGAGTTTATAAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                                           <  2:796550/150‑1 (MQ=255)                                            aaaaaggcttgggttGATATTTTCACTTTCGCCTACAGGCGAGAGGGTATTCCCCGGCCCGTGGTCGCGCGGCGGTTTATGAATATTTTGTTCATGCGGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                          <  1:816972/135‑1 (MQ=255)                                            aaaaaaGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                                        >  2:876401/1‑151 (MQ=255)                                             aaaagGCTTTCGGTATTATTTCCCCTTCCGACTAAAAGGGAGCGGTTATCCCCCGCCCCCTGGCCGCGCGGGAGTTTATGAATTTTCTTTTCAGGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                                       <  2:9996/151‑1 (MQ=255)                                             aaaaaGCTTTGGGTGATTTTTCCCCTTCCGACTCCAAGCGAGCCGTTTTTCCCCGCCCCCTGGCCGCGCGGAAGTTTATAAATTTTCTTTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                                       <  1:316050/151‑1 (MQ=255)                                                     ttCGGTGATTTTTCCCCTTCCGACTTCAGGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATAAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGc                                                                                                <  1:657380/150‑1 (MQ=255)                                                     ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                               >  1:848517/1‑151 (MQ=255)                                                       gggTAATATTTCCCCTTTCGACTCCAAGGGAGAGGGTATTCCCCGGCCCCGGGCCGGGCGGCAGGTTATGAATATTCTTTTCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa                                                                                             <  1:256908/150‑1 (MQ=255)                                                        ggTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAAt                                                                                            >  2:681596/1‑151 (MQ=255)                                                          tGATTTTTTCACTTCCGACTCAAAGCGAGACGTTTTTCCCCGGCCCCGGGGCGGGGGGCGGGTTATAAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATa                                                                                           <  2:411343/150‑1 (MQ=255)                                                             tatTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                                       >  1:187152/1‑151 (MQ=255)                                                              atTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACt                                                                                      >  1:30078/1‑151 (MQ=255)                                                               ttttCTTTTTCGACTTCAAGCGAGACGGTATTCCCCGCCCCCTGGTCGCGCGGCGGTTTATAAATTTTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                                       <  1:480845/149‑1 (MQ=255)                                                               tttccccTTTCTCCTTCAAGCGAGCCGTTATCCCCCGCCCCTGGGCCGCGCGGCAGTTTATAAATTTTCT‑TTTATGAAGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACt                                                                                      <  2:486424/144‑1 (MQ=255)                                                               tttcCACTTCCGCCTAAAAGCGAGCGGTTATTCCCCGCCCCCGGGCCGCGCGGAAGTTTATAAATATTTTGTGCATGAAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa                                                                                         <  2:741849/147‑1 (MQ=255)                                                                tttCTCTTTCGTCTACAAGCGAGGCGTTATTACCGGCCCACTGGGCGCGGGGCAGGTTAAAAATATTTTGTTCATGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACtt                                                                                     <  2:67071/150‑1 (MQ=255)                                                                tttCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtc                                                                                    >  2:10383/1‑151 (MQ=255)                                                                          gACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATc                                                                          >  1:636036/1‑151 (MQ=255)                                                                             ttCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGc                                                                       >  1:309059/1‑151 (MQ=255)                                                                             tccAAGCGAGCCGTTTTTCCCCGGCCCCTGGTCGCGCGGCAGTTTAAGAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa                                                                            <  1:460191/144‑1 (MQ=255)                                                                                aaGCGATACGTTATCCCCGGCCCCCTGCCCGCGGGGCAGTTTATAAATTTTTTGTCCAAGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                              <  2:349105/141‑1 (MQ=255)                                                                                  gCGAGACGGTATCCCCCGGCCCCTGGTCGCGCGGCAGTTTATGAATATTCTGTTCAGGCAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcaca                                                                    <  1:312659/149‑1 (MQ=255)                                                                                     agccgttATTCCCCGGCCCCTGGCCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                                <  1:413758/144‑1 (MQ=255)                                                                                         ggTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                           >  1:766276/1‑151 (MQ=255)                                                                                           tATTCCCCGCCCCGGGGTCGCGCGGCAGTTTATAAATATTCTGTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTcc                                                             <  1:532509/147‑1 (MQ=255)                                                                                           tATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgcg                                                         >  1:753025/1‑151 (MQ=255)                                                                                            aTTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCgcgc                                                        >  2:747817/1‑151 (MQ=255)                                                                                              tACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                                      >  2:90574/1‑151 (MQ=255)                                                                                                 ccGCCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                                   <  2:903168/151‑1 (MQ=255)                                                                                                  cGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                                   >  2:507984/1‑150 (MQ=255)                                                                                                     ccccTGGTCGCGCGGAAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGt                                               >  1:615021/1‑151 (MQ=255)                                                                                                      cccGGGTCGCGCGGCGGGTTAAAAATATTCTGTCCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGtt                                              <  2:463660/151‑1 (MQ=255)                                                                                                                 cgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt                                    >  1:269248/1‑150 (MQ=255)                                                                                                                 cgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt                                    >  2:873700/1‑150 (MQ=255)                                                                                                                    cAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                                >  1:706641/1‑151 (MQ=255)                                                                                                                        ttATGAATTTTCTGTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACg                               <  1:217929/148‑1 (MQ=255)                                                                                                                          aTGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCCCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc                          >  2:628357/1‑151 (MQ=255)                                                                                                                             aaTATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg                       >  1:103935/1‑151 (MQ=255)                                                                                                                                  tCTTTTCATCCGGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCg                   <  2:850558/150‑1 (MQ=255)                                                                                                                                              ggaaaaac‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGt         <  2:194602/145‑1 (MQ=255)                                                                                                                                              ggaaaaac‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTg       <  2:320469/147‑1 (MQ=255)                                                                                                                                                  aaac‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGccacc  <  2:605846/151‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                  GTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACC  >  NC_000913/3556057‑3556352

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