Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,182,059	G→T	16.8%	intergenic (‑117/+37)	PP_4560 ← / ← csbD	ribonuclease BN/stress response protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,182,059	0	G	T	16.8%	188.2 / 30.2	83	intergenic (‑117/+37)	PP_4560/csbD	ribonuclease BN/stress response protein
Reads supporting (aligned to +/- strand):  ref base G (40/29);  new base T (8/6);  total (48/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.64e-01

GGCAAGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGC  >  NC_002947/5181914‑5182205                                                                                                                                                    |                                                                                                                                                                gtccaGCCGCGCAGGTCGGGGAAACTCATAGCTTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCgg                                                                                                                                                              <  4:301697/145‑1 (MQ=255) ggCAAGCCGCGCAGGTCGGGGAAAATTATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCgg                                                                                                                                                              >  1:9745/1‑149 (MQ=255) ggCAAGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCgg                                                                                                                                                              <  2:146608/149‑1 (MQ=255)  gCAAGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCggg                                                                                                                                                             >  2:287830/1‑149 (MQ=255)       ccGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAc                                                                                                                                                        >  5:332512/1‑149 (MQ=255)        cgcgCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCgg                                                                                                                                                              >  8:85324/1‑142 (MQ=255)        cgcgCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCgg                                                                                                                                                              <  7:85324/142‑1 (MQ=255)         gcgcAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGAGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAccc                                                                                                                                                      >  1:194494/1‑149 (MQ=255)            cAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCgag                                                                                                                                                   <  6:332512/149‑1 (MQ=255)              ggTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGcgaa                                                                                                                                                 <  4:74546/149‑1 (MQ=255)              ggTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGagaa                                                                                                                                                 <  3:115327/149‑1 (MQ=255)              ggTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGGCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGagaa                                                                                                                                                 <  2:263903/149‑1 (MQ=255)               gTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCgag                                                                                                                                                   >  4:121007/1‑146 (MQ=255)               gTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCgag                                                                                                                                                   >  3:67753/1‑146 (MQ=255)                tCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGcgaaga                                                                                                                                               <  3:328955/149‑1 (MQ=255)                   gggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAccc                                                                                                                                                      >  1:77148/1‑139 (MQ=255)                   gggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAccc                                                                                                                                                      >  1:82670/1‑139 (MQ=255)                   gggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAccc                                                                                                                                                      <  2:82670/139‑1 (MQ=255)                   gggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAccc                                                                                                                                                      <  2:77148/139‑1 (MQ=255)                    ggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑c                                                                                                                                           <  1:287830/149‑1 (MQ=255)                             tAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGc                                                                                                                                  >  8:86918/1‑148 (MQ=255)                                 gTCTCCTTTCGCCGCTTGATCTTAACACCGCAACGTATAATACACCGGCCCCGCGCTGTGCTGCTCCCGGCCGTGGTTT‑AAC‑AATTTGCTTACCAAAAAGGGCCGGCCTCTTCGCGGGTCCACCCTAGAAGTGG‑CCGGCACCGCCtcaa                                                                                                                             <  2:9745/149‑1 (MQ=255)                                       tttCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAAGAGG‑CCgg                                                                                                                                        <  5:218634/132‑1 (MQ=255)                                       tttCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAAGAGG‑CCgg                                                                                                                                        >  6:218634/1‑132 (MQ=255)                                                   aTCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAC‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCAcc                                                                                                            <  5:345190/148‑1 (MQ=255)                                                       ttACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATc                                                                                                       >  3:25062/1‑149 (MQ=255)                                                           acaGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTc                                                                                                   >  5:214011/1‑149 (MQ=255)                                                           acaGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CAGGCACCGCCTCAAGTCAAGGCTTCTTCACCCCATCCtt                                                                                                    >  5:208858/1‑148 (MQ=255)                                                             aGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAGAAATTTGCTTACGAATAAGTGCCGGCCTATTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGATTCTTCACCGCATCCTTc                                                                                                   >  1:103558/1‑148 (MQ=255)                                                             aGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCAACCTTca                                                                                                  >  2:171003/1‑148 (MQ=255)                                                                    tCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGTTTG‑AAG‑AATTTGCTTACTAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCCGGCACCGCCTCAAGTCAAGGCTTCTTCAc                                                                                                             <  1:225781/131‑1 (MQ=255)                                                                    tCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCACGTCAAGGCTTCTTCACCGCATCCTTCACATTTccc                                                                                          >  3:258511/1‑149 (MQ=255)                                                                      tAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCtt                                                                                        >  6:23330/1‑149 (MQ=255)                                                                           gACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑Cct                                                                         >  3:283038/1‑149 (MQ=255)                                                                              cGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAagag                                                                                                                                              >  4:72113/1‑88 (MQ=37)                                                                              cGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAagag                                                                                                                                              <  3:72113/88‑1 (MQ=37)                                                                              cGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑Cctgct                                                                      >  7:124131/1‑149 (MQ=255)                                                                                      aGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTACTGcg‑c‑ctcg                                                            >  3:66555/1‑149 (MQ=255)                                                                                          gtgCTGCTCCCGGCCGTGGTTG‑ATG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCtt                                                        >  4:50948/1‑149 (MQ=255)                                                                                           tgctgcTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCttt                                                       >  7:104055/1‑149 (MQ=255)                                                                                            gctgctCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGcg‑c‑ctc                                                             >  1:39173/1‑142 (MQ=255)                                                                                            gctgctCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGcg‑c‑ctc                                                             <  2:39173/142‑1 (MQ=255)                                                                                             ctgctCCCGGCCGTGGTTG‑ATG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTg                                                      <  3:50948/148‑1 (MQ=255)                                                                                             ctgctCCCGGCCGTGGTTC‑AAT‑ATTTTGCTTACGAATACGGGCCGGCCTCTTCGCGGGTGCCCCCGAGAAGTGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGa                                                     <  4:67753/149‑1 (MQ=255)                                                                                              tgctCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAg                                                    >  8:203652/1‑149 (MQ=255)                                                                                                 tCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCg                                                                                                                                         <  4:319652/73‑1 (MQ=255)                                                                                                 tCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCg                                                                                                                                         >  3:319652/1‑73 (MQ=255)                                                                                                 tCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTc                                                 >  7:103334/1‑149 (MQ=255)                                                                                                   ccGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCt                                               >  7:209023/1‑149 (MQ=255)                                                                                                     ggCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTgc                                             <  7:86918/149‑1 (MQ=255)                                                                                                        cGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCCTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCgcc                                          >  2:283126/1‑149 (MQ=255)                                                                                                               g‑aag‑aaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTGCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑ATCGCCTTTGAGCTCCTGAgcctttcctt                                   >  7:129235/1‑149 (MQ=255)                                                                                                               g‑aag‑aaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCgccttgcctt                                   >  7:134404/1‑149 (MQ=255)                                                                                                                 aag‑aaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCgccttgcctt                                   <  4:283038/148‑1 (MQ=255)                                                                                                                   acaaTTTCCTAACGAATACGGGCCGGCCTCTCCGCGGGTCCACCCGACAAGCGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCg                                 <  8:55574/147‑1 (MQ=255)                                                                                                                              aCGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGc                        <  1:105196/147‑1 (MQ=255)                                                                                                                                 aaTAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                   <  6:214011/149‑1 (MQ=255)                                                                                                                                 aaTAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                   <  3:121007/149‑1 (MQ=255)                                                                                                                                  atcctggCCGGCCTCTTCGCGGGTTCACCGTATAAGTGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                   <  6:208858/142‑1 (MQ=255)                                                                                                                                  aTAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCg                  >  6:150598/1‑149 (MQ=255)                                                                                                                                     aGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTg               <  3:281633/149‑1 (MQ=255)                                                                                                                                       tGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTc             <  2:103558/149‑1 (MQ=255)                                                                                                                                        gCCGGCCTCTTCTCGGGTGAACCCGCGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCg            >  3:280100/1‑149 (MQ=255)                                                                                                                                        gCCGGCCTCTGCGCGGGTCCCCCCCACATGTGC‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCg            <  1:171003/149‑1 (MQ=255)                                                                                                                                         ccGGCCTCTTCGCGGGTTCACCCGAGAATTGG‑CCGGCACCGCCTCAACTCAATGCTTCTTCACCGCATCCTTCACATTTACCCTTA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCC‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGt              <  1:283126/147‑1 (MQ=255)                                                                                                                                         ccGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCgg           >  3:366756/1‑149 (MQ=255)                                                                                                                                         ccGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCg            <  4:25062/148‑1 (MQ=255)                                                                                                                                          cGGCCTCTTCTCGGGTGAACCCGCGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGtt                <  1:280321/143‑1 (MQ=255)                                                                                                                                          cGGCCTCTTCTCGGGTGAACCCGCGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGtt                >  2:280321/1‑143 (MQ=255)                                                                                                                                          cGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGt          >  7:143919/1‑149 (MQ=255)                                                                                                                                           ggCCTCTTCTCGGGTGAACCCGCGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGt          >  7:345754/1‑148 (MQ=255)                                                                                                                                           ggCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CATGCTGCG‑C‑ATCGCCTTTGAGCTCCTGCGCCTTGACTTCGGCCCGCAGCTTTTCTTTGTCTGt          >  8:256212/1‑148 (MQ=255)                                                                                                                                            gCCTCTTCTCGGGTGAACCCGCGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCGTTTCGTTGTAGGTCa        >  1:360316/1‑149 (MQ=39)                                                                                                                                               tctTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTgcg                                            <  2:320823/110‑1 (MQ=255)                                                                                                                                               tctTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTgcg                                            >  1:320823/1‑110 (MQ=255)                                                                                                                                               tctTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCt     >  7:66302/1‑149 (MQ=255)                                                                                                                                               tctTCGCGGGTTCACCCGAGAAGAGG‑CCGCCACCGCCACAAGTCAAGGCTTCTTCACCGCCTCCTTCACATTTACCATGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCT‑C‑ATCGCATTAGAACTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCt     >  2:186581/1‑149 (MQ=255)                                                                                                                                                ctTCTCGGGTGAACCCGCGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTgc                                             <  8:104402/108‑1 (MQ=38)                                                                                                                                                ctTCTCGGGTGAACCCGCGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTgc                                             >  7:104402/1‑108 (MQ=38)                                                                                                                                                 ttCTCGGGTGAACCCGCGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTg   >  6:360708/1‑149 (MQ=255)                                                                                                                                                 ttCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTg   >  4:128061/1‑149 (MQ=255)                                                                                                                                                 ttCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTg   >  7:203200/1‑149 (MQ=255)                                                                                                                                                 ttCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTg   >  1:6700/1‑149 (MQ=255)                                                                                                                                                   cgcgGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGc  <  8:103334/148‑1 (MQ=255)                                                                                                                                                    gcgGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGAAGGTGAAGGGGATGATGAG‑CACACCACTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCgg           >  1:122949/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                                GGCAAGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTG‑AAG‑AATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGG‑CCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCG‑C‑CTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGC  >  NC_002947/5181914‑5182205

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