Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,968,143	T→G	9.7%	intergenic (‑39/+36)	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,143	0	T	G	9.7%	202.0 / 13.0	72	intergenic (‑39/+36)	PP_4377/fliC	flagellin FlaG/flagellin
Reads supporting (aligned to +/- strand):  ref base T (28/37);  new base G (2/5);  total (30/42)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

ACGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGGCA  >  NC_002947/4968001‑4968287                                                                                                                                                |                                                                                                                                                  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                                                                                 <  6:300086/148‑1 (MQ=255)                                                             gcgggcgCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATgg                                                                                 <  3:155996/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                                                                                                                   >  2:331101/1‑73 (MQ=255)                                                                                                       cATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑t                                                                                                                   <  1:331101/73‑1 (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)                                                                                                                                                |                                                                                                                                                  ACGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGGCA  >  NC_002947/4968001‑4968287

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