Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,886,595	A→T	7.2%	D421V (GAC→GTC)	rmuC →	DNA recombination limiting inversions protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,886,595	0	A	T	7.2%	200.3 / 5.8	69	D421V (GAC→GTC)	rmuC	DNA recombination limiting inversions protein
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base T (3/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

ATCGTGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCA  >  NC_002947/1886450‑1886743                                                                                                                                                   |                                                                                                                                                      aTCGTGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCt                                                                                                                                                    >  5:246428/1‑149 (MQ=255)  tCGTGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTg                                                                                                                                                   <  1:239641/149‑1 (MQ=255)      gATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGACTCTACGACAAGTTCGTGCTGTTAATCCAGGACCTGGAc                                                                                                                                                >  1:269515/1‑148 (MQ=255)           tCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTg                                                                                                                                          >  4:203117/1‑149 (MQ=255)             aGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTggg                                                                                                                                        >  3:112635/1‑149 (MQ=255)             aGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTggg                                                                                                                                        <  2:136062/149‑1 (MQ=255)              gCCCTACCACGTTGCTGGCCACCTTTCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGACCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTggg                                                                                                                                        >  4:161470/1‑148 (MQ=255)               cccTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTgg                                                                                                                                         <  6:246428/146‑1 (MQ=255)               cccTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCa                                                                                                                                      >  4:126975/1‑149 (MQ=255)                  tACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGAc                                                                                                                                                <  5:87766/136‑1 (MQ=255)                  tACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGAc                                                                                                                                                >  6:87766/1‑136 (MQ=255)                     cACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCGTCCAGGTCCTGGATGAAca                                                                                                                                           <  3:159077/138‑2 (MQ=255)                     cACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCGTCCAGGTCCTGGATGAAca                                                                                                                                           >  4:159077/1‑137 (MQ=255)                      aCGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCg                                                                                                                                  <  6:15745/146‑1 (MQ=255)                      aCGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCg                                                                                                                                  >  5:15745/1‑146 (MQ=255)                       cGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTg                                                                                                                              <  2:269515/149‑1 (MQ=255)                         ttGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTgca                                                                                                                            <  4:112635/149‑1 (MQ=255)                                ccACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGACCGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTg                                                                                                                     >  7:50388/1‑149 (MQ=255)                                ccACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTg                                                                                                                     <  1:319484/149‑1 (MQ=255)                                    cTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGAc                                                                                                                  >  1:27885/1‑148 (MQ=255)                                              aTCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATg                                                                                                        <  2:27885/148‑1 (MQ=255)                                                     gCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAgcg                                                                                                >  5:18670/1‑149 (MQ=255)                                                          tgGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGTCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAgcgcgc                                                                                             <  3:60100/147‑1 (MQ=255)                                                            gAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGca‑a                                                                                          >  1:229046/1‑148 (MQ=255)                                                              aGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAg                                                                                      <  6:30056/149‑1 (MQ=255)                                                               gCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCGTCCAGGTCCTGGATGAAca                                                                                                                                           <  3:244516/96‑2 (MQ=255)                                                               gCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCGTCCAGGTCCTGGATGAAca                                                                                                                                           >  4:244516/1‑95 (MQ=255)                                                                   gAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAg                                                                                      <  7:52578/144‑1 (MQ=255)                                                                   gAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGATGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAg                                                                                      >  8:52578/1‑144 (MQ=255)                                                                       gTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGaa                                                                             <  5:47694/149‑1 (MQ=255)                                                                        tCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGAACTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAg                                                                            <  3:203117/149‑1 (MQ=255)                                                                          aaGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAggg                                                                          >  7:197236/1‑149 (MQ=255)                                                                           aGGGCAGCACGCCCGCGAACATGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGc                                                                         <  8:266326/149‑1 (MQ=255)                                                                              gCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTg                                                                      >  7:137956/1‑149 (MQ=255)                                                                               cAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTgg                                                                     >  4:39079/1‑149 (MQ=255)                                                                                 gAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGc                                                                    >  2:231964/1‑148 (MQ=255)                                                                                        cgcgAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGt                                                            >  8:251240/1‑149 (MQ=255)                                                                                        cgcgAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGt                                                            >  6:237633/1‑149 (MQ=255)                                                                                        cgcgAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGt                                                            >  6:233090/1‑149 (MQ=255)                                                                                          cgAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTa                                                          <  8:50388/149‑1 (MQ=255)                                                                                            aaaTTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAAc                                                                 >  1:158059/1‑140 (MQ=255)                                                                                            aaaTTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAAc                                                                 <  2:158059/140‑1 (MQ=255)                                                                                            aaaTTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGc                                                        <  1:231964/149‑1 (MQ=255)                                                                                            aaaTTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGc                                                        <  6:18670/149‑1 (MQ=255)                                                                                             aaTTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGcc                                                       <  3:161470/149‑1 (MQ=255)                                                                                                  ccGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGc                                                  <  3:126975/149‑1 (MQ=255)                                                                                                                      ctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAgctgctcg                               <  3:39079/149‑3 (MQ=255)                                                                                                                          gACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTggg                              <  2:8150/146‑1 (MQ=255)                                                                                                                          gACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTggg                              >  1:8150/1‑146 (MQ=255)                                                                                                                          gACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGgcgc                           <  7:251240/149‑1 (MQ=255)                                                                                                                          gACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGgcgc                           <  2:229046/149‑1 (MQ=255)                                                                                                                           aCAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGc                                                                                <  6:159992/95‑1 (MQ=255)                                                                                                                           aCAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGc                                                                                >  5:159992/1‑95 (MQ=255)                                                                                                                            cAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGccc                         >  1:273399/1‑149 (MQ=255)                                                                                                                             aaGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCg                        >  4:78527/1‑149 (MQ=255)                                                                                                                               gTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGgc                             <  8:184153/142‑1 (MQ=255)                                                                                                                               gTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGgc                             <  8:179351/142‑1 (MQ=255)                                                                                                                               gTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGgc                             >  7:184153/1‑142 (MQ=255)                                                                                                                               gTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGgc                             >  7:179351/1‑142 (MQ=255)                                                                                                                                  cGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAgctgc                                  <  2:185603/134‑1 (MQ=255)                                                                                                                                  cGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAgctgc                                  >  1:185603/1‑134 (MQ=255)                                                                                                                                     ccagTTCGTCCAGGTCCTGGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTCGCAACCTGGTTAGCCTTAGCGAGCAGCTTAAGCTGATGGGCGCCCGTGCCAGc                 >  3:164294/4‑148 (MQ=255)                                                                                                                                         ttCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAaga             <  2:273399/148‑1 (MQ=255)                                                                                                                                              ccAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAGACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTg        >  1:258148/1‑149 (MQ=255)                                                                                                                                              ccAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGc       >  3:173662/1‑149 (MQ=255)                                                                                                                                               cAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGcc      >  2:13250/1‑149 (MQ=255)                                                                                                                                               cAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGc       <  7:300535/148‑1 (MQ=255)                                                                                                                                                aGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCg     <  5:237633/149‑1 (MQ=255)                                                                                                                                                aGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCg     <  5:233090/149‑1 (MQ=255)                                                                                                                                                   aCCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCGGCAAGCTCTTGCTGGCa  <  4:236961/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      ATCGTGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCA  >  NC_002947/1886450‑1886743

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