Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,243,742	T→A	5.9%	intergenic (+31/+68)	tsaA → / ← rnt	peroxiredoxin/ribonuclease T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,243,742	0	T	A	5.9%	205.4 / 2.7	68	intergenic (+31/+68)	tsaA/rnt	peroxiredoxin/ribonuclease T
Reads supporting (aligned to +/- strand):  ref base T (30/34);  new base A (2/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01

CCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCGGT  >  NC_002947/1243597‑1243890                                                                                                                                                  |                                                                                                                                                      ccACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGc                                                                                                                                                    >  4:160873/1‑149 (MQ=255)  cACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGcc                                                                                                                                                   >  7:46896/1‑149 (MQ=255)      aaGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGtt                                                                                                                                               <  6:156388/149‑1 (MQ=255)       agagCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGttt                                                                                                                                              >  7:116490/1‑149 (MQ=255)         agCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGttttt                                                                                                                                            <  8:116490/149‑1 (MQ=255)             cGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt                                                                                                                                        <  1:23118/149‑1 (MQ=255)              ggCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATg                                                                                                                                       >  6:66890/1‑149 (MQ=255)                cGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGcc                                                                                                                                     >  3:19242/1‑149 (MQ=255)                   aGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTc                                                                                                                                  <  7:109461/149‑1 (MQ=255)                       tgCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCt                                                                                                                                    <  2:142957/143‑1 (MQ=255)                       tgCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCt                                                                                                                                    >  1:142957/1‑143 (MQ=255)                         cccTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGc                                                                                                                            <  2:96040/149‑1 (MQ=255)                         cccTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGc                                                                                                                            >  4:277412/1‑149 (MQ=255)                            tGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCgg                                                                                                                         >  5:128410/1‑149 (MQ=255)                              ccAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGc                                                                                                                        <  8:46896/148‑1 (MQ=255)                                aaCTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCctc                                                                                                                     >  7:84780/1‑149 (MQ=255)                                     gAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTcgc                                                                                                                >  2:84836/1‑149 (MQ=255)                                     gAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTcgc                                                                                                                >  6:51555/1‑149 (MQ=255)                                     gAAAAAAGGCGACAAGGGCAGGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTcg                                                                                                                 <  4:19242/148‑1 (MQ=255)                                          aaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGgcg                                                                                                           >  6:112197/1‑149 (MQ=255)                                             gCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGgcgcgc                                                                                                        <  6:279647/149‑1 (MQ=255)                                               gACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTGTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGccc                                                                                                      >  3:35554/1‑149 (MQ=255)                                                   aGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGa                                                                                                 <  5:51555/149‑1 (MQ=255)                                                        aTGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGcc                                                                                            >  4:196759/1‑149 (MQ=255)                                                              gCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTCCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAg                                                                                       >  1:336341/1‑148 (MQ=255)                                                                  cGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAAGCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGa                                                                                   <  6:207777/148‑1 (MQ=255)                                                                  cGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTCCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGa                                                                                   <  2:336341/148‑1 (MQ=255)                                                                   tccTGAAGGCGTTGCTGCCTCCCTGCGCGCGTACGCGGGCATGCTGTACTTGCTGCGTGCATGAAACAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATc                                                                                 <  4:43573/148‑1 (MQ=255)                                                                   gCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAAt                                                                                               >  7:200767/1‑135 (MQ=255)                                                                   gCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAAt                                                                                               <  8:200767/135‑1 (MQ=255)                                                                    ccTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATc                                                                                 >  8:171849/1‑148 (MQ=255)                                                                     cTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAg                                                                                       >  7:239476/1‑141 (MQ=255)                                                                     cTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGAt                                                                                  <  8:239476/146‑1 (MQ=255)                                                                     cTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCaa                                                                               >  7:300530/1‑149 (MQ=255)                                                                     cTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCa                                                                                >  4:315478/1‑148 (MQ=255)                                                                        aaGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATc                                                                             >  6:209882/1‑148 (MQ=255)                                                                        aaGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATc                                                                             <  8:96958/148‑1 (MQ=255)                                                                        aaGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATc                                                                             <  5:209882/148‑1 (MQ=255)                                                                        aaGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATc                                                                             <  3:315478/148‑1 (MQ=255)                                                                          ggCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGtt                                                                          <  8:84780/149‑1 (MQ=255)                                                                          ggCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGtt                                                                          <  3:196759/149‑1 (MQ=255)                                                                          ggCGCTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGtt                                                                          <  3:155850/149‑1 (MQ=255)                                                                           gCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTg                                                                         <  5:112197/149‑1 (MQ=255)                                                                              tttctgcCTACCTGATCTAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATTAAAAAACCGGCCTTCGTGGCCGTTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGaaa                                                                      <  3:185189/146‑1 (MQ=255)                                                                                  tgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATcgc                                                                  >  4:70597/1‑149 (MQ=255)                                                                                   gcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGAt                                                                                  <  4:79112/132‑1 (MQ=255)                                                                                    ccTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCgcgc                                                                <  3:160873/149‑1 (MQ=255)                                                                                         cTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgcc                                                           <  5:66890/149‑1 (MQ=255)                                                                                                   aCGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttct                                                 >  5:19616/1‑149 (MQ=255)                                                                                                   aCGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttct                                                 >  6:26044/1‑149 (MQ=255)                                                                                                     gCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctt                                                <  1:84836/149‑1 (MQ=255)                                                                                                         gCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCa                                            <  5:26044/148‑1 (MQ=255)                                                                                                          cAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttcttt                                               <  4:4696/144‑1 (MQ=255)                                                                                                          cAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttcttt                                               >  3:4696/1‑144 (MQ=255)                                                                                                               tGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTATTTCCacatcat                                      >  1:287619/1‑142 (MQ=255)                                                                                                                 ttaTTGCTGCGTGCATGAAAAAACCGGGGTTCGTGGCCGGTTTTTTTATGCATTCCTGTGCGGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCATTCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGa                                    <  8:30016/146‑1 (MQ=255)                                                                                                                 tAATTGCTGCTTGCATGAAAAAACCGGCCTTCGTTGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGttta                                    <  6:128410/148‑3 (MQ=255)                                                                                                                  aaTTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTAAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGAc                                   >  6:293659/1‑148 (MQ=255)                                                                                                                   aTTGCTGCGTTCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACg                                  <  4:202456/148‑1 (MQ=255)                                                                                                                   aTTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCTGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGa                                 <  3:277412/149‑1 (MQ=255)                                                                                                                               ataaaaaaaCCGGCCACGAAGGCCGGTTTTTTCATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAt                                                     >  1:42987/4‑117 (MQ=37)                                                                                                                               ataaaaaaaCCGGCCACGAAGGCCGGTTTTTTCATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAt                                                     <  2:42987/114‑1 (MQ=37)                                                                                                                                 aaaaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCa                                                              <  3:88401/105‑1 (MQ=38)                                                                                                                                 aaaaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCa                                                              >  4:88401/2‑106 (MQ=38)                                                                                                                                  aaaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAAGCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGCACAGCt                  >  5:137770/1‑149 (MQ=255)                                                                                                                                    aaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGCTGCCGCAGAACAGCt                  >  7:146786/1‑147 (MQ=255)                                                                                                                                          ggCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGtt          >  7:30365/1‑149 (MQ=255)                                                                                                                                                 gTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCg    <  3:70597/148‑1 (MQ=255)                                                                                                                                                  tGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCGgt  <  7:171849/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                      CCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCG‑CGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCGGT  >  NC_002947/1243597‑1243890

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