Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,968,149	C→T	5.8%	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	5.8%	263.7 / 3.3	86	intergenic (‑45/+30)	PP_4377/fliC	flagellin FlaG/flagellin
Reads supporting (aligned to +/- strand):  ref base C (45/36);  new base T (2/3);  total (47/39)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.99e-01

CGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGGCAGCAAAGTCG  >  NC_002947/4968002‑4968296                                                                                                                                                    |                                                                                                                                                    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                                                                                                                     <  3:96444/149‑1 (MQ=255)                                tGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTa                                                                                                                     <  4:78977/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                                                              <  7:15601/140‑1 (MQ=255)                                                                                                atatCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCg                                                              >  8:15601/1‑140 (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                                   <  6:385224/147‑1 (MQ=255)                                                                                                                    gTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgctgc                                   <  5:241526/147‑1 (MQ=255)                                                                                                                    gTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgctgc                                   <  5:247543/147‑1 (MQ=255)                                                                                                                    gTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgctgc                                   <  1:44526/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                                          >  7:117246/1‑136 (MQ=255)                                                                                                                        aaTCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGt                                          <  8:117246/136‑1 (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)                                                                                                                                                  cGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCa                                                                                        <  3:205791/64‑1 (MQ=255)                                                                                                                                                  cGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCa                                                                                        >  4:205791/1‑64 (MQ=255)                                                                                                                                                  cGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGGCAGCAAAGTc   >  7:196294/1‑149 (MQ=255)                                                                                                                                                   gCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGGCAGCAAAGTCg  <  6:408638/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGGCAGCAAAGTCG  >  NC_002947/4968002‑4968296

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