Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,886,602	C→T	7.2%	D423D (GAC→GAT)	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,602	0	C	T	7.2%	203.9 / 4.5	69	D423D (GAC→GAT)	rmuC	DNA recombination limiting inversions protein
Reads supporting (aligned to +/- strand):  ref base C (30/34);  new base T (3/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

GATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGA  >  NC_002947/1886455‑1886745                                                                                                                                                     |                                                                                                                                                 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                                                                                                                                                  >  6:87766/1‑136 (MQ=255)             tACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGAc                                                                                                                                                  <  5:87766/136‑1 (MQ=255)                cACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCGTCCAGGTCCTGGATGAAca                                                                                                                                             <  3:159077/138‑2 (MQ=255)                cACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCGTCCAGGTCCTGGATGAAca                                                                                                                                             >  4:159077/1‑137 (MQ=255)                 aCGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCg                                                                                                                                    >  5:15745/1‑146 (MQ=255)                 aCGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCg                                                                                                                                    <  6:15745/146‑1 (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                                                                                                                                             >  4:244516/1‑95 (MQ=255)                                                          gCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCGTCCAGGTCCTGGATGAAca                                                                                                                                             <  3:244516/96‑2 (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:233090/1‑149 (MQ=255)                                                                                   cgcgAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGt                                                              >  6:237633/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                                                                                  >  5:159992/1‑95 (MQ=255)                                                                                                                      aCAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGc                                                                                  <  6:159992/95‑1 (MQ=255)                                                                                                                       cAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGccc                           >  1:273399/1‑149 (MQ=255)                                                                                                                        aaGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCg                          >  4:78527/1‑149 (MQ=255)                                                                                                                          gTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGgc                               >  7:179351/1‑142 (MQ=255)                                                                                                                          gTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGgc                               >  7:184153/1‑142 (MQ=255)                                                                                                                          gTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGgc                               <  8:184153/142‑1 (MQ=255)                                                                                                                          gTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGgc                               <  8:179351/142‑1 (MQ=255)                                                                                                                             cGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAgctgc                                    >  1:185603/1‑134 (MQ=255)                                                                                                                             cGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAgctgc                                    <  2:185603/134‑1 (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)                                                                                                                                               ccTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCa    <  4:173662/148‑1 (MQ=255)                                                                                                                                                cTGGACGAACTGGGCAGCCGCTTGCAGTAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGa  >  6:159296/1‑149 (MQ=255)                                                                                                                                                   gACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCt            <  7:173894/136‑1 (MQ=255)                                                                                                                                                   gACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCt            >  8:173894/1‑136 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGC‑GGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCA‑ACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGA  >  NC_002947/1886455‑1886745

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