Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	532,945	C→G	6.1%	Y70* (TAC→TAG)	tufA →	elongation factor Tu‑A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	532,945	0	C	G	6.1%	230.9 / 5.0	81	Y70* (TAC→TAG)	tufA	elongation factor Tu‑A
Reads supporting (aligned to +/- strand):  ref base C (40/36);  new base G (2/3);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GACCACGGTAAGACCACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACAAACCCGTGAGCA  >  NC_002947/532799‑533091                                                                                                                                                    |                                                                                                                                                   gACCACGGTAAGGCCACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGCCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCCA‑GTACaa                                                                                                                                                  >  3:301793/1‑149 (MQ=255)  aCCACGGTAAGACCACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACaa                                                                                                                                                  <  4:363175/148‑1 (MQ=34)   ccACGGTAAGACCACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACt                                                                                                                                                <  1:61905/149‑1 (MQ=34)        gTAAGACCACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCTAAGTTTTCTGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTTGTACTCGAcg                                                                                                                                           <  6:210316/149‑2 (MQ=11)        gTAAGACCACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGAc                                                                                                                                            >  5:339336/1‑148 (MQ=35)         tAAGACCACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTTGTACTCGAcgt                                                                                                                                          <  2:97041/149‑3 (MQ=255)            gACCACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCGAGGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTc                                                                                                                                       <  5:268576/149‑1 (MQ=37)                acgctGCCCGGAGCTCTGACTCGCGTTTGCTCCGAAGGTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACAGCGCTCAGGTCGA‑GTAGAACTCGACCATTCGTCa                                                                                                                                   <  4:301793/146‑1 (MQ=11)                        cGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATAGAATAGGCTACGGAAGAAAAAGCGCGCGGTATCACCATCAACAACGCTCACGTCGA‑ATACAACTAGACCATTCGTCACTACgat                                                                                                                            >  5:295473/1‑146 (MQ=11)                           aGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACg                                                                                                                        <  7:44056/149‑1 (MQ=37)                           aGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACg                                                                                                                        <  3:147232/149‑1 (MQ=37)                                  aCTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAATTCGACCATTCGTCACTACGCTCACGTTGa                                                                                                                    <  1:362591/146‑1 (MQ=37)                                  aCTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACt                                                                                                                  <  5:194655/148‑1 (MQ=37)                                   cTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTg                                                                                                                 <  6:339336/148‑1 (MQ=37)                                    tCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAACAGCGCGCGGTATCACCATCAACACCTCTCACGTCGA‑GTACAACTAGACCATTCGTCACTACGCTCACGTTGCCGGcc                                                                                                               >  2:253314/1‑149 (MQ=11)                                    tCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGcc                                                                                                               >  7:125068/1‑149 (MQ=37)                                      gcgTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCa                                                                                                             <  6:265058/149‑1 (MQ=37)                                           tGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACTTTGACTGCCCAGGTCa                                                                                                        >  8:116679/1‑149 (MQ=21)                                                 gAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGa                                                                                                  >  7:78843/1‑149 (MQ=34)                                                      tttCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTAc                                                                                              <  7:195949/148‑1 (MQ=34)                                                        tcggttcggCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGt                                                                                            <  5:150615/148‑1 (MQ=34)                                                          ggttcggCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTgaa                                                                                         <  7:116679/149‑1 (MQ=34)                                                          ggttcggCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTgaa                                                                                         <  7:99552/149‑1 (MQ=34)                                                           gttcggCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTg                                                                                           >  5:309456/1‑146 (MQ=34)                                                           gttcggCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTg                                                                                           >  6:351764/1‑146 (MQ=34)                                                            ttcggCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTgaaga                                                                                       <  8:125068/149‑1 (MQ=34)                                                               ggCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTc                                                                                                                                    <  5:191977/101‑1 (MQ=35)                                                               ggCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTc                                                                                                                                    >  6:191977/1‑101 (MQ=35)                                                                  aaTCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGa                                                                                 >  7:211333/1‑149 (MQ=33)                                                                   aTCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGAt                                                                                <  4:101409/149‑1 (MQ=32)                                                                         gAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCgg                                                                          >  2:131766/1‑149 (MQ=32)                                                                          aGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGa                                                                                                  <  7:224301/124‑1 (MQ=32)                                                                          aGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGa                                                                                                  >  8:224301/1‑124 (MQ=32)                                                                                aCAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGcc                                                                   >  1:114520/1‑149 (MQ=32)                                                                                aCAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGcc                                                                   <  5:351764/149‑1 (MQ=32)                                                                                 cAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGccc                                                                  <  8:78843/149‑1 (MQ=32)                                                                                      tCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGAt                                                              <  3:9899/148‑1 (MQ=32)                                                                                       cGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATg                                                             >  7:262234/1‑148 (MQ=32)                                                                                        gactcgactcCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGa                                                           >  5:354075/1‑149 (MQ=17)                                                                                        gACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGa                                                           >  5:361322/1‑149 (MQ=32)                                                                                          cTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGcc                                                                   >  4:144365/1‑139 (MQ=32)                                                                                          cTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGcc                                                                   <  3:144365/139‑1 (MQ=32)                                                                                             ggCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGgcg                                                      <  2:114520/149‑1 (MQ=32)                                                                                             ggCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACAGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGgcg                                                      >  7:258383/1‑149 (MQ=17)                                                                                              gCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGCTGGACGgcgc                                                     <  1:253314/149‑1 (MQ=17)                                                                                              gCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGgcgc                                                     >  6:218726/1‑149 (MQ=32)                                                                                              gCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGgcg                                                      >  7:37008/1‑148 (MQ=32)                                                                                              gCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGgcg                                                      <  8:37008/148‑1 (MQ=32)                                                                                                  cGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGGCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGa                                                   <  4:98962/147‑1 (MQ=17)                                                                                                    gaagaaAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACg                                                         <  5:154703/139‑1 (MQ=32)                                                                                                    gaagaaAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACg                                                         >  6:154703/1‑139 (MQ=32)                                                                                                    gaagaaAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCt                                               <  8:262234/149‑1 (MQ=32)                                                                                                       gaaAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGt                                            >  4:222897/1‑149 (MQ=32)                                                                                                       gaaAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGt                                            <  3:222897/149‑1 (MQ=32)                                                                                                        aaaaaGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCGCGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGtt                                           >  2:205781/1‑149 (MQ=17)                                                                                                        aaaaaGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGtt                                           >  8:332364/1‑149 (MQ=32)                                                                                                             gcgcgcGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTc                                      <  7:332364/149‑1 (MQ=32)                                                                                                               gcgcgGTATCACCATCAACACCGCTCACGTCGAGGTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCg                                     >  1:251285/1‑149 (MQ=17)                                                                                                                 gcgGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGgcc                                  <  6:361322/149‑1 (MQ=32)                                                                                                                     tATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGgccgcc                               >  5:272708/1‑148 (MQ=32)                                                                                                                     tATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGgccgcc                               >  7:58545/1‑148 (MQ=32)                                                                                                                          ccATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGt                         >  7:229434/1‑149 (MQ=32)                                                                                                                           cATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTc                        >  6:285939/1‑149 (MQ=32)                                                                                                                           cATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGATCGGCCGCCGATGGTc                        >  7:91356/1‑149 (MQ=17)                                                                                                                            aTCAAAACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTcc                       >  7:249216/1‑149 (MQ=17)                                                                                                                             tCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTcc                       >  3:139586/1‑148 (MQ=32)                                                                                                                               aaCACCGCTCACGTCGA‑GTACAACTCGACGATTCTTCACTAGGCTCACGTGTACTGCCCACGTGACGCTGACTACGTGAAGAACATTATCACCGGTGCTGCCCATATTGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGa                     <  6:309456/148‑1 (MQ=11)                                                                                                                                acacCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATg                   >  6:75933/1‑149 (MQ=32)                                                                                                                                  acCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGcc                 <  5:63515/149‑1 (MQ=32)                                                                                                                                      gaatggTCGA‑GTTGTACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACaa             >  8:205112/6‑149 (MQ=255)                                                                                                                                         tggTCGA‑GTTGTACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACAAAcc          >  2:48238/3‑149 (MQ=255)                                                                                                                                         aCGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACaaa            <  5:285939/147‑1 (MQ=34)                                                                                                                                         aCGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACAAAcc          >  6:161571/1‑149 (MQ=34)                                                                                                                                         aCGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACAAAcc          >  3:342829/1‑149 (MQ=34)                                                                                                                                           gTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACAAACCCg        >  6:30733/1‑149 (MQ=34)                                                                                                                                            tCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACAAACCCGt       <  8:211333/149‑1 (MQ=34)                                                                                                                                              gA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACAAACCCGTGa     <  6:272708/149‑1 (MQ=34)                                                                                                                                               a‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACAAACCCGTGAg    >  3:199834/1‑149 (MQ=34)                                                                                                                                                 gTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCg                      <  4:36064/130‑1 (MQ=32)                                                                                                                                                 gTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCg                      >  3:36064/1‑130 (MQ=32)                                                                                                                                                  tACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGCTCACGCTGACTACGTGAAGAACATGATCACCGGTGTTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACAAACCCGTGAGca  >  4:188326/1‑149 (MQ=21)                                                                                                                                                    |                                                                                                                                                   GACCACGGTAAGACCACTCTGACCGCAGCTCTGACTCGCGTTTGCTCCGAAGTTTTCGGTTCGGCAATCGTTGAGTTCGACAAGATCGACTCGGCTCCGGAAGAAAAAGCGCGCGGTATCACCATCAACACCGCTCACGTCGA‑GTACAACTCGACCATTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCACAAACCCGTGAGCA  >  NC_002947/532799‑533091

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