Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,968,140	A→G	9.7%	intergenic (‑36/+39)	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,140	0	A	G	9.7%	197.7 / 9.3	72	intergenic (‑36/+39)	PP_4377/fliC	flagellin FlaG/flagellin
Reads supporting (aligned to +/- strand):  ref base A (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 = 1.00e+00

TTGGTCGGCACGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGGCA  >  NC_002947/4967992‑4968287                                                                                                                                                      |                                                                                                                                                     ttGGTCGGCACGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTa                                                                                                                                                      >  4:155996/1‑149 (MQ=255) ttGGTCGGCACGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTa                                                                                                                                                      >  5:300086/1‑149 (MQ=255)   ggTCGGCACGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTa                                                                                                                                                      <  3:370226/147‑1 (MQ=255)          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                                                                                                                           >  2:197185/1‑129 (MQ=255)                                                tgCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATtcatc                                                                                                                           <  1:197185/129‑1 (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                                                                                               <  5:318709/149‑1 (MQ=255)                                                       cTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTc                                                                                               <  8:80627/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                                                                                                     <  7:200183/134‑1 (MQ=255)                                                                tGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGagc                                                                                                     >  8:200183/1‑134 (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                                                                                                                   >  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                    <  7:104093/133‑1 (MQ=255)                                                                                                                                                   ggTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTc                    >  8:104093/1‑133 (MQ=255)                                                                                                                                                     tACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGgca  <  8:35537/149‑1 (MQ=255)                                                                                                                                                     tACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGgca  <  3:353334/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     TTGGTCGGCACGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGGCA  >  NC_002947/4967992‑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: ‑