Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,948,667:1	+A	5.6%	coding (181/435 nt)	fliO ←	flagellar protein FliO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,948,667	1	.	A	5.6%	287.2 / 11.0	71	coding (181/435 nt)	fliO	flagellar protein FliO
Reads supporting (aligned to +/- strand):  ref base . (38/29);  new base A (2/2);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

TGCAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCGTGGCTGG  >  NC_002947/4948522‑4948813                                                                                                                                                    |                                                                                                                                                    tGCAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑gcg                                                                                                                                                  >  5:128511/1‑149 (MQ=255)  gCAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑gcgc                                                                                                                                                 >  1:347870/1‑149 (MQ=255)  gCAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑gcgc                                                                                                                                                 <  4:165402/149‑1 (MQ=255)   cAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑gcgc                                                                                                                                                 <  2:347870/148‑1 (MQ=255)            gATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑Gcagcag                                                                                                                                        <  7:210191/147‑1 (MQ=255)            gATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAgcc                                                                                                                                      >  4:93600/1‑149 (MQ=255)             ctGGTGCCGGGGGTATGGCCGATCAGGCTCTGCTCCTTGCCCCCCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCGCTGCGGCGC‑CTGCAGAGCACC‑gcagtgaagg                                                                                                                                     <  5:211844/148‑7 (MQ=255)               ggTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAgccagg                                                                                                                                   >  2:261645/1‑149 (MQ=255)                gTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAgccaggc                                                                                                                                  <  6:10125/149‑1 (MQ=255)                  gCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGTGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAgccaggc                                                                                                                                  >  1:289871/1‑148 (MQ=255)                     gggggTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCagaa                                                                                                                             <  6:273787/149‑1 (MQ=255)                      ggggTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCaga                                                                                                                              >  2:211289/1‑147 (MQ=255)                      ggggTATGGCCGATCAGGATCTGATCCTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCTGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCaga                                                                                                                              >  8:92808/1‑147 (MQ=255)                       gggTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCagaaag                                                                                                                           <  3:313725/149‑1 (MQ=255)                          tATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAaagaag                                                                                                                        <  6:10268/149‑1 (MQ=255)                               ccGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGAGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGAC‑GCAGCAGCCAGGCCAGAAAGAAGATCAg                                                                                                                   >  1:268836/1‑149 (MQ=255)                                cGATCAGGATCTGCTCCTTGCCCAACTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAgg                                                                                                                  <  8:244039/149‑1 (MQ=255)                                  aTCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAgccaggcca                                                                                                                                >  6:106242/1‑133 (MQ=255)                                  aTCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAgccaggcca                                                                                                                                <  5:106242/133‑1 (MQ=255)                                      ggATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCccacc                                                                                                            >  5:22943/1‑149 (MQ=255)                                                cTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACcagcagc                                                                                                   <  1:351133/148‑1 (MQ=255)                                                 ttGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAgcca                                                                                                                                     <  4:22827/113‑1 (MQ=255)                                                 ttGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAgcca                                                                                                                                     >  3:22827/1‑113 (MQ=255)                                                    cccACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGcc                                                                                              <  5:49921/149‑1 (MQ=255)                                                     ccACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGcc                                                                                              >  7:29723/1‑148 (MQ=255)                                                     ccACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGcc                                                                                              >  6:140488/1‑148 (MQ=255)                                                      cACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGccc                                                                                             >  4:257364/1‑148 (MQ=255)                                                        ccTGCACCAGCAGCAGGCG‑GTCACGCGGTCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACcag                                                                                                       <  4:126623/136‑1 (MQ=255)                                                        ccTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAg                                                                                          >  3:259287/1‑149 (MQ=255)                                                        ccTGCACCAGCAGCAGGCG‑GTCACGCGGACCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACcag                                                                                                       >  3:126623/1‑136 (MQ=255)                                                         cTGCACCAGCAGCAGGCG‑GTCACGCGGGCCTATGGCACGGCTGCAGACAATCTCGATCACCTGCCCACCCTTCACTGCTGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGTCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGa                                                                                         <  2:268836/149‑1 (MQ=255)                                                          tGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGAc                                                                                        >  5:368838/1‑149 (MQ=255)                                                           gCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGAc                                                                                        >  3:80356/1‑148 (MQ=255)                                                            cACCAGCAGCAGGCG‑GTCACGCGGGCTGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCa                                                                                      <  5:16151/149‑1 (MQ=255)                                                             aCCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCAt                                                                                     >  8:215018/1‑149 (MQ=255)                                                               cagcagcagGCG‑GTCACGCGGGCCGATGGCACGGCTTCCGACAATCTCGATCACCTGCCCACCATTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCCGCAGCAGGCCCAAGACCATCt                                                                                   >  5:19694/1‑149 (MQ=255)                                                                agcagcagGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTg                                                                                  >  2:295706/1‑149 (MQ=255)                                                                      agGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAg                                                                            >  5:342616/1‑149 (MQ=255)                                                                       ggCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCccac                                                                                                             <  8:87081/115‑1 (MQ=255)                                                                             tCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGccc                                                                      >  6:77108/1‑149 (MQ=255)                                                                                 gcgGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACcagcagc                                                                                                   >  3:298646/1‑116 (MQ=255)                                                                                 gcgGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACcagcagc                                                                                                   <  4:298646/116‑1 (MQ=255)                                                                                             cACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGt                                                      <  4:259287/149‑1 (MQ=255)                                                                                                ggCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGa                                                   <  6:368838/149‑1 (MQ=255)                                                                                                  cTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCACGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAgc                                                 <  2:289871/149‑1 (MQ=255)                                                                                                          aaTCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCgggg                                         >  6:367389/1‑149 (MQ=255)                                                                                                           aTCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGt                                        <  1:295706/149‑1 (MQ=255)                                                                                                              tcGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGgc                                     >  5:275605/1‑149 (MQ=255)                                                                                                              tcGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGgc                                     >  3:363912/1‑149 (MQ=255)                                                                                                                  tCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGTCCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGgctgct                                 <  4:246497/149‑1 (MQ=255)                                                                                                                   cACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGACAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGctgctg                                <  8:121156/149‑1 (MQ=255)                                                                                                                     ccTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGc                              >  2:315992/1‑149 (MQ=255)                                                                                                                       tGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCTGCAGCAGCCAGGCTAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCg                             >  1:239105/1‑149 (MQ=255)                                                                                                                         cccACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTgg                          >  2:372963/1‑149 (MQ=255)                                                                                                                           cACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGcc                        >  1:39997/1‑149 (MQ=255)                                                                                                                           cACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGcc                        >  6:186034/1‑149 (MQ=255)                                                                                                                           cACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGcc                        >  4:128756/1‑149 (MQ=255)                                                                                                                            aCCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCg                       <  3:128756/149‑1 (MQ=255)                                                                                                                               cTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAg                    <  1:52184/149‑1 (MQ=255)                                                                                                                                  cACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCca                 >  7:266119/1‑149 (MQ=255)                                                                                                                                   gcTGCGGCGC‑CTGCAGAGCACC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTg                           >  5:52108/2‑138 (MQ=255)                                                                                                                                   gcTGCGGCGC‑CTGCAGAGCACC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTg                           <  6:52108/137‑1 (MQ=255)                                                                                                                                   aCTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCcag                <  6:22943/149‑1 (MQ=255)                                                                                                                                   aCTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCcag                >  6:187968/1‑149 (MQ=255)                                                                                                                                   aCTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCcag                <  5:339880/149‑1 (MQ=255)                                                                                                                                    cTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCcagg               >  4:202913/1‑149 (MQ=255)                                                                                                                                    cTGCGGCGC‑CTGCAGAGCACC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCcagg               >  6:345352/1‑149 (MQ=255)                                                                                                                                     tGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCcaggc              >  3:225482/1‑149 (MQ=255)                                                                                                                                         gTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCaggcaggc          <  6:342616/149‑1 (MQ=255)                                                                                                                                          tGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCg         >  4:352097/1‑149 (MQ=255)                                                                                                                                            ctctGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCGt        <  4:200957/148‑1 (MQ=255)                                                                                                                                             tctGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCGtgg      >  3:66743/1‑149 (MQ=255)                                                                                                                                               tGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCGtggct    >  1:351112/1‑149 (MQ=255)                                                                                                                                               tGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCGtgg      <  5:335843/147‑1 (MQ=255)                                                                                                                                                 cAG‑GCGCC‑GCAGCAGCCAGGCTAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCGtggctgg  <  2:239105/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TGCAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCG‑GTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCC‑GCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCGTGGCTGG  >  NC_002947/4948522‑4948813

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