Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,968,149	C→T	9.3%	intergenic (‑45/+30)	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,149	0	C	T	9.3%	209.6 / 9.7	75	intergenic (‑45/+30)	PP_4377/fliC	flagellin FlaG/flagellin
Reads supporting (aligned to +/- strand):  ref base C (30/38);  new base T (2/5);  total (32/43)
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.01e-01

TCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGGCAGCAAAGTCGG  >  NC_002947/4968004‑4968297                                                                                                                                                   |                                                                                                                                                      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                                                                                                                                   >  8:45661/1‑139 (MQ=255)                            gATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATcac                                                                                                                                   <  7:45661/139‑1 (MQ=255)                            gATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAAt                                                                                                                        >  2:13533/1‑149 (MQ=255)                              tGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTa                                                                                                                      <  1:13533/149‑1 (MQ=255)                                 cTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGcc                                                                                                                   >  7:77059/1‑149 (MQ=255)                                 cTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGcc                                                                                                                   <  2:91827/149‑1 (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                                                                                                         <  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                                                                                                               >  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                                                        >  3:226899/1‑149 (MQ=255)                                                                                             tCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAg                                                        >  4:326575/1‑149 (MQ=255)                                                                                             tCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAg                                                        >  6:17897/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            <  3:353334/149‑1 (MQ=255)                                                                                                                                         tACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGgca            <  8:35537/149‑1 (MQ=255)                                                                                                                                               gCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCGGGTTCGCCTGCGACAGGAGAGCGGTGGACGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCCGGTTCGGCAGCAAAg      >  8:352262/1‑149 (MQ=255)                                                                                                                                                ccGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGGCAGCAAAGt     >  2:325279/1‑149 (MQ=255)                                                                                                                                                 cGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGGCAGCAAAGTc    >  4:369067/1‑149 (MQ=255)                                                                                                                                                  gCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGGCAGCAAAGTCg   <  3:326575/149‑1 (MQ=255)                                                                                                                                                   cACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGGCAGCAAAGTCgg  >  2:65607/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                      TCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCA‑GCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAA‑TCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGGCAGCAAAGTCGG  >  NC_002947/4968004‑4968297

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