Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,361,072	A→T	8.5%	intergenic (‑79/+22)	PP_t49 ← / ← PP_t48	tRNA‑Met/tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,361,072	0	A	T	8.5%	271.0 / 12.2	94	intergenic (‑79/+22)	PP_t49/PP_t48	tRNA‑Met/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base A (49/37);  new base T (4/4);  total (53/41)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

GCTGGATTTGAACCAACGACCTTCGGGTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGTTACTTCTGCTCTTGA  >  NC_002947/5360928‑5361216                                                                                                                                                  |                                                                                                                                                   gCTGGATTTGAACCAACGACCTTCGGGTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTTTAGCg                                                                                                                                                >  4:100694/1‑149 (MQ=255)  cTGGATTTGAACCAACGACCTTCGGGTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCgg                                                                                                                                               >  5:17590/1‑149 (MQ=255)   tGGATTTGAACCAACGACCTTCGGGTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCggg                                                                                                                                              >  8:309989/1‑149 (MQ=255)   tGGATTTGAACCAACGACCTTCGGGTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCggg                                                                                                                                              >  6:300356/1‑149 (MQ=255)     gATTTGAACCAACGAACTTCGGGTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCggg                                                                                                                                              <  6:231595/147‑1 (MQ=255)      aTTTGACCCAACGACCTTCGGGTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGctc                                                                                                                                           <  8:256060/149‑1 (MQ=255)       tttGAACCAACGACCTTCGGGTTATGAGCCCGACG‑AGCTACTAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGctc                                                                                                                                           >  2:70382/1‑148 (MQ=255)         tGAACCAACGACCTTCGGGTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTTTAGCGGGCtttc                                                                                                                                         <  5:164024/148‑1 (MQ=255)          gAACCAACGACCTTCGGGTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTTTAGCGGGCtttct                                                                                                                                        <  3:100694/148‑2 (MQ=255)           aaCCAACGACCTTCGGGTTATGAGCCCGACGAAGCTACAAAGCTGCACCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCATGTTGAGGGTCAAATCCAAACATTCCCCTACAAGAAAGCCCGCTAAAGCGGGCtgtc                                                                                                                                         >  1:231602/1‑147 (MQ=39)            aCCAACGACCTTCGGGTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTTTAGCGGGCtttcttgt                                                                                                                                     <  2:297799/149‑5 (MQ=255)            aCCAACGACCTTCGGGTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTTTAGCGGGCtttcttgt                                                                                                                                     <  7:134919/149‑5 (MQ=255)              cAACGACCTTCGGGTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCtt                                                                                                                                   <  3:128500/149‑1 (MQ=255)              cAACGACCTTCGGGTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCtt                                                                                                                                   >  1:338620/1‑149 (MQ=255)               aaCGACCTTCGGGTTATGAGCCCGGCG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCtt                                                                                                                                   >  8:18778/1‑148 (MQ=255)               aaCGACCTTCGGGTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCAGCTAAAGCGGGCTCTGAAGCtt                                                                                                                                   >  8:67369/1‑148 (MQ=255)                     cTTCCGGTTATCAGCCCGACG‑ACCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAAACCAGACCCTGTTGAGGCTCAAATCCAAACAATCACCAACAAAAAAGCCCGCAAAAGCGAGCTCACAAGCTTTAAtatg                                                                                                                            >  8:232624/1‑146 (MQ=255)                      ttCGGGTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTgg                                                                                                                            >  6:120368/1‑148 (MQ=255)                        cGGGTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTAc                                                                                                                         >  6:97294/1‑149 (MQ=255)                         gggTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTAcc                                                                                                                        <  2:35074/149‑1 (MQ=255)                         gggTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTAc                                                                                                                         <  2:338620/148‑1 (MQ=255)                          ggTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTAcc                                                                                                                        >  6:295154/1‑148 (MQ=255)                           gTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTAcc                                                                                                                        >  7:25304/1‑147 (MQ=255)                           gTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCga                                                                                                                      <  6:17590/149‑1 (MQ=255)                                gAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGAGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCAGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAgt                                                                                                                 >  5:184527/1‑148 (MQ=255)                                 aGCCCGACG‑AGCTACCAAGCTTCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAgg                                                                                                                 <  7:309989/148‑1 (MQ=255)                                 aGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAgg                                                                                                                 <  3:71346/148‑1 (MQ=255)                                  gCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAgggg                                                                                                               <  5:220673/149‑1 (MQ=255)                                                ccAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGTGACTCGAAC‑‑CCCTa                                                                                                >  6:105036/1‑149 (MQ=255)                                                ccAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTa                                                                                                >  7:285075/1‑149 (MQ=255)                                                 cAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTa                                                                                                >  6:103510/1‑148 (MQ=255)                                                 cAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTa                                                                                                >  3:153979/1‑148 (MQ=255)                                                  aaGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTaca                                                                                              >  4:159605/1‑149 (MQ=255)                                                  aaGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTaca                                                                                              >  6:300913/1‑149 (MQ=255)                                                        ctCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTAt                                                                                        >  3:32127/1‑149 (MQ=255)                                                         tCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGaa                                                                                                        >  8:124921/1‑134 (MQ=255)                                                         tCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGaa                                                                                                        <  7:124921/134‑1 (MQ=255)                                                           cACCCCGCGACAAAGCTGGGGCAAAAGTATTAGACCGAACCCTTTTTAGTGTCAAATCTAAACATTCACCTACAACAATGCGCTATAAAGCGGGATGTCAAGCTTTATTTTGTACCGATAATGGGACTCGAAC‑‑CCCTACACCCTATgg                                                                                      <  6:184527/148‑1 (MQ=255)                                                            aCCACGCGACAAAGCTGTGGCAAAAGTATAATACAGAACTCAGTTCATGGTCAAATACAAACATTCACCTACAAGAAAGCCCGATAAACCGGGCTCTCAAGCTTTAATAGGTACAGAGAAGGGGACTCGACC‑‑CCCTACACCCTATcgg                                                                                     >  6:296089/1‑148 (MQ=255)                                                              cccGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGca                                                                                   >  5:17343/1‑148 (MQ=255)                                                                          cTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAg                                                                      <  1:313529/149‑1 (MQ=255)                                                                            ggggCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAgga                                                                    >  2:231124/1‑148 (MQ=255)                                                                             gggCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAgg                                                                     <  8:285075/147‑1 (MQ=255)                                                                             gggCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGt                                                                    <  5:103510/148‑1 (MQ=255)                                                                               gCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGc                                                                 <  6:17343/149‑1 (MQ=255)                                                                                cAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCg                                                                >  4:64329/1‑149 (MQ=255)                                                                                 cactgTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCgt                                                               <  1:325220/145‑1 (MQ=255)                                                                                 aaaaGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCgt                                                               <  5:300913/149‑1 (MQ=255)                                                                                   aaGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCgtgt                                                             >  4:145016/1‑149 (MQ=255)                                                                                   aaGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCgtgt                                                             <  1:70382/149‑1 (MQ=255)                                                                                        taAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTAc                                                         <  5:300356/148‑1 (MQ=255)                                                                                           gACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAAt                                                     >  1:243315/1‑149 (MQ=255)                                                                                            aCCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCAAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAAtt                                                    <  6:294425/149‑1 (MQ=255)                                                                                             ccGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTc                                                   <  4:32127/149‑1 (MQ=255)                                                                                                aaCCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTccac                                                <  5:120368/149‑1 (MQ=255)                                                                                                   ccTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAgg                                                                                                                 <  4:288768/83‑1 (MQ=255)                                                                                                   ccTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAgg                                                                                                                 >  3:288768/1‑83 (MQ=255)                                                                                                    cTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCaccacc                                            >  3:125245/1‑149 (MQ=255)                                                                                                    cTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCaccacc                                            >  5:204161/1‑149 (MQ=255)                                                                                                           gggTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCaa                                     >  2:154502/1‑149 (MQ=255)                                                                                                              tCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCTGCAAta                                   >  6:102396/1‑148 (MQ=255)                                                                                                              tCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAAtac                                  >  5:10730/1‑149 (MQ=255)                                                                                                                   tCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCACAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAAtactact                              >  3:342845/1‑148 (MQ=255)                                                                                                                    ccAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTAACGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGAGTGTCTACCAATTACACCACCTCGGCAAtactacta                             >  2:285828/1‑148 (MQ=255)                                                                                                                         cATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGTGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGaa                       <  5:97294/149‑1 (MQ=255)                                                                                                                              aCCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGCCTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCg                  >  2:350572/1‑149 (MQ=255)                                                                                                                               ccTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGt                 <  7:67369/149‑1 (MQ=255)                                                                                                                               ccTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGt                 <  4:342845/149‑1 (MQ=255)                                                                                                                                cTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGtt                <  3:145016/149‑1 (MQ=255)                                                                                                                                cTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGtt                <  3:159605/149‑1 (MQ=255)                                                                                                                                 tACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGTTa               >  5:339552/1‑149 (MQ=255)                                                                                                                                  aCAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAAc                     >  1:204027/1‑142 (MQ=255)                                                                                                                                  aCAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAAc                     <  2:204027/142‑1 (MQ=255)                                                                                                                                  aCAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGTTAc              >  2:295497/1‑149 (MQ=255)                                                                                                                                   cAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAgg                                                                                                                 >  4:113400/1‑51 (MQ=255)                                                                                                                                    ttgagaGCCCGCTTTAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAACAGGCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGTTAc              >  1:109151/6‑149 (MQ=255)                                                                                                                                    ttgagaGCCCGCTTTAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGTTACtt            >  1:5278/6‑149 (MQ=255)                                                                                                                                    aaGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAAc                     >  2:164775/1‑140 (MQ=255)                                                                                                                                    aaGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAAc                     <  1:164775/140‑1 (MQ=255)                                                                                                                                    aaGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGTTACtt            >  8:47014/1‑149 (MQ=255)                                                                                                                                    aaGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGTTACtt            <  8:25304/149‑1 (MQ=255)                                                                                                                                    aaGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGTTACt             >  2:363126/1‑148 (MQ=255)                                                                                                                                      gAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCGCAACCACCTCAAGGTAGCGTGTCTACCAATTCGACCACCTCGGCAATACTACTACTTTCAACCCGt                 >  7:178110/1‑142 (MQ=255)                                                                                                                                      gAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGAGACTCGCAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGTTACTTCt          >  1:87634/1‑149 (MQ=255)                                                                                                                                      gAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGCGGACTCGAAC‑‑CCCTACCCCCTAACGGCACAACCACCTCAAGCTTGCGTGTCTACCAAAACCACCACCTCGGCAATACTACTACAAGAAACCCCTAACATCt          >  3:306988/1‑149 (MQ=255)                                                                                                                                       agaGCCCGCTTTAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAgg                                                                                                                 >  2:128489/3‑47 (MQ=25)                                                                                                                                       aaaGCACGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGCGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGTTACTTCTg         <  6:313202/149‑1 (MQ=255)                                                                                                                                        aaGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGTTACTTCTGc        <  3:64329/149‑1 (MQ=255)                                                                                                                                        aaGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGTTACTTCTGc        <  7:18778/149‑1 (MQ=255)                                                                                                                                          gCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGTTACTTCTGctc      >  1:33230/1‑149 (MQ=255)                                                                                                                                           cccGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGTTACTTCTGctct     >  8:179265/1‑149 (MQ=255)                                                                                                                                            ccGCTAAAGCGGGCTCTCAAGCTTTAATTGGAACCTAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACAACCTCAAGGTTGCGTGTCTACCAATTCCACAACCTCGGCAATCCTACTACTTGAAACCCGTTACTTCTGCTCtt    <  1:231124/149‑1 (MQ=255)                                                                                                                                            ccGCTAAAGCGGGCTCTCAAGATTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGTTACTTCTGCTCtt    <  1:154502/149‑1 (MQ=255)                                                                                                                                              gCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGTTGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGTTACTTCTGCTCTTGa  <  5:105036/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GCTGGATTTGAACCAACGACCTTCGGGTTATGAGCCCGACG‑AGCTACCAAGCTGCTCCACCCCGCGACAAAGCTGGGGCAAAAGTATAAGACCGAACCCTGTTGAGGGTCAAATCCAAACATTCACCTACAAGAAAGCCCGCTAAAGCGGGCTCTCAAGCTTTAATTGGTACCGAGAAGGGGACTCGAAC‑‑CCCTACACCCTATGGGCACAACCACCTCAAGGTTGCGTGTCTACCAATTCCACCACCTCGGCAATACTACTACTTGAAACCCGTTACTTCTGCTCTTGA  >  NC_002947/5360928‑5361216

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