Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,391,520	T→C	10.7%	W244R (TGG→CGG) ‡	PP_2994 →	FMN‑binding oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,391,520	0	T	C	10.7%	172.6 / 9.3	65	W244R (TGG→CGG) ‡	PP_2994	FMN‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (26/32);  new base C (5/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTC  >  NC_002947/3391376‑3391663                                                                                                                                                 |                                                                                                                                                  gAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGc                                                                                                                                              <  4:369052/148‑1 (MQ=255)       cGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGa                                                                                                                                       <  1:69418/149‑1 (MQ=255)       cGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGC‑‑CACTCAACGa                                                                                                                                       <  5:274911/149‑1 (MQ=255)            gcgcCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTgg                                                                                                                                  <  5:364634/149‑1 (MQ=255)                  gTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCg                                                                                                                               >  8:240514/1‑146 (MQ=255)                     gCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACa                                                                                                                         >  5:47649/1‑149 (MQ=255)                       tCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACAtg                                                                                                                       >  1:118447/1‑149 (MQ=255)                         gAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACAtggt                                                                                                                     >  2:107155/1‑149 (MQ=255)                          aaaTCGTCAAGGCGGGCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACAtggtg                                                                                                                    <  6:267582/149‑1 (MQ=255)                           aaTCGTCAAGGCGGTCCTGGCCGTGGTTGCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCTCGGAGGGTTCACCACCGAAGACGCGAAAAAGGGCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACAtggtgg                                                                                                                   <  8:51707/149‑1 (MQ=255)                            aTCGTCAAGGCGGTCCGGTCCGTGGTTTCTGTGGATTTTGCAGTGGCGCTCAAACTCAACTCCGCCGATTTCCAGCGCGGTGGGTTCACCACCGAAGACGCGAAAAAGGGCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGa                                                                                                                  <  6:293836/149‑1 (MQ=255)                            aTCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGa                                                                                                                  <  1:107155/149‑1 (MQ=255)                            aTCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGa                                                                                                                  <  3:296578/149‑1 (MQ=255)                             tCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAg                                                                                                                 >  6:311998/1‑149 (MQ=255)                                tCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCTCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGC‑‑CACTCAACGAc                                                                                                                                      >  1:317935/1‑125 (MQ=255)                                tCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGC‑‑CACTCAACGAc                                                                                                                                      <  2:317935/125‑1 (MQ=255)                                      cGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACAtgg                                                                                                                      >  1:54830/1‑135 (MQ=255)                                      cGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACAtgg                                                                                                                      <  2:54830/135‑1 (MQ=255)                                       ggTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAg                                                                                                       <  7:33812/149‑1 (MQ=255)                                                  tGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACAtg                                                                                                                       >  1:352957/1‑122 (MQ=255)                                                  tGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACAtg                                                                                                                       <  2:352957/122‑1 (MQ=255)                                                  tGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGaa                                                                                            <  1:10225/149‑1 (MQ=255)                                                           tgGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGaa                                                                                            >  8:105268/1‑140 (MQ=255)                                                           tgGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGaa                                                                                            <  7:105268/140‑1 (MQ=255)                                                           tgGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGc                                                                                    >  2:246343/1‑148 (MQ=255)                                                           tgGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGCAGGCAGCTACGAAGTCCCCGcc                                                                                   <  7:240514/149‑1 (MQ=255)                                                              aTTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATg                                                                                >  5:187018/1‑149 (MQ=255)                                                               ttttGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGc                                                                               <  5:311998/149‑1 (MQ=255)                                                                 ttGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGACT‑‑GGCGCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCa                                                                              <  4:2056/148‑1 (MQ=255)                                                                  tgaaGTGCCGGTAAAACTCAACTCCGCCTATTTCCAGCCCGGTGGGTTAACCACCGATTACCCTAAAAAGGTCGTTTTTTCCTGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGc                                                                               <  1:113071/145‑1 (MQ=255)                                                                  tGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAg                                                                             <  2:130605/148‑1 (MQ=255)                                                                    cAGTGGCTGTCAAACTCAACTCCGGCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAATAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGt                                                                          <  7:342972/149‑1 (MQ=255)                                                                    cAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGt                                                                          <  2:118447/149‑1 (MQ=255)                                                                               aaaCTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGa                                                               <  1:246343/149‑1 (MQ=255)                                                                                    caactcCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATgg                                                            >  2:57140/1‑147 (MQ=255)                                                                                    caactcCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATgg                                                            <  1:57140/147‑1 (MQ=255)                                                                                         ccgccgATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGTTGACGCGCGCGATGGTCGCAcc                                                     >  8:58589/1‑149 (MQ=255)                                                                                            ccgATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTg                                                  <  6:47649/149‑1 (MQ=255)                                                                                               aTTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGa                                                                                                                  >  5:221765/1‑82 (MQ=255)                                                                                               aTTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGa                                                                                                                  <  6:221765/82‑1 (MQ=255)                                                                                               aTTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCt                                               <  3:213719/149‑1 (MQ=255)                                                                                                 ttCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTcg                                             >  6:297484/1‑149 (MQ=255)                                                                                                          ggAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTa                                    >  4:92254/1‑149 (MQ=255)                                                                                                              ggTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTtc                                >  4:214791/1‑149 (MQ=255)                                                                                                              ggTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTtc                                >  6:214620/1‑149 (MQ=255)                                                                                                                ttCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCAACGCCATGCAGTGTGACGCGCGCGATCGAAGCACACTGGCTCGCGAAGCTTATTacgt                              >  6:248494/1‑149 (MQ=255)                                                                                                                  caccacCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTtcgtc                             <  3:214791/148‑1 (MQ=255)                                                                                                                   accaccGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGaa                                         >  2:307204/1‑135 (MQ=255)                                                                                                                   accaccGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGaa                                         <  1:307204/135‑1 (MQ=255)                                                                                                                     caccGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGt                         <  5:214620/149‑1 (MQ=255)                                                                                                                       ccGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTc                       >  6:35773/1‑149 (MQ=255)                                                                                                                        cGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCg                      >  1:178673/1‑149 (MQ=255)                                                                                                                        cGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCg                      <  8:173901/149‑1 (MQ=255)                                                                                                                            gACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCg                  >  5:379990/1‑149 (MQ=255)                                                                                                                                cgAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGAt              >  5:140767/1‑149 (MQ=255)                                                                                                                                 gAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGAtc             <  3:92254/149‑1 (MQ=255)                                                                                                                                   aaaaGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGAtctc           >  5:283680/1‑149 (MQ=255)                                                                                                                                    ccaGGTCGTTGAGC‑‑CACTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGAtctc           >  2:77394/3‑148 (MQ=255)                                                                                                                                      aGGTCGTTGAGC‑‑CACTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTACCAGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCa         >  5:272698/1‑148 (MQ=255)                                                                                                                                        gTCGTTGAGC‑‑CACTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAAc      >  7:187372/1‑149 (MQ=255)                                                                                                                                        gTCGTTGAGC‑‑CACTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAAc      >  2:116430/1‑149 (MQ=255)                                                                                                                                         tCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAAcc     >  8:46931/1‑149 (MQ=255)                                                                                                                                            ttGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTc  >  2:362380/1‑149 (MQ=255)                                                                                                                                            ttGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTc  <  1:362380/149‑1 (MQ=255)                                                                                                                                            ttGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTc  <  2:178673/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                  GAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGT‑‑GGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTC  >  NC_002947/3391376‑3391663

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