Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	709,197	C→A	10.8%	intergenic (+60/‑59)	murJ → / → ribF	lipid II flippase/bifunctional riboflavin kinase/FMN adenylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	709,197	0	C	A	10.8%	178.2 / 11.4	65	intergenic (+60/‑59)	murJ/ribF	lipid II flippase/bifunctional riboflavin kinase/FMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (30/27);  new base A (4/3);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

ATCCTGGCGGGCGTCGTGACGTATTTTGGCTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTGGCAACTTCGACGGGGTTCACCGCGGC  >  NC_002947/709051‑709345                                                                                                                                                     |                                                                                                                                                     aTCCTGGCGGGCGTCGTGACGTATTTTGGCTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCt                                                                                                                                                    <  3:41486/149‑1 (MQ=255)       gcgggcgTCGTGACGTATTTTGGCTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCAc                                                                                                                                               <  2:182103/148‑1 (MQ=255)            cgtcgtGACGTATTTTGGCTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGGTGACAGGCAGCAgcgc                                                                                                                                         <  4:309677/149‑1 (MQ=255)               cgtGACGTATTTTGGCTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGGTGACAGGCAGCAGCGccc                                                                                                                                       >  1:127338/1‑148 (MQ=255)                gtGACGTATTTTGGCTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGGTGACAGGCAGCAGCGcccaa                                                                                                                                     <  2:127338/149‑3 (MQ=255)                gtGACGTATTTTGGCTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCgg                                                                                                                                     <  4:265237/149‑1 (MQ=255)                     gTATTTTGGCTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGtgtg                                                                                                                                >  2:98113/1‑149 (MQ=255)                     gTATTTTGGCTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTAACCAGCGCCCGGGtgtg                                                                                                                                >  7:161631/1‑149 (MQ=255)                     gTATTTTGGCTGCCTGTACCTGTGTGGCTTCCGCCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGGTTTTGGCGGGTGAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGtgtg                                                                                                                                <  2:190664/149‑1 (MQ=255)                             gCTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAAt                                                                                                                        <  7:329904/149‑1 (MQ=255)                              cTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGTCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATc                                                                                                                       <  5:263912/149‑1 (MQ=255)                                gCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTA‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTta                                                                                                                            >  7:144082/1‑142 (MQ=255)                                gCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTA‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTta                                                                                                                            <  8:144082/142‑1 (MQ=255)                                 ccTGTACCTGTGTGGCTTCCGGCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCgg                                                                                                                     <  1:98113/148‑1 (MQ=255)                                  cTGTACCTGTGTGGCTTCCGGCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGc                                                                                                                    >  2:291080/1‑148 (MQ=255)                                    gTACCTGTGTGGCTTCCGGCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCAc                                                                                                                 >  1:5870/1‑149 (MQ=255)                                         tgtgtgGCTTCCGGCCCCGGCATTTCAGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCAGCCACTTTa                                                                                                             >  1:4478/1‑149 (MQ=255)                                           tgtgGCTTCCGGCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTAt                                                                                                            <  3:31342/147‑1 (MQ=255)                                             tgGCTTCCGGGCCCGGCATTTC‑GCTCGCAAGGCGTTGCACTGAGGCTTTTGGCTGGTCAGGGGTCGGTTTTTCTCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTCTCACCAGCGCCCGGGTGTGGTTATTATCGGGCACTTTATGa                                                                                                          <  6:163659/147‑1 (MQ=255)                                             tgGCTTCCGGCCCCGGCATTTC‑GCTCCCAAGGCCTTGCAGTGAGGCTTTTGCCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGGCTGTCACCAGCGCCCGGGTGTGGTTTTATTCGGGCACTTTATGa                                                                                                          >  5:163659/1‑147 (MQ=255)                                               gCTTCCGGCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGc                                                                                                        >  6:179288/1‑147 (MQ=255)                                                 ttCCGGCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCaag                                                                                                     >  8:75951/1‑148 (MQ=255)                                                   ccGGCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCaagaag                                                                                                  >  2:69164/1‑149 (MQ=255)                                                       ccccGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGcgcg                                                                                              <  6:158263/149‑1 (MQ=255)                                                                 tc‑gctcgcAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCt                                                                                    >  8:17868/1‑149 (MQ=255)                                                                  c‑gctcgcAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTg                                                                                   >  6:353708/1‑149 (MQ=255)                                                                    gctcgcAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTgg                                                                                  >  8:172509/1‑149 (MQ=255)                                                                    gctcgcAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTgg                                                                                  >  8:167574/1‑149 (MQ=255)                                                                             gCCTTGCACTCAGGCTTATGGCGGGTCAGGCGTAGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGACACCAGCGCCCGGGTGTGGTTATAATCGGACACTTTATGAGCATGAAGCGCGAGATGCAGCTGGTTCGAGGTc                                                                         >  6:384751/1‑149 (MQ=255)                                                                              ccTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCt                                                                        <  5:353708/149‑1 (MQ=255)                                                                                  gCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTc                                                                      <  1:291080/147‑1 (MQ=255)                                                                                     cTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCAGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAAc                                                                 <  2:4478/149‑1 (MQ=255)                                                                                        aGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTg                                                              >  6:253351/1‑149 (MQ=255)                                                                                          gCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTgc                                                             >  1:91997/1‑149 (MQ=255)                                                                                           cTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTgcgc                                                           <  2:40549/149‑1 (MQ=255)                                                                                            ttttGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGcc                                                          >  7:200219/1‑149 (MQ=255)                                                                                            ttttGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGcc                                                          >  3:120734/1‑149 (MQ=255)                                                                                                ggCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGcccc                                                        >  8:99521/1‑147 (MQ=255)                                                                                                   gggTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCAGAACCTGCGCCCCGAGCa                                                   >  2:374180/1‑149 (MQ=255)                                                                                                     gTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTgg                                                                                                                               >  5:307963/1‑71 (MQ=255)                                                                                                     gTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTgg                                                                                                                               <  6:307963/71‑1 (MQ=255)                                                                                                      tCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCg                                                <  8:200219/149‑1 (MQ=255)                                                                                                                gTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTc                                      >  4:311491/1‑149 (MQ=255)                                                                                                                   tttCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGc                                    >  6:275403/1‑148 (MQ=255)                                                                                                                        cATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCgggg                                             <  3:29898/134‑1 (MQ=255)                                                                                                                        cATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCgggg                                             >  4:29898/1‑134 (MQ=255)                                                                                                                        cATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCAt                              >  5:173819/1‑149 (MQ=255)                                                                                                                         aTTCCACGCCGTCCTT‑GTGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTTGTTATAATCGGCCACTTTATGAGCAAGAATCGCGTTATTCAGCTGGTTCGATGTCTTCACAACCTGCGCCACGAGCACCGGGGCTGTGTCGCCACCAtt                             <  2:142650/149‑1 (MQ=255)                                                                                                                         aTTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCAtt                             <  5:179288/149‑1 (MQ=255)                                                                                                                         aTTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCACCGCCCGGGTGTGGTTATAATCGGCCACTGTATGAGCAAGAAGCGCGTTATTCCGCTGGTTCGACCTCTACCCAACCTTCGCCCAGAGCACCGGGGCTGTGTCGCCACCAtt                             >  1:142650/1‑149 (MQ=255)                                                                                                                          ttCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTg                            >  5:210602/1‑149 (MQ=255)                                                                                                                            ccACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTgg                           <  4:324718/148‑1 (MQ=255)                                                                                                                            ccACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTgg                           <  4:319445/148‑1 (MQ=255)                                                                                                                                gccgccCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTGTTCACAACCTGCGCCCCGAGCAGCGGGGCTGTGTCGCCACCATTGGCAACt                      >  2:299088/1‑149 (MQ=255)                                                                                                                                gccgccCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTGGCAACt                      >  2:303683/1‑149 (MQ=255)                                                                                                                                  cgccCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTGGCAACTTc                    <  2:91997/149‑1 (MQ=255)                                                                                                                                     ccTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCAGAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTGGCAACTTCGAc                 <  1:374180/149‑1 (MQ=255)                                                                                                                                     ccTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTGGCAACTTCGAc                 <  7:17868/149‑1 (MQ=255)                                                                                                                                      acccgggCGCTGGTGACAGGCAGCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCtt                                                                       <  5:167568/91‑1 (MQ=37)                                                                                                                                      acccgggCGCTGGTGACAGGCAGCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCtt                                                                       >  6:167568/5‑95 (MQ=37)                                                                                                                                       tt‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTGGCAACTTCGACgg               >  3:298266/1‑149 (MQ=255)                                                                                                                                         cgggCGCTGGTGACAGGCAGCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTGGCAACTTCGACggg              >  6:376016/2‑149 (MQ=255)                                                                                                                                          gggCGCTGGTGACAGGCAGCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTGGCAACTTCGACgggg             >  8:71372/1‑149 (MQ=255)                                                                                                                                          gggCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTGGCAACTTCGACgg               <  4:120734/147‑1 (MQ=255)                                                                                                                                             cgcTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATACTCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTGGCAACTTCGACGGGGTTc          <  8:161631/149‑1 (MQ=255)                                                                                                                                                     ccTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTGGCAACTTCGACGGGGTTCACCGCGGc  >  7:20640/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     ATCCTGGCGGGCGTCGTGACGTATTTTGGCTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTC‑GCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTT‑GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTGGCAACTTCGACGGGGTTCACCGCGGC  >  NC_002947/709051‑709345

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