Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,964,646	T→C	6.8%	N238D (AAC→GAC)	fleQ ←	transcriptional regulator FleQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,964,646	0	T	C	6.8%	166.7 / 2.5	59	N238D (AAC→GAC)	fleQ	transcriptional regulator FleQ
Reads supporting (aligned to +/- strand):  ref base T (25/30);  new base C (2/2);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

CGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGG  >  NC_002947/4964500‑4964791                                                                                                                                                   |                                                                                                                                                      cttCGCAGCGATGATGAGCACGTCGATGCGTGGGGTCTTGGTGCTACCCACTCGCTCGACGGTACGCGCCTGCAGCACAAGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCTTCGAGGAACAGGGTACCGCCGTTgg                                                                                                                                                     <  7:164172/147‑1 (MQ=255) cGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTgg                                                                                                                                                     >  2:106436/1‑149 (MQ=255) cGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTg                                                                                                                                                      <  2:5502/148‑1 (MQ=255)                  gcACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCg                                                                                                                                    >  5:312279/1‑149 (MQ=255)                    aCGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCgg                                                                                                                                   >  7:179455/1‑148 (MQ=255)                         gATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGc                                                                                                                             >  4:144782/1‑149 (MQ=255)                                      ttgttgCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGcccc                                                                                                                >  6:317745/1‑149 (MQ=255)                                        gttgCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGccc                                                                                                                 <  6:312279/146‑1 (MQ=255)                                        gttgCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTCGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCgg                                                                                                              >  1:190561/1‑149 (MQ=255)                                           gCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACAGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTa                                                                                                            >  8:26099/1‑148 (MQ=255)                                               cccACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑g                                                                                                        <  8:179455/147‑1 (MQ=255)                                               cccACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑Ggc                                                                                                      <  5:308056/149‑1 (MQ=255)                                               cccACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTCGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑Ggc                                                                                                      <  2:190561/149‑1 (MQ=255)                                                ccACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑Ggcg                                                                                                     >  8:257188/1‑149 (MQ=255)                                                ccACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTCGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑Ggcg                                                                                                     >  8:186798/1‑149 (MQ=255)                                                 cACGGGCTCGAAGGGACGCTCCTGCAGCACACGGAGCAGCTTGTCCTGCATCGGCAACGGCATGTCGACGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCATTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGTTAAA‑Ggcgc                                                                                                    <  8:253545/149‑1 (MQ=255)                                                      gcTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTtc                                                                                               <  1:106436/149‑1 (MQ=255)                                                        tCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGt                                                                                                                         <  3:140055/122‑1 (MQ=255)                                                        tCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGt                                                                                                                         >  4:140055/1‑122 (MQ=255)                                                          gAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTtctc                                                                                             <  7:26099/147‑1 (MQ=255)                                                          gAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTCGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGt                                                                                           <  7:186798/149‑1 (MQ=255)                                                                        gcagcaCACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTtctc                                                                                             <  4:119742/133‑1 (MQ=255)                                                                        gcagcaCACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCAGGCACGCCTCGTGATGGCCCCGGTAAA‑cgcg                                                                                                     >  3:182776/1‑125 (MQ=255)                                                                         cagcaCACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAACGGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑c                                                                           <  1:9878/148‑1 (MQ=255)                                                                             acacGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTc                                                                     <  3:144782/149‑1 (MQ=255)                                                                               acGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCca                                                                   >  8:194608/1‑149 (MQ=255)                                                                                 gcagcaGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCcag                                                                  >  4:117279/1‑148 (MQ=255)                                                                                 gcagcaGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCcag                                                                  >  8:230852/1‑148 (MQ=255)                                                                                  cagcagCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGa                                                                                     >  1:308239/1‑131 (MQ=255)                                                                                  cagcagCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGa                                                                                     <  2:308239/131‑1 (MQ=255)                                                                                    gcagcTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGaa                                                                                    >  7:23323/1‑130 (MQ=255)                                                                                    gcagcTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGAACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGaa                                                                                    <  8:23323/130‑1 (MQ=255)                                                                                        cTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGc                                                                                        >  8:188675/1‑122 (MQ=255)                                                                                        cTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGc                                                                                        <  7:188675/122‑1 (MQ=255)                                                                                        cTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCgg                                                          <  3:157677/149‑1 (MQ=255)                                                                                              cTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGa                                                    >  2:347951/1‑149 (MQ=255)                                                                                              cTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGa                                                    <  7:230852/149‑1 (MQ=255)                                                                                                      gCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATCgcgc                                              <  5:312733/147‑1 (MQ=255)                                                                                                      gCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATCgcgc                                              >  6:312733/1‑147 (MQ=255)                                                                                                       cAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATCGCGCcac                                           <  7:257188/149‑1 (MQ=255)                                                                                                        aaCGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATCGCGCca                                            >  7:238098/1‑147 (MQ=255)                                                                                                         aCGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAg                                         <  3:117279/149‑1 (MQ=255)                                                                                                              aTGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGAc                                    <  7:194608/149‑1 (MQ=255)                                                                                                                   gCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCa                               >  1:93596/1‑149 (MQ=255)                                                                                                                   gCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCa                               >  2:278682/1‑149 (MQ=255)                                                                                                                   gCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCAAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCa                               >  8:272200/1‑149 (MQ=255)                                                                                                                    ccGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCa                               <  1:317290/148‑1 (MQ=255)                                                                                                                    ccGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCAc                              >  6:189692/1‑149 (MQ=255)                                                                                                                    ccGATTTCGTCGAGGAACAGGGCACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCAc                              <  5:317745/149‑1 (MQ=255)                                                                                                                          tcgtcgAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATcgc                                                <  2:129800/125‑1 (MQ=255)                                                                                                                          tcgtcgAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATcgc                                                >  1:129800/1‑125 (MQ=255)                                                                                                                          tcgtcgAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCACGaatg                         <  5:189692/148‑3 (MQ=255)                                                                                                                               gAGGAACGGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCtt                   >  2:37714/1‑149 (MQ=255)                                                                                                                                        ggTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTCCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATcgc                                                <  6:171959/111‑1 (MQ=255)                                                                                                                                        ggTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTCCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATcgc                                                >  5:171959/1‑111 (MQ=255)                                                                                                                                             cgccgTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTTCCCCTCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTg        <  1:278682/149‑1 (MQ=255)                                                                                                                                             cgccgTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCgg                                                      <  4:315617/100‑1 (MQ=255)                                                                                                                                               ccgTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATg   <  7:272200/149‑1 (MQ=255)                                                                                                                                                cgTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATgg  <  1:347951/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      CGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAA‑GGCGCCCTTCTCGTGGCCGAACAGTT‑‑‑CGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGG  >  NC_002947/4964500‑4964791

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