Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,556,200	0	G	A	43.7%	82.7 / 34.8	71	S217F (TCC→TTC)	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base G (35/5);  new base A (0/31);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 3.39e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGC  >  NC_000913/3556051‑3556348                                                                                                                                                      |                                                                                                                                                     ccAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAggg                                                                                                                                                      <  1:400729/150‑1 (MQ=255)  cAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGa                                                                                                                                                     >  1:721959/1‑150 (MQ=255)   aaCTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGaaa                                                                                                                                                   <  2:933971/151‑1 (MQ=255)       gtgCCGCAACCACCTCGGCACTGACGCGCTGTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACGTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGaaaaa                                                                                                                                                 <  2:97225/149‑1 (MQ=255)         gCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTcc                                                                                                                                             >  1:719232/1‑151 (MQ=255)          ccGCAACCACCTCGGCACTGACGCGCTTTTCATCGACGACAAACAAGCGTTCGGTGATATTTTCACTTTCGACTTCACGCGCGTCGGAATTAGCCGGCCCCTGGTCGCCCGGCAGGTTATGAATGTTCTGTTCATGCAGGGAGAAACTccc                                                                                                                                            >  2:106378/1‑151 (MQ=255)          ccGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTccc                                                                                                                                            >  2:985564/1‑151 (MQ=255)          ccGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTccc                                                                                                                                            <  2:368508/151‑1 (MQ=255)           cGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTcccc                                                                                                                                           >  1:445575/1‑151 (MQ=255)            gCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTcccc                                                                                                                                           >  2:748016/1‑150 (MQ=255)            gCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCg                                                                                                                                          <  1:266557/151‑1 (MQ=255)              aaCCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGCTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGcc                                                                                                                                        >  2:912493/1‑151 (MQ=255)                  accTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAgtg                                                                                                                                    >  1:382541/1‑151 (MQ=255)                         cACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATAAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGAt                                                                                                                             >  1:1094517/1‑151 (MQ=255)                             gACGCGTTTTTCACCGACGAAAAAAAAGCGTTGGGTAATTTTTCCATTTCCGTCTACAAGCGGGACGTTATTCCCCGCCCCCGGGGCGCGCGGCGGTTAATAAATATCCTGTTCATGCAGGAAAAAACTCCCCGCCAGTGTAGCGATTTCa                                                                                                                         <  2:525887/151‑1 (MQ=255)                               cgcgCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACg                                                                                                                       >  1:533214/1‑151 (MQ=255)                                  gcTTTTCACCGACGACAAAAAAGAGTTCGGTGATATTTTCACTTTCGCCTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                                    >  2:881285/1‑151 (MQ=255)                                    ttttCACCGACGAAAAAAAAGCGTTCGGTGATATTTTCTCTTTCGACTTCAAGCGAGACGTTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATAAATATTCTTTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                                    <  1:1042746/149‑1 (MQ=255)                                       tCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCa                                                                                                                >  2:9405/1‑150 (MQ=255)                                          ccGACGACAAAAAGGCGTTGGTTGTTTTTTCCACTTTCGCCTTAAAGCGAGACGGTATCCCCCGCCCCGTGGCCGGGGGGTGGGTTATAAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCa                                                                                                                   <  2:890867/144‑1 (MQ=255)                                          ccGACGAAAAAAAGGCGTTCGGTGATTTTTTCACTTTCGACTTCAAGCGAGACGTTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATGAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAAc                                                                                                              <  1:627691/149‑1 (MQ=255)                                              cgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcg                                                                                                        >  2:1308034/1‑151 (MQ=255)                                              cgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcg                                                                                                        >  2:1042160/1‑151 (MQ=255)                                                acAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                      >  1:359308/1‑151 (MQ=255)                                                aaaaaaagGCGTTGGGTGATATTTTCACTTGCGCCTACAAGCGAGGCGTTATTCCCCGCCCCCTGGTCGCGCGGGAGTTTATAAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                                       <  1:258711/148‑1 (MQ=255)                                                aaaaaaaaGCGTTGGGTGATTTTTTCACTTTCGACTTAAGGCGAGACGGTATTCCCCGCCCCCTGGTCGCGCGGCGGTTTATAAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                      <  1:774561/149‑1 (MQ=255)                                                 aaaaaaaGCGTTGGGTGATTTTTTCCCTTTCGACTTAAAGCGAGACGTTATTCCCCGCCCCCGGGCCGCGCGGCAGGTTATAAATATTCTGTTAATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACa                                                                                                             <  1:681376/142‑1 (MQ=255)                                                  aaaaaaGCTTTGGGTGATTTTTTCCCTTTCGACTCCAAGCGAGAGGTTTTTCCCCGGCCCTTGGGCGCGCGGCGGTTTATGAATTTTCTGTTCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                                    <  1:529180/151‑1 (MQ=255)                                                  aaaaaaGCGTTCGTTGATTTTTCCCCTTGCAACTACAAGCGAGACGGTATTACCCGCCCCCTGGTCGCGGGGAAGGTTATAAATATTCTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc                                                                                                         <  1:298054/146‑1 (MQ=255)                                                       aGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACAcc                                                                                               >  2:1098753/1‑151 (MQ=255)                                                         cGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACAcc                                                                                               >  1:1254758/1‑149 (MQ=255)                                                           ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                           >  2:131458/1‑151 (MQ=255)                                                              ggTAATTTTTTCATTTTCGACTTCAGGGGAGACGGTATTCCCCGGCCCCTGGCCGCGCGGCAGTTTATAATTATTTTGTTCATGGAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa                                                                                              <  2:1290834/145‑1 (MQ=255)                                                               gTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAAt                                                                                        >  2:1108251/1‑150 (MQ=255)                                                                tGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAAt                                                                                        >  2:1032423/1‑149 (MQ=255)                                                                   ttttttCACTTTGGCCTTAAAGGGGGGCGGTATTCCCCGGCCCGGGGTCGCGCGGAAGGTTATGAATATTCTGTTCAGGAAGAAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGaa                                                                                    <  2:1204486/148‑1 (MQ=255)                                                                     ttttCGTTTTGGACTTCAGGGGAGAGGTTGTCACCCGCCCGCAGTTCGCGGGGGTGGTTGTGAATATTCTGTTCAAGGAGAAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAg                                                                                             <  1:890614/139‑1 (MQ=255)                                                                     ttttCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACtt                                                                                 >  1:960188/1‑151 (MQ=255)                                                                      tttCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtc                                                                                >  1:34442/1‑151 (MQ=255)                                                                         cACTTTCGACTCAAAGCGAGCCGGTATTCCCCGCCCCCGGGTCGCGCGGCAGGTTATAAATATTCTGTTCATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtct                                                                               <  2:478148/149‑1 (MQ=255)                                                                             tccGACTTCAAGCGAGACGGTATTCCCGGGCCCCTGGTCGCGCGGCGGTTTATAATTTTTCTGTTCATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                          <  1:782095/148‑1 (MQ=255)                                                                                gACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATc                                                                      >  1:1089383/1‑151 (MQ=255)                                                                                   ttCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGc                                                                   >  1:635853/1‑151 (MQ=255)                                                                                   ttCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGc                                                                   >  2:290120/1‑151 (MQ=255)                                                                                      acgCGAGAGGGTTTTCCCCGCCCCCGGGCCGCGCGGCGGGTTATGAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcaca                                                                <  2:520767/149‑1 (MQ=255)                                                                                      aaGCGAGACGTTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                                    <  1:1327201/147‑1 (MQ=255)                                                                                      aaGCGAGACGGTATACCCCGGCCCCTGGCCGCGCGGCAGGTTATAAATATTTTGTTCAGGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGc                                                                   <  1:512667/148‑1 (MQ=255)                                                                                        gCGAGACGTTATTCCCCGCCCCCTGGTCGCGCGGGGGTTTATGAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACaa                                                               <  1:1188300/150‑1 (MQ=255)                                                                                          gagaCGTTATTCCCCGGCCCCGGGCCGCGCGGCAGTTTATAAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACaa                                                               <  1:891563/148‑1 (MQ=255)                                                                                            gaggttATTCCCGGGCCCCTGGTCGCGGGGCAGTTTATGAATTTTTTGTTCATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGc                                                                   <  1:826591/137‑1 (MQ=255)                                                                                            gaggttATTCCCCGGCCCCTGGGCGCGGGGCAGGTTATAAATTTTCTGTTCATGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACaa                                                               <  1:330429/141‑1 (MQ=255)                                                                                                gTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                      >  2:1228173/1‑151 (MQ=255)                                                                                                 tATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgcg                                                     >  2:1111300/1‑151 (MQ=255)                                                                                                  aTTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCgcgc                                                    >  1:573804/1‑151 (MQ=255)                                                                                                  aTTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCgcgc                                                    >  2:642934/1‑151 (MQ=255)                                                                                                    ccccgcgCCCTTGGCCGGGCGGCATTTTATGATTTGTCTGTTCGCGGAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCt                                                   <  2:909226/144‑1 (MQ=255)                                                                                                      cccGCCCCCGGGTCGGGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCg                                                 <  1:949652/150‑1 (MQ=255)                                                                                                        cGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                               >  2:547645/1‑150 (MQ=255)                                                                                                           ccccGGGTCGCGCGGCAGGTTATAAATATTCTGTTCATGCGGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGt                                           <  1:1102005/151‑1 (MQ=255)                                                                                                                        ggCAGTTTTTAAATATTCTGTTCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTa                               <  1:830841/150‑1 (MQ=255)                                                                                                                        ggCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTa                               >  2:486842/1‑150 (MQ=255)                                                                                                                          cAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                            >  2:1266754/1‑151 (MQ=255)                                                                                                                            ggTTATGAATTTTCTTTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaaa                             <  1:1224611/148‑1 (MQ=255)                                                                                                                              ttATGAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                            <  1:369526/147‑1 (MQ=255)                                                                                                                              ttATGAATTTTCTGTTCATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACg                           <  1:1108274/148‑1 (MQ=255)                                                                                                                               tATGAATTTTCTTTTCAGACGGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATg                        <  2:874885/150‑1 (MQ=255)                                                                                                                                 tGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc                      >  1:1081307/1‑150 (MQ=255)                                                                                                                                   aaTTTTCTTTTCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCAc                    <  2:975459/150‑1 (MQ=255)                                                                                                                                     tatTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg                 >  2:325406/1‑151 (MQ=255)                                                                                                                                        ttttttCATGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCg               <  2:880308/146‑1 (MQ=255)                                                                                                                                              cATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTc        >  2:90003/1‑151 (MQ=255)                                                                                                                                              cATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTc        <  1:1097273/151‑1 (MQ=255)                                                                                                                                                    gggAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGc  >  1:115194/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGC  >  NC_000913/3556051‑3556348

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