Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,950,972	T→G	7.6%	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	7.6%	177.0 / 4.2	64	intergenic (‑122/+35)	fliL/fliK	flagellar protein FliL/flagellar hook‑length control protein FliK
Reads supporting (aligned to +/- strand):  ref base T (27/32);  new base G (3/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCAC  >  NC_002947/4950843‑4951114                                                                                                                                  |                                                                                                                                                ttCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagccagc                                                                                                                              >  8:275128/1‑149 (MQ=255) ttCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagcca                                                                                                                                <  6:254643/147‑1 (MQ=255) ttCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTAGCCCGCGAAGAagccagc                                                                                                                              >  7:360293/1‑149 (MQ=255)         gCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGcc                                                                                                                       >  2:35536/1‑148 (MQ=255)         gCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGcc                                                                                                                       >  7:19607/1‑148 (MQ=255)           cAATAAGCCGTCGTCCATCGGGGGTTCAAAGGCGAATGGCTTGGCTTGAGCAAGTGTTTTGCCAGATTTTGCAGGGCATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAc                                                                                                                    <  7:352211/149‑1 (MQ=255)                  ccGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTAt                                                                                                             >  3:138017/1‑148 (MQ=255)                  ccGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTAGCCCGCGAAGAAGCCAGACCAGGCCAAC‑AGTTAt                                                                                                             >  1:208298/1‑148 (MQ=255)                  ccGGCGCCCCTCGGGGGTCCAAAGCAGAATGGCTTGGCTGGAGCAAGTGTTATCCCAGATTTTGCCGGGGATGGAGTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATc                                                                                                            <  4:194584/149‑1 (MQ=255)                    gtcgtcCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCGGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAg                                                                                                          <  5:92639/149‑1 (MQ=255)                         gcATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGTTCAAGTGTTATGCCAGATTTTGCAGGGGAGTGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCgta                                                                                                     <  7:370278/148‑1 (MQ=255)                          gatggggggTTCATAGTATATTGGCTTGGGTTGAGAATGTGTTATGCCAGATTTTGCAGGGTTTGGTTTCCAGATGGTTTTGGGGTGCTGGCTTCTTCGCGGGCTCGCCCGCTAATAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCgtag                                                                                                    <  2:208298/145‑1 (MQ=255)                           aTCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAg                                                                                                          >  6:203310/1‑142 (MQ=255)                           aTCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAg                                                                                                          <  5:203310/142‑1 (MQ=255)                           aTCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGACCAGGCCAAC‑AGTTATCATGCgtag                                                                                                    >  5:354888/1‑148 (MQ=255)                               ggggTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATc                                                                                               <  4:136629/149‑1 (MQ=255)                                    tCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGc                                                                                                        <  3:300092/135‑1 (MQ=255)                                    tCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGc                                                                                                        >  4:300092/1‑135 (MQ=255)                                      aaaGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGAcc                                                                                           <  7:103969/146‑1 (MQ=255)                                      aaaGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGAcc                                                                                           >  8:103969/1‑146 (MQ=255)                                      aaaGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCa                                                                                        >  7:263536/1‑149 (MQ=255)                                      aaaGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCa                                                                                        >  3:66501/1‑149 (MQ=255)                                          gagaATGGCTTGGCTTGAGCAAGTGATATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAg                                                                                                          >  7:69152/1‑127 (MQ=255)                                          gagaATGGCTTGGCTTGAGCAAGTGATATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAg                                                                                                          <  8:69152/127‑1 (MQ=255)                                            gaATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTc                                                                                   <  4:131001/148‑1 (MQ=255)                                            gaATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTc                                                                                   >  3:131001/1‑148 (MQ=255)                                              aTGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGcc                                                                                <  7:51357/149‑1 (MQ=255)                                                 gcttgggttgAGCATGTGTTATGCCATATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCgtag                                                                                                    <  3:408252/126‑1 (MQ=255)                                                 gcttggcttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGa                                                                              <  8:19607/148‑1 (MQ=255)                                                 gcttggcttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAACCCAGGCCAAC‑AGTTATCAGGCGTAGTAAAAGACAGCAATGTCGCCGa                                                                              >  8:58507/1‑148 (MQ=255)                                                      gcttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAacc                                                                         <  8:263536/148‑1 (MQ=255)                                                        ttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAg                                                                                                          >  3:75535/1‑113 (MQ=255)                                                        ttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAg                                                                                                          <  4:75535/113‑1 (MQ=255)                                                              aaGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCtgctgc                                                                >  5:96480/1‑149 (MQ=255)                                                               aGTGTTATGCCAGATTTTGAAGGGTATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCTAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTgctgct                                                               <  1:160832/149‑1 (MQ=255)                                                                gtgtTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTc                                                              >  8:326606/1‑149 (MQ=255)                                                                     aTGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCaaa                                                          <  3:167227/148‑1 (MQ=255)                                                                     aTGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCaaa                                                          <  4:138017/148‑1 (MQ=255)                                                                         cAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCg                                                     >  7:289959/1‑149 (MQ=255)                                                                          aGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCgg                                                    >  4:350622/1‑149 (MQ=255)                                                                          aGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCgg                                                    <  6:96480/149‑1 (MQ=255)                                                                                 gCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTAt                                             >  5:24989/1‑149 (MQ=255)                                                                                  cAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTAt                                             <  7:326606/148‑1 (MQ=255)                                                                                    ggggATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGa                                          <  7:275128/149‑1 (MQ=255)                                                                                     gggATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACctgctg                                                                 <  8:27388/125‑1 (MQ=255)                                                                                                  gcagGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCa                            <  8:289959/149‑1 (MQ=255)                                                                                                            cGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCgg                  <  6:24989/149‑1 (MQ=255)                                                                                                            cGGTGATGGCTTCTTCTCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCg                   <  8:360293/148‑1 (MQ=255)                                                                                                                gatatattcttcGCTGGCGAGCCCGCGAAGAAGCCATCCCAGGCCAAC‑AGTTATCAGTCGTAGTATTCGACGGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGCCCGCTTCGGCCCCACCCTGCTCGGCCCg              <  5:299288/143‑1 (MQ=255)                                                                                                                 cTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCgg             <  2:56997/149‑1 (MQ=255)                                                                                                                  tGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACGAGTTATCAGGCGTAGTAATCGACCGCACTGTCg                                                                                  >  1:90885/1‑80 (MQ=255)                                                                                                                  tGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCg                                                                                  <  2:90885/79‑1 (MQ=255)                                                                                                                   ggCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTATTCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCcccccc                          >  1:162377/1‑134 (MQ=255)                                                                                                                    gCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCAcc                          >  7:313187/1‑133 (MQ=255)                                                                                                                    gCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCAcc                          <  8:313187/133‑1 (MQ=255)                                                                                                                      ttcttcGCGGGCTCGCCCGCGAAGAAGCCCGCCCAGGCCAAC‑AGATATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCAGCTCCACCAAAGGCCCGGCGCTATCGAACGCTTCGGCCACAACCTGCTCGGCCCGccgcgc        >  7:176508/1‑149 (MQ=255)                                                                                                                      ttcttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGgcgcgc        >  2:348381/1‑149 (MQ=255)                                                                                                                      ttcttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCCGGACAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCACAAGGCCGCGCGCTATCGACAGCTTCGGCCCAACCCTGcacggcacggagaaa        >  6:299288/1‑144 (MQ=255)                                                                                                                      ttcttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGgcgcgc        >  1:273988/1‑149 (MQ=255)                                                                                                                        cttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgcc       <  4:66501/148‑1 (MQ=255)                                                                                                                         ttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgccgc     >  2:260567/1‑149 (MQ=255)                                                                                                                          tcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCGCGCCGggccgag    >  3:318682/1‑147 (MQ=255)                                                                                                                            gcgGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGNCCCACCCTGCTCGGCCCGGCGCGCCGCCAc  <  1:35536/149‑1 (MQ=255)                                                                                                                            ccgGGCGCGCCCGCTAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCAc  <  5:66895/148‑1 (MQ=255)                                                                                                                                  |                                                                                                                                                TTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAC‑AGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCAC  >  NC_002947/4950843‑4951114

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