Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,388,336	C→G	11.2%	intergenic (+73/+71)	pyrD → / ← rmf	dihydroorotate dehydrogenase/ribosome modulation factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,388,336	0	C	G	11.2%	176.9 / 10.5	62	intergenic (+73/+71)	pyrD/rmf	dihydroorotate dehydrogenase/ribosome modulation factor
Reads supporting (aligned to +/- strand):  ref base C (29/26);  new base G (3/4);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACCGC  >  NC_002947/2388189‑2388480                                                                                                                                                      |                                                                                                                                                  taCTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTc                                                                                                                                                   <  8:141244/148‑1 (MQ=255)  aCTCGGGCTTCATCTACAAAGGGCCGGAGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTc                                                                                                                                                   >  1:125788/1‑147 (MQ=255)   cTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCaa                                                                                                                                                 <  1:66253/148‑1 (MQ=255)     cGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGCTGGATGCTATCGCGGCAATGCCGCGGTTACCAGACCGCCCATGT‑GGGGGGGCGGCTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAggg                                                                                                                                              >  7:77592/1‑149 (MQ=255)      gggCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAgg                                                                                                                                               <  2:109933/147‑1 (MQ=255)       cgattCATATACAAAGGGCCGGGGCTGATTCGCGAAGGGGTGGATTCTATCGCGGCAAAGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGGCGGGTGGG‑TGCAGCTTATGAGGTTTGCCAATTCATAAAAAAGGGCCCCTTGAAgggg                                                                                                                                             <  2:62978/146‑1 (MQ=255)             accTACAACGGGCGGCCGC‑GCTTCCCGAACCGGGGGCAGCTATCGCGGCAAAGCCGCGGTGACCACCCCGCCCCTGG‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTGAAGGGGcccttttt                                                                                                                                     <  4:318165/147‑6 (MQ=255)             aTCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTGAAgggg                                                                                                                                             <  7:255787/142‑1 (MQ=255)             aTCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTGAAgggg                                                                                                                                             >  8:255787/1‑142 (MQ=255)             aTCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                                                                                                      <  2:284415/149‑1 (MQ=255)              tCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTggg                                                                                                                                     >  1:10717/1‑149 (MQ=255)                   aaaGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑GGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTc                                                                                                                                >  5:268140/1‑149 (MQ=255)                        gCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTcgcccg                                                                                                                           >  4:161434/1‑149 (MQ=255)                            gcgcTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGACATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCcg‑ccgc                                                                                                                      >  7:131133/1‑149 (MQ=255)                                   ttCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGACATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGAt                                                                                                                <  8:131133/148‑1 (MQ=255)                                     cgcgAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACGAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATgg                                                                                                              <  3:18931/148‑1 (MQ=255)                                      gcgAAGCGGTGTATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTTT‑GGGGGGGCGGGTGGG‑GGCCGCTGATGAGGTTTGCCAACGCATAAAAAAGTGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATggg                                                                                                             <  5:11143/148‑1 (MQ=255)                                      gcgAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATgggg                                                                                                            >  5:368535/1‑149 (MQ=255)                                      gcgAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATggg                                                                                                             <  7:205289/148‑1 (MQ=255)                                       cgAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGc                                                                                                           <  2:125788/149‑1 (MQ=255)                                                    gCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCccg‑ccg                                                                                                                       >  7:140145/1‑124 (MQ=255)                                                    gCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCccg‑ccg                                                                                                                       <  8:140145/124‑1 (MQ=255)                                                     tcaccgcgGCATTGCCGCGATGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATc                                                                                             >  4:135513/5‑149 (MQ=255)                                                     cTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATc                                                                                             <  2:10717/149‑1 (MQ=255)                                                      caccgcgGCATTGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATc                                                                                             >  2:107961/4‑148 (MQ=255)                                                        tCGCTGCATTGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGaa                                                                                          <  3:203466/149‑1 (MQ=255)                                                        tCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGaa                                                                                          >  3:69027/1‑149 (MQ=255)                                                               aaTGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAgtg                                                                                   >  2:50956/1‑149 (MQ=255)                                                                     gcgGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGGATGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTc                                                                              >  1:258394/1‑149 (MQ=255)                                                                         tGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCa                                                                         <  6:368535/149‑1 (MQ=255)                                                                          gACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAg                                                                        >  8:264558/1‑149 (MQ=255)                                                                             cAGCCCGCCCCTGT‑‑GGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTGaa                                                                                                                                                 >  2:153648/1‑73 (MQ=255)                                                                             cAGCCCGCCCCTGT‑‑GGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTGaa                                                                                                                                                 <  1:153648/73‑1 (MQ=255)                                                                               gcccgcccCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAAc                                                                   <  3:161434/149‑1 (MQ=255)                                                                                     cccTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGTCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGt                                                              <  1:107961/149‑1 (MQ=255)                                                                                     cccTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGT‑CTCGTTCAGCCAACGGCGTGa                                                            >  3:220763/1‑149 (MQ=255)                                                                                               ggggCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGtt                                                    <  1:50956/149‑1 (MQ=255)                                                                                                 ggCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGa                                                  <  3:135513/149‑1 (MQ=255)                                                                                                 ggCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGa                                                  <  5:148217/149‑1 (MQ=255)                                                                                                      gtggg‑tgCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTGAGCCAACg                                                                  >  5:288584/1‑128 (MQ=255)                                                                                                       tggg‑tgCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCTCCCGGATGGGGGGTGCATGGTAAATCTAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTAt                                            <  8:77592/149‑1 (MQ=255)                                                                                                        ggg‑tgCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATg                                           <  4:69027/149‑1 (MQ=255)                                                                                                            tgCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGAt                                        <  6:39662/149‑1 (MQ=255)                                                                                                                   gatgaGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGTATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCg                                 >  2:44484/1‑149 (MQ=255)                                                                                                                   gatgaGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCg                                 >  7:261459/1‑149 (MQ=255)                                                                                                                      gaGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGt                               <  7:264558/148‑1 (MQ=255)                                                                                                                        ggTTTGCCCAGGCATAAAAAAGGCCCCCTTCAAGGGGCCGCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGc                             <  6:268140/148‑1 (MQ=255)                                                                                                                         gTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGc                             >  4:144264/1‑147 (MQ=255)                                                                                                                          tttGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAg                          <  2:258394/149‑1 (MQ=255)                                                                                                                               cAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCAtt                     >  6:11050/1‑149 (MQ=255)                                                                                                                               cAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCAtt                     >  1:196415/1‑149 (MQ=255)                                                                                                                               cAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCAtt                     >  2:356476/1‑149 (MQ=255)                                                                                                                                aaGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCa                                                                         <  7:158369/96‑1 (MQ=255)                                                                                                                                aaGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCa                                                                         >  8:158369/1‑96 (MQ=255)                                                                                                                                   gCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCTTCTCATTCCGt                 >  1:319516/1‑149 (MQ=255)                                                                                                                                    cATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTc                >  2:315896/1‑149 (MQ=255)                                                                                                                                     aTAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTcc               >  5:371489/1‑149 (MQ=255)                                                                                                                                          aaaGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGt           >  2:355859/1‑148 (MQ=255)                                                                                                                                            ggggCCCCTTGAAGGGGCCCTTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTa         >  6:212279/2‑148 (MQ=255)                                                                                                                                            aGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCATCCAAGTTAt        >  2:364393/1‑149 (MQ=255)                                                                                                                                               gCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCAc     >  7:211270/1‑149 (MQ=255)                                                                                                                                               gCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGT‑CTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCAcc    <  4:220763/149‑1 (MQ=255)                                                                                                                                                 cccTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACcg   >  2:140933/1‑149 (MQ=255)                                                                                                                                                  ccTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACcgc  <  2:319516/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  TACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGT‑GGGGGGGCGGGTGGG‑TGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCG‑CCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACCGC  >  NC_002947/2388189‑2388480

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