Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,950,972	T→G	15.8%	intergenic (‑122/+35)	fliL ← / ← fliK	flagellar protein FliL/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,950,972	0	T	G	15.8%	122.3 / 16.5	62	intergenic (‑122/+35)	fliL/fliK	flagellar protein FliL/flagellar hook‑length control protein FliK
Reads supporting (aligned to +/- strand):  ref base T (28/24);  new base G (5/5);  total (33/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 = 6.59e-01

TCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACC  >  NC_002947/4950825‑4951115                                                                                                                                                    |                                                                                                                                                tCTTTGACTGCTTCGCACTTAGCCATGCCAATAATACGTCGTCAATCGCGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATATTGCAGAGGATGGATTCCAGCAGGTTATGCGGTACTGGCTTCTTCGCGGGCTc                                                                                                                                                >  6:256461/1‑149 (MQ=255)     tGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGccc                                                                                                                                            >  4:38145/1‑149 (MQ=255)     tGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGc                                                                                                                                              <  8:13073/147‑1 (MQ=255)       aCTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCcgc                                                                                                                                          >  1:179990/1‑149 (MQ=255)       aCTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCcgc                                                                                                                                          <  3:38145/149‑1 (MQ=255)         tGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCcgcg                                                                                                                                         >  6:320195/1‑148 (MQ=255)           cTTCGCTCTTCCGCATGCCAATCATCCGGCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaa                                                                                                                                       <  8:265800/148‑1 (MQ=255)              cGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaagaa                                                                                                                                    <  3:138912/148‑1 (MQ=255)              cGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATCGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaagaa                                                                                                                                    <  1:262942/148‑1 (MQ=255)              cGCTCTTCGCCATGCCAATAATCCGGCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaagaa                                                                                                                                    <  2:118818/148‑1 (MQ=255)                 tctTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagc                                                                                                                                  >  6:149431/1‑147 (MQ=255)                  ctTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagccag                                                                                                                               >  2:174530/1‑149 (MQ=255)                          tGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGc                                                                                                                        >  4:24984/1‑148 (MQ=255)                             cAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGccc                                                                                                                                            <  1:17011/125‑1 (MQ=255)                             cAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGccc                                                                                                                                            >  2:17011/1‑125 (MQ=255)                             cAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCaa                                                                                                                     >  1:121651/1‑148 (MQ=255)                              aataatCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCaa                                                                                                                     >  6:11893/1‑147 (MQ=255)                              aataatCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCaa                                                                                                                     >  5:110098/1‑147 (MQ=255)                              aataatCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaa                                                                                                                                       >  3:7037/1‑129 (MQ=255)                              aataatCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaa                                                                                                                                       <  4:7037/129‑1 (MQ=255)                                  atCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTCCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGAGGGCTAGCCCGCGAACAAGCCAGCCCACGCCAACACTTa                                                                                                               >  7:161597/1‑149 (MQ=255)                                         gtcCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCCGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGc                                                                                                         <  1:250100/148‑1 (MQ=255)                                           ccATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgta                                                                                                      <  3:24984/149‑1 (MQ=255)                                            cATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgta                                                                                                      <  2:179990/148‑1 (MQ=255)                                            cATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgta                                                                                                      <  7:24490/148‑1 (MQ=255)                                                 ggggTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTAGCGGGGGGGCCCCCGGAAAAACCAGCCCCGGCCCACAATTATAAGGGGGAGTTAt                                                                                                 >  2:136572/1‑148 (MQ=255)                                                        aaaGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAg                                                                                                           <  2:238480/131‑1 (MQ=255)                                                        aaaGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCACGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAg                                                                                                           >  1:238480/1‑131 (MQ=255)                                                           ggAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCa                                                                                         >  6:15251/1‑146 (MQ=255)                                                           ggAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCa                                                                                         <  5:15251/146‑1 (MQ=255)                                                              gaATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAATCAAGCCCAGGCCAACAGTTATAAGGCGTAGTAATCGACCGCACTGTCg                                                                                   >  4:95613/1‑149 (MQ=255)                                                              gaATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCg                                                                                   <  1:174530/149‑1 (MQ=255)                                                               aaTGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGc                                                                                  <  5:11893/149‑1 (MQ=255)                                                                aTGGCTTGGCTTTACAATGTTTTATGCCTTATTTTGCAGGGTTTGCATACCACCAGTTTTGGCGGGGCGGGCTTCTTCGCGGCCGCCCCCGCTAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGcc                                                                                 <  1:311931/149‑1 (MQ=255)                                                                   gcttggcttgAGAAAGTGTTATGCCATATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGa                                                                               <  1:136572/148‑1 (MQ=255)                                                                      tggattgAGCAAGTGTTATGCCATATTTTGGAGGTGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGCGCAATCGAAGAAGCCAGCCCTGGCCAACAGTTATCAGGCGTAGTTATCGACCGCACTGTCGCCGACaa                                                                            <  5:256461/148‑1 (MQ=255)                                                                        gcttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAacca                                                                         >  1:35718/1‑149 (MQ=255)                                                                             aGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTCTCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACctgc                                                                    <  5:320195/149‑1 (MQ=255)                                                                                       aTGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCTCCGACAACCACCTGCTGCTCCaaa                                                           <  3:95613/148‑1 (MQ=255)                                                                                          ccAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGcc                                                       >  4:22862/1‑149 (MQ=255)                                                                                           cAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGcc                                                       <  6:110098/148‑1 (MQ=255)                                                                                            aGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTAGCCGACAACCACCTGCAGCTCCAAAGGCCgc                                                     >  6:61348/1‑148 (MQ=255)                                                                                                    cAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATc                                             >  4:48751/1‑149 (MQ=255)                                                                                                     aGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCg                                                      >  1:117753/1‑139 (MQ=255)                                                                                                     aGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCg                                                      <  2:117753/139‑1 (MQ=255)                                                                                                       gggATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGAc                                          <  3:22862/149‑1 (MQ=255)                                                                                                          atggatTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCtgctgc                                                                 <  4:96811/123‑1 (MQ=255)                                                                                                          atggatTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCtgctgc                                                                 >  3:96811/1‑123 (MQ=255)                                                                                                          atggatTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGc                                       >  5:7987/1‑149 (MQ=255)                                                                                                             gatTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAG‑TATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCg                                   >  6:240464/1‑149 (MQ=255)                                                                                                                 ccAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAG‑TATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGcc                                <  5:240464/148‑1 (MQ=255)                                                                                                                    gcagGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTTATGGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCa                             >  5:294227/1‑149 (MQ=255)                                                                                                                            tGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCt                      <  6:7987/148‑1 (MQ=255)                                                                                                                                gggCTGGCTTCTTCGCGGGCGAGCCCGCGAAGAAGCCAGCCCCGGCCAACAGTGATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGcc                 >  5:97221/3‑149 (MQ=255)                                                                                                                                gTGCTGGCTTCTTAGCGGGCTCGCCCCCGAAGAAGCCAGCCCAGGCCCACAGTTATAAGGCGTAGTAATCGACCCCACAGTAGCCGACAAACACCAGCTGCTGCAAAGGCCGGGCGCTATAGACCCCATAGGCCCCCCCCAGCTAGGcc                 >  4:153945/1‑149 (MQ=255)                                                                                                                                      gCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGgcgc           >  5:138988/1‑149 (MQ=255)                                                                                                                                        ttcttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCgg                                  >  6:100975/1‑124 (MQ=255)                                                                                                                                        ttcttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCgg                                  <  5:100975/124‑1 (MQ=255)                                                                                                                                        ttcttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTGTCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGgcgcgc         >  1:71547/1‑149 (MQ=255)                                                                                                                                           ttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGgcg                                                  <  2:80885/105‑1 (MQ=255)                                                                                                                                           ttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGgcg                                                  >  1:80885/1‑105 (MQ=255)                                                                                                                                               cgGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCAcc  >  5:301987/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                TCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACC  >  NC_002947/4950825‑4951115

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