Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,968,144	T→C	5.9%	intergenic (‑40/+35)	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,144	0	T	C	5.9%	249.0 / 4.0	85	intergenic (‑40/+35)	PP_4377/fliC	flagellin FlaG/flagellin
Reads supporting (aligned to +/- strand):  ref base T (44/36);  new base C (2/3);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

GGCACGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGG  >  NC_002947/4967998‑4968285                                                                                                                                                   |                                                                                                                                              ggCACGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGAACGTTGc                                                                                                                                             >  6:219350/1‑149 (MQ=255) ggCAAGTCTGGTAGCTCAGCATTCGATCTTCTCGATGTCTTTGCACTCCTGTCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGc                                                                                                                                             <  5:311059/149‑1 (MQ=255)   cACGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCg                                                                                                                                           <  5:50804/149‑1 (MQ=255)     cGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCAGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCa                                                                                                                                         <  1:194988/149‑1 (MQ=255)           gTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCttt                                                                                                                                   <  7:117085/149‑1 (MQ=255)           gTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCtt                                                                                                                                    <  2:388717/148‑1 (MQ=255)            tAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTc                                                                                                                                  >  3:78977/1‑149 (MQ=255)              gCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGGAAGGGTACGTTGCCGCACCCtt                                                                                                                                    >  1:372760/1‑145 (MQ=255)              gCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTGCGGCAACGTACCCTTTcg                                                                                                                                 <  3:367207/148‑1 (MQ=255)              gCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTcgc                                                                                                                                >  5:346414/1‑149 (MQ=255)                tCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGa                                                                                                                              >  7:362222/1‑149 (MQ=255)                  cgCATTCGATCTGCGCGTTGGCTGTGCACTCCTGGCTTCGTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAACACGAAACGCTACGTTGCCGCACCCTTTCGCGa                                                                                                                              <  6:130606/146‑1 (MQ=255)                  aGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGAt                                                                                                                             >  8:11401/1‑148 (MQ=255)                   gCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATtc                                                                                                                           <  7:11401/149‑1 (MQ=255)                     aTTCTATATGCTCGATGGCTGTGCACTCCTGGATTCTTGTGTTGCGTGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCATAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATtcat                                                                                                                         <  2:14593/149‑1 (MQ=255)                               cTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTGCGGCAACGTACCCTTTcgt                                                                                                                                <  7:69727/132‑2 (MQ=255)                               cTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTGCGGCAACGTACCCTTTcgt                                                                                                                                >  8:69727/1‑131 (MQ=255)                               cTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATcaca                                                                                                                     <  7:30914/143‑1 (MQ=255)                               cTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATc                                                                                                                >  4:271522/1‑148 (MQ=255)                                 cGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGaa                                                                                                             >  4:285971/1‑149 (MQ=255)                                    tGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTa                                                                                                          <  4:78977/149‑1 (MQ=255)                                    tGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTa                                                                                                          <  3:96444/149‑1 (MQ=255)                                           gCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGag                                                                                                    >  6:247543/1‑148 (MQ=255)                                           gCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGAGATTCATCACAGAATCGAATTAGCCGag                                                                                                    >  6:241526/1‑148 (MQ=255)                                             aCTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGagcag                                                                                                 >  3:261531/1‑149 (MQ=255)                                             aCTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGagcag                                                                                                 >  6:229101/1‑149 (MQ=255)                                                  tGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCa                                                                                            >  2:44526/1‑149 (MQ=255)                                                         ttGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCg                                                                                      >  5:10645/1‑148 (MQ=255)                                                           gAGTGGCTGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGTGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGaa                                                                                    <  2:311048/148‑1 (MQ=255)                                                           gAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGaa                                                                                    <  3:285971/148‑1 (MQ=255)                                                           gAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGaa                                                                                    <  5:229101/148‑1 (MQ=255)                                                            aGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGa                                                                                  <  4:261531/149‑1 (MQ=255)                                                             gTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTTCCGAAGAt                                                                                 <  8:362222/149‑1 (MQ=255)                                                             gTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGAt                                                                                 <  6:346414/149‑1 (MQ=255)                                                                gcgggcgCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGCAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGc                                                                              <  2:372760/149‑1 (MQ=255)                                                                       cAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGt                                                                       >  7:317698/1‑149 (MQ=255)                                                                          agagaACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATTGCATCTGGTTCg                                                                    >  3:51198/1‑149 (MQ=255)                                                                            agaACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGtt                                                                      >  5:345717/1‑145 (MQ=255)                                                                            agaACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCATCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGtt                                                                      <  6:345717/145‑1 (MQ=255)                                                                              aaCAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTg                                                                >  7:187797/1‑149 (MQ=255)                                                                                   agTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGa                                                             >  8:81618/1‑147 (MQ=255)                                                                                    gTTCAGCTTGACACGCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAg                                                          >  7:367731/1‑149 (MQ=255)                                                                                     ttCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCCAATTAGCAGAGCAGCTTCAGAACCGAAGATCGCAGATTGTTCGCATGCGACac                                                          >  5:292078/1‑147 (MQ=255)                                                                                          cTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGc                                                    <  2:219395/149‑1 (MQ=255)                                                                                            tGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCgg                                                  >  4:135898/1‑149 (MQ=255)                                                                                             gACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGGAGCTGGTTCGCCTGCGACAGGATAGCGGt                                                 >  5:385224/1‑149 (MQ=255)                                                                                             gACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCgg                                                  <  3:135898/148‑1 (MQ=255)                                                                                             gACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCgg                                                  <  3:271522/148‑1 (MQ=255)                                                                                              acacTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGt                                                 <  8:317698/148‑1 (MQ=255)                                                                                                acTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTgg                                               >  2:402389/1‑148 (MQ=255)                                                                                                  tCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAg                                            >  1:167676/1‑149 (MQ=255)                                                                                                  tCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAg                                            >  5:306936/1‑149 (MQ=255)                                                                                                    atatCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCg                                                   >  8:15601/1‑140 (MQ=255)                                                                                                    atatCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCg                                                   <  7:15601/140‑1 (MQ=255)                                                                                                    atatCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCt                                          >  5:156096/1‑149 (MQ=255)                                                                                                        ccATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTtgct                                      >  5:408638/1‑149 (MQ=255)                                                                                                          atgatgCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTtgctg                                     >  2:27170/1‑148 (MQ=255)                                                                                                              tgCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAgg                                 <  6:10645/148‑1 (MQ=255)                                                                                                              tgCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCAGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAgg                                 <  6:274841/148‑1 (MQ=255)                                                                                                               gCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGt                               >  7:351793/1‑149 (MQ=255)                                                                                                                 ttGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGttt                             >  6:17650/1‑149 (MQ=255)                                                                                                                 ttGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTATCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGttt                             >  1:95901/1‑149 (MQ=255)                                                                                                                    ccTCGTAGAATCAGAAAACGAAAGGGGACGTTGCCGCACCCTTTCGCGCTTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGGCTGCGTCAGGAGAGCGGTGGACGCTTGCTGCCGGGTTtgct                          >  2:266786/1‑149 (MQ=255)                                                                                                                      tCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgctgc                        >  6:10375/1‑149 (MQ=255)                                                                                                                        gTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgctgc                        <  1:402389/147‑1 (MQ=255)                                                                                                                        gTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgctgc                        <  1:44526/147‑1 (MQ=255)                                                                                                                        gTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgctgc                        <  5:247543/147‑1 (MQ=255)                                                                                                                        gTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgctgc                        <  5:241526/147‑1 (MQ=255)                                                                                                                        gTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgctgc                        <  6:385224/147‑1 (MQ=255)                                                                                                                        gTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCtgctt                      <  2:357533/149‑1 (MQ=255)                                                                                                                         tAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGaa                                             <  4:306040/125‑1 (MQ=255)                                                                                                                         tAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGaa                                             >  3:306040/1‑125 (MQ=255)                                                                                                                         tAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCtgcttg                     >  5:256617/1‑149 (MQ=255)                                                                                                                           taaTCATAAAACGAAATGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACATGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCtgcttgct                   <  8:367731/148‑1 (MQ=255)                                                                                                                            aaTCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGt                               <  8:117246/136‑1 (MQ=255)                                                                                                                            aaTCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGt                               >  7:117246/1‑136 (MQ=255)                                                                                                                              tCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAg                >  2:10688/1‑149 (MQ=255)                                                                                                                                 gAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTc             >  3:125263/1‑149 (MQ=255)                                                                                                                                 aatcgcGAAAGGGTGCGGCAACGTACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgctg                         <  5:219973/132‑1 (MQ=38)                                                                                                                                 aatcgcGAAAGGGTGCGGCAACGTACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgctg                         >  6:219973/6‑137 (MQ=38)                                                                                                                                    aaCGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGTGTTTGCTGCTTGCTCAGCTcag           >  5:116373/1‑146 (MQ=255)                                                                                                                                      cGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAg        >  4:297872/1‑149 (MQ=255)                                                                                                                                      cGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTACGCAg        >  4:292261/1‑149 (MQ=255)                                                                                                                                          aGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTc    <  6:256617/149‑1 (MQ=255)                                                                                                                                            ggTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCgg  >  1:50378/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                              GGCACGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGG  >  NC_002947/4967998‑4968285

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑