Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,951,279	T→C	6.7%	I332V (ATC→GTC)	fliK ←	flagellar hook‑length control protein FliK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,951,279	0	T	C	6.7%	173.5 / 2.8	61	I332V (ATC→GTC)	fliK	flagellar hook‑length control protein FliK
Reads supporting (aligned to +/- strand):  ref base T (30/26);  new base C (2/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

CTGCGCCTGCCCCTGCTGTTGTTGCTGCCCGCGCGACTGGTCAGCCACGTTGACATCCGGCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTCTGGCTGGACAGGTACATGACCCGGTTGACCAGGCCCTCGGCCCAAGCG  >  NC_002947/4951133‑4951427                                                                                                                                                   |                                                                                                                                                      cTGCGCCTGCCCCTGCTGTTGTTGCTGCCCGCGCGACTGGTCAGCCACGTTGACATCCGGCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCtgatg                                                                                                                                                      >  5:329520/1‑148 (MQ=255)     gcCTGGCCCTGCTGGTGTTGCTGCACGCGCGACTGGTCAGCCACGTTGACATCCGGCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGt                                                                                                                                                 <  4:289250/149‑1 (MQ=255)       cTGCCCCTGCTGTTGTTGCTGCCCGCGCGACTGGTCAGCCACGTTGACATCCGGCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGa                                                                                                                                               >  2:158681/1‑149 (MQ=255)       cTGCCCCTGCTGTTGTTGCTGCCCGCGCGACTGGTCAGCCACGTTGACATCCGGCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGACGAAGGTGa                                                                                                                                               >  7:162437/1‑149 (MQ=255)       cTGCCCCTGCTGTTGTTGCTGCCCGCGCGACTGGTCAGCCACGTTGACATCCGGCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGACGAAGGTGa                                                                                                                                               >  3:198541/1‑149 (MQ=255)          ccccTGCTGTTGTTGCTGCCCGCGCGACTGGTCAGCCACGTTGACATCCGGCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCt                                                                                                                                            <  2:180679/149‑1 (MQ=255)            ccTGCTGTTGTTGCTGCCCGCGCGACTGGTCAGCCACGTTGACATCCGGCTGCGCCAACCCCTGCTCGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTgg                                                                                                                                          >  7:274995/1‑149 (MQ=255)                ctgTTGTTGCTGCCCGCGCGACTGGTCAGCCACGTTGACATCCGGCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTg                                                                                                                                       <  6:92692/148‑1 (MQ=255)                  gttgttGCTGCCCGCGCGACTGGTCAGCCACGTTGACATCCGGCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGcc                                                                                                                                    >  7:302916/1‑149 (MQ=255)                        gctgcCCGCGCGACTGGTCAGCCACGTTGACATCCGGCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTc                                                                                                                               <  6:329520/148‑1 (MQ=255)                           gcCCGCGCGACTGGTCAGCCACGTTGACATCCGGCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTg                                                                                                                                  >  6:106986/1‑142 (MQ=255)                           gcCCGCGCGACTGGTCAGCCACGTTGACATCCGGCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTg                                                                                                                                  <  5:106986/142‑1 (MQ=255)                                      tGGTCAGCCACGTTGACATCCGGCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGAc                                                                                                                >  4:32088/1‑149 (MQ=255)                                             ccACGTTGACATCCGGCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATg                                                                                                         >  1:356384/1‑149 (MQ=255)                                                    gACATCCGGCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAgg                                                                                                   <  2:356384/148‑1 (MQ=255)                                                       aTCCGGCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGAc                                                                                                                                              <  4:77090/102‑1 (MQ=255)                                                       aTCCGGCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGAc                                                                                                                                              >  3:77090/1‑102 (MQ=255)                                                         ccGGCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑cc                                                                                             >  2:91778/1‑148 (MQ=255)                                                            gCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCaa                                                                                                                       <  5:315436/120‑1 (MQ=255)                                                            gCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCaa                                                                                                                       >  6:315436/1‑120 (MQ=255)                                                             cTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGtt                                                                                        <  1:158681/149‑1 (MQ=255)                                                                 gcCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTctgc                                                                                    >  3:235642/1‑149 (MQ=255)                                                                      ccccTGCTGGGCAAACAGTTCGCGCGGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCt                                                                               >  4:360302/1‑149 (MQ=255)                                                                        ccTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTc                                                                              >  3:198273/1‑148 (MQ=255)                                                                          tgctgGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAg                                                                           >  7:319636/1‑149 (MQ=255)                                                                               ggCAAACAGTTCGCGCAGGCGGTGCACCTTGCTGTCGAGGGCGTCACGTACGTCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGAcgc                                                                                                              <  1:107598/110‑1 (MQ=255)                                                                               ggCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGa                                                                      <  3:363197/149‑1 (MQ=255)                                                                               ggCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGa                                                                      >  7:292127/1‑149 (MQ=255)                                                                               ggCAAACAGTTCGCGCAGGCGGTGCACCTCGCTGTCGAGGGCGTCACGTACGTCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGAcgc                                                                                                              >  2:107598/1‑110 (MQ=255)                                                                                 cAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGt                                                                                                                                >  1:234426/1‑90 (MQ=255)                                                                                 cAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGt                                                                                                                                <  2:234426/90‑1 (MQ=255)                                                                                  aaaCAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACAt                                                                                                                    >  1:53522/1‑101 (MQ=255)                                                                                  aaaCAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACAt                                                                                                                    <  2:53522/101‑1 (MQ=255)                                                                                  aaaCAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTAACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGatat                                                                   >  5:294282/1‑149 (MQ=255)                                                                                   aaCAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGatat                                                                   >  3:62451/1‑148 (MQ=255)                                                                                   aaCAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATc                                                                  <  2:330307/149‑1 (MQ=255)                                                                                       gTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGACGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGcc                                                              <  8:162437/149‑1 (MQ=255)                                                                                        ttCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGACGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCGGTTGGATATCAGcc                                                              <  4:198541/148‑1 (MQ=255)                                                                                           gcgcAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACt                                                          >  5:330543/1‑149 (MQ=255)                                                                                              cAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGa                                                       <  8:302916/149‑1 (MQ=255)                                                                                                   ggTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTc                                                  <  3:32088/149‑1 (MQ=255)                                                                                                       cACGTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGGCGCTGAGGAAGGTGACCTGGGTGGCCTGTTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCCGTTGGATATCAGCCGACTTGAGGTTctggc                                              <  8:274995/149‑1 (MQ=255)                                                                                                         ccTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTctggctg                                            >  7:36693/1‑149 (MQ=255)                                                                                                          cTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGACCCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTctggctg                                            <  1:91778/149‑1 (MQ=255)                                                                                                          cTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTctggctg                                            >  8:251278/1‑148 (MQ=255)                                                                                                             gCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTCTGGCTGGa                                          <  8:319636/147‑1 (MQ=255)                                                                                                             gCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTCTGGCTGGAc                                         <  4:198273/148‑1 (MQ=255)                                                                                                              cTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTCTGGCTGGACAg                                       <  6:330543/149‑1 (MQ=255)                                                                                                                gTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTCTGGCTGGACAgg                                      >  3:143976/1‑148 (MQ=255)                                                                                                                 tCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTCTGGCTGGACAGGTa                                    <  7:251278/149‑1 (MQ=255)                                                                                                                     gggCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTCTGGCTGGACAGGTACAt                                 >  7:248277/1‑148 (MQ=255)                                                                                                                     gggCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTCTGGCTGGACAGGTACAt                                 >  3:151181/1‑148 (MQ=255)                                                                                                                        cGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTAACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTCTGGCTGGACAGGTACATGAcc                             <  6:294282/149‑1 (MQ=255)                                                                                                                                aCGCCGGCGTGGCCGCTGATGAAGCTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATCTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGTTTCTGGCTGGACAGGTACATGACCCGGTTGAc                     >  5:60673/1‑149 (MQ=255)                                                                                                                                  gCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTCTGGCTGGACAGGTACATGACCCGGTTGAc                     >  6:198690/1‑147 (MQ=255)                                                                                                                                  gCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTCTGGCTGGACAGGTACATGACCCGGTTGAc                     >  6:125487/1‑147 (MQ=255)                                                                                                                                    cGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTCTGGCTGGACAGGTACATGAcc                             <  1:109545/137‑1 (MQ=255)                                                                                                                                    cGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTCTGGCTGGACAGGTACATGAcc                             >  2:109545/1‑137 (MQ=255)                                                                                                                                           gCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTCTGGCTGGACAGGTACATGACCCGGTTGACCAGGCCCTCgg          <  4:151181/149‑1 (MQ=255)                                                                                                                                            ccGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTCTGGCTGGACAGGTACATGACCCGGTTGACCAGGCCCTCGGc         <  5:125487/149‑1 (MQ=255)                                                                                                                                              gCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTCTGGCTGGACAGGTACATGACCCgg                          <  8:110988/130‑1 (MQ=255)                                                                                                                                              gCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTCTGGCTGGACAGGTACATGACCCgg                          >  7:110988/1‑130 (MQ=255)                                                                                                                                                 gatgaAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAg                                                      >  5:121005/1‑99 (MQ=255)                                                                                                                                                 gatgaAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAg                                                      <  6:121005/99‑1 (MQ=255)                                                                                                                                                   tgaAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTCTGGCTGGACAGGTACATGACCCGGTTGACCAGGCCCTCGGCCCAAGCg  <  4:62451/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      CTGCGCCTGCCCCTGCTGTTGTTGCTGCCCGCGCGACTGGTCAGCCACGTTGACATCCGGCTGCGCCAACCCCTGCTGGGCAAACAGTTCGCGCAGGCGGTGCACCTGGCTGTCGAGGGCGTCACGTACGCCGGCGTGGCCGCTGATGAAGGTGACCTGGGTGGCCTGGTCCGCCGCAACATTGACGCGAATGTCGAGGCGA‑CCGAGTTCTGCTGGCTCCAGTTGGATATCAGCCGACTTGAGGTTCTGGCTGGACAGGTACATGACCCGGTTGACCAGGCCCTCGGCCCAAGCG  >  NC_002947/4951133‑4951427

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