Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,968,146	C→A	9.6%	intergenic (‑42/+33)	PP_4377 ← / ← fliC	flagellin FlaG/flagellin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,968,146	0	C	A	9.6%	217.6 / 10.6	73	intergenic (‑42/+33)	PP_4377/fliC	flagellin FlaG/flagellin
Reads supporting (aligned to +/- strand):  ref base C (29/37);  new base A (2/5);  total (31/42)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

ACGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGGCAGCAAAGT  >  NC_002947/4968001‑4968294                                                                                                                                                   |                                                                                                                                                      aCGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCg                                                                                                                                                     >  7:47666/1‑148 (MQ=255)    tCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCAcc                                                                                                                                                 <  7:168686/149‑1 (MQ=255)       ggTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCtt                                                                                                                                              >  8:181787/1‑149 (MQ=255)         tAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTGGAATCAGAAAACGAAAGGGTGCGGCAACGTACCCTTTc                                                                                                                                            <  6:15788/149‑1 (MQ=255)         tAGCTCAGCATTCGATCTGCTCGATGGCAGTGCACTCCTGGCGTCTTGAGTTGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTGCGGCAACGTACCCTTTc                                                                                                                                            <  6:315100/149‑1 (MQ=255)         tAGCTCAGCATTCGATCTGATCGATGGCTGTGCACTCCTGTCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTGCGGCAACGTACCCTTTc                                                                                                                                            <  1:244041/149‑1 (MQ=255)              cAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTGCGGCAACGTACCCTTTcgtttt                                                                                                                                       >  4:139743/1‑145 (MQ=255)               aGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTGCGGCAACGTACCCTTTcgttttc                                                                                                                                      <  3:30322/149‑6 (MQ=255)               aGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGa                                                                                                                                        <  5:293825/147‑1 (MQ=255)                gCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTGCGGCAACGTACCCTTTcgttttct                                                                                                                                     <  3:139743/149‑7 (MQ=255)                gCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATtc                                                                                                                                     <  7:352808/149‑1 (MQ=255)                 cATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATtca                                                                                                                                    <  4:84272/149‑1 (MQ=255)                   ttCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATtcatc                                                                                                                                  <  8:47666/149‑1 (MQ=255)                       aTCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATtcatca                                                                                                                                 >  4:261493/1‑146 (MQ=255)                               gATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATcac                                                                                                                                <  7:45661/139‑1 (MQ=255)                               gATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATcac                                                                                                                                >  8:45661/1‑139 (MQ=255)                               gATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAAt                                                                                                                     >  2:13533/1‑149 (MQ=255)                                 tGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTa                                                                                                                   <  1:13533/149‑1 (MQ=255)                                    cTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGcc                                                                                                                <  2:91827/149‑1 (MQ=255)                                    cTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGcc                                                                                                                >  7:77059/1‑149 (MQ=255)                                      gtgCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTATCGCGATTCATCACAGAA‑TCGAATTAGcc                                                                                                                >  6:348390/1‑147 (MQ=255)                                       tgCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACTGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGa                                                                                                              <  7:355296/148‑1 (MQ=255)                                       tgCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATtcatc                                                                                                                                  <  1:197185/129‑1 (MQ=255)                                       tgCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATtcatc                                                                                                                                  >  2:197185/1‑129 (MQ=255)                                       tgCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGa                                                                                                              <  1:322932/148‑1 (MQ=255)                                          aCTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGATCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGagcag                                                                                                          <  7:340145/149‑1 (MQ=255)                                            tCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGagcagc                                                                                                         >  1:71505/1‑149 (MQ=255)                                              cTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCtt                                                                                                       >  3:329788/1‑148 (MQ=255)                                              cTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTc                                                                                                      <  8:80627/149‑1 (MQ=255)                                              cTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTc                                                                                                      <  5:318709/149‑1 (MQ=255)                                                 gCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAg                                                                                                    <  4:173189/148‑1 (MQ=255)                                                 gCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGt                                                                                                   <  5:235167/149‑1 (MQ=255)                                                  cttcttGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTa                                                                                                  >  3:17247/1‑149 (MQ=255)                                                     cttGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATTCGAATTAGCCGAGCAGCTTCAGTAc                                                                                                 <  1:3609/148‑1 (MQ=255)                                                       tGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGagc                                                                                                            >  8:200183/1‑134 (MQ=255)                                                       tGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGagc                                                                                                            <  7:200183/134‑1 (MQ=255)                                                         aGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGa                                                                                           <  4:329788/149‑1 (MQ=255)                                                             gcgggcgCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATTCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATgg                                                                                        <  5:379204/148‑1 (MQ=255)                                                             gcgggcgCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATgg                                                                                        <  3:155996/148‑1 (MQ=255)                                                             gcgggcgCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATgg                                                                                        <  6:300086/148‑1 (MQ=255)                                                                   gcAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTgg                                                                                 <  7:181787/149‑1 (MQ=255)                                                                     agagagAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGtt                                                                               >  4:235369/1‑149 (MQ=255)                                                                           aaCAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTg                                                                         >  3:264499/1‑149 (MQ=255)                                                                             cAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGc                                                                        >  7:218011/1‑148 (MQ=255)                                                                                 gTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAg                                                                   >  3:184218/1‑149 (MQ=255)                                                                                            acacTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGt                                                          <  4:17247/148‑1 (MQ=255)                                                                                                tCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGc                                                                                       <  1:361594/115‑1 (MQ=255)                                                                                                tCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGc                                                                                       >  2:361594/1‑115 (MQ=255)                                                                                                tCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAg                                                     >  6:17897/1‑149 (MQ=255)                                                                                                tCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAg                                                     >  3:226899/1‑149 (MQ=255)                                                                                                tCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAg                                                     >  4:326575/1‑149 (MQ=255)                                                                                                  atatCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGc                                                    <  3:235369/148‑1 (MQ=255)                                                                                                  atatCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCt                                                   >  7:35537/1‑149 (MQ=255)                                                                                                     tCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTtgc                                                <  3:261493/149‑1 (MQ=255)                                                                                                       cATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑t                                                                                                                          <  1:331101/73‑1 (MQ=255)                                                                                                       cATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑t                                                                                                                          >  2:331101/1‑73 (MQ=255)                                                                                                       cATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGa                                                                                                                       <  8:301233/76‑1 (MQ=255)                                                                                                         tgatgCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCa                                            <  4:226899/149‑1 (MQ=255)                                                                                                                  ccTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgct                                   <  8:77059/149‑1 (MQ=255)                                                                                                                      gTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCtgctt                               >  2:134114/1‑149 (MQ=255)                                                                                                                        aGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCtgcttgc                             >  7:168255/1‑149 (MQ=255)                                                                                                                         gAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCtgcttgct                            <  8:218011/149‑1 (MQ=255)                                                                                                                           aTCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCa                          <  2:71505/149‑1 (MQ=255)                                                                                                                              aGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTAGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgct                                   >  7:60145/1‑137 (MQ=255)                                                                                                                              aGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTAGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgct                                   <  8:60145/137‑1 (MQ=255)                                                                                                                                    cGAAAGGGTGCGGCAACGTACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAg                 >  1:371750/1‑149 (MQ=255)                                                                                                                                       aaGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGttt              >  5:67762/1‑149 (MQ=255)                                                                                                                                         gggTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCg            >  4:353334/1‑149 (MQ=255)                                                                                                                                          ggTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTc                           >  8:104093/1‑133 (MQ=255)                                                                                                                                          ggTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTc                           <  7:104093/133‑1 (MQ=255)                                                                                                                                            tACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGgca         <  8:35537/149‑1 (MQ=255)                                                                                                                                            tACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGgca         <  3:353334/149‑1 (MQ=255)                                                                                                                                                  gCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCGGGTTCGCCTGCGACAGGAGAGCGGTGGACGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCCGGTTCGGCAGCAAAg   >  8:352262/1‑149 (MQ=255)                                                                                                                                                   ccGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGGCAGCAAAGt  >  2:325279/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                      ACGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGGCAGCAAAGT  >  NC_002947/4968001‑4968294

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