Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	602,117	C→T	6.3%	H97Y (CAC→TAC)	ribE →	riboflavin synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	602,117	0	C	T	6.3%	191.5 / 2.2	63	H97Y (CAC→TAC)	ribE	riboflavin synthase
Reads supporting (aligned to +/- strand):  ref base C (33/26);  new base T (2/2);  total (35/28)
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.00e-01

CTCACGGCCGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCACCGTCGACGGCACCAGC  >  NC_002947/601973‑602263                                                                                                                                                 |                                                                                                                                                     cTCACGGCCGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCt                                                                                                                                                  >  6:21764/1‑149 (MQ=255)    aCGGCCGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑Gt                                                                                                                                              >  8:150943/1‑149 (MQ=255)      ggCCGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGTACCTG‑Gt                                                                                                                                              >  7:172165/1‑147 (MQ=255)      ggCCGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGTACCTG‑Gt                                                                                                                                              <  8:172165/147‑1 (MQ=255)          gTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCgg                                                                                                                                        <  6:37383/149‑1 (MQ=255)               aCTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCAc                                                                                                                                    <  4:99199/148‑1 (MQ=255)               aCTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GACAGCGGCCACg                                                                                                                                   >  6:131108/1‑149 (MQ=255)                   ccGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTc                                                                                                                                 <  8:219653/147‑1 (MQ=255)                          aCGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGtt                                                                                                                        <  3:241295/149‑1 (MQ=255)                          aCGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GACAGCGGCCACGTCGACGGGGtt                                                                                                                        <  5:131108/149‑1 (MQ=255)                            ggCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTc                                                                                                                                             <  5:87996/126‑1 (MQ=255)                            ggCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTc                                                                                                                                             >  6:87996/1‑126 (MQ=255)                             gCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGTTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGc                                                                                                                     >  5:241564/1‑149 (MQ=255)                                   gggCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAAt                                                                                                                >  6:24039/1‑148 (MQ=255)                                   gggCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAAt                                                                                                                >  4:678/1‑148 (MQ=255)                                     gCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCAt                                                                                                             <  2:222878/149‑1 (MQ=255)                                        gATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTc                                                                                                          <  3:678/149‑1 (MQ=255)                                            tCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCgcgc                                                                                                      <  6:314450/149‑1 (MQ=255)                                               gCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGc                                                                                                   <  7:150943/149‑1 (MQ=255)                                                    gAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATaa                                                                                              >  7:245442/1‑149 (MQ=255)                                                     aaaCCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAAc                                                                                             >  3:16124/1‑149 (MQ=255)                                                     aaaCCCTCAAGCGCACCGCCTTCATCGACCTCAACAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTAACCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGACCACGGGGTTGGCGAAATCATTACGCGCAGCGATAAc                                                                                             >  6:210649/1‑149 (MQ=255)                                                        cccTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGcc                                                                                          <  5:45568/149‑1 (MQ=255)                                                         ccTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGccc                                                                                         >  2:313146/1‑149 (MQ=255)                                                         ccTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCCCGCAGCGATAACGccc                                                                                         >  2:241283/1‑149 (MQ=255)                                                          cTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GACAGCGGCCACGTCGACGGGGATGGCGAAATCATTTCGCGCAGCGATAACGCCcg                                                                                        >  3:169791/1‑149 (MQ=255)                                                           tCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCAtt                                                                                                            >  6:109869/1‑128 (MQ=255)                                                           tCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCAtt                                                                                                            <  5:109869/128‑1 (MQ=255)                                                                 gcACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAAc                                                                                             <  6:116270/137‑1 (MQ=255)                                                                 gcACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAAc                                                                                             >  5:116270/1‑137 (MQ=255)                                                                 gcACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATc                                                                                 <  4:310265/149‑1 (MQ=255)                                                                           cATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCc‑g                                                                        <  6:15781/147‑1 (MQ=255)                                                                           cATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCc‑g                                                                        <  1:163615/147‑1 (MQ=255)                                                                                aCCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑gcgtgcgt                                                                 >  3:285016/1‑149 (MQ=255)                                                                                     aaGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCa                                                              <  5:24039/147‑1 (MQ=255)                                                                                        agCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCg                                                           >  3:54302/1‑147 (MQ=255)                                                                                         gCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAg                                                        <  1:313146/149‑1 (MQ=255)                                                                                          cGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAg                                                        >  6:195377/1‑148 (MQ=255)                                                                                                caAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTgg                                                 >  4:237400/1‑149 (MQ=255)                                                                                                  aGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGc                                                >  3:243326/1‑148 (MQ=255)                                                                                                    gTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCaa                                             <  5:195377/149‑1 (MQ=255)                                                                                                     tCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGTGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGcc                                               >  1:80103/1‑147 (MQ=255)                                                                                                         ccTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACa                                        >  7:84316/1‑149 (MQ=255)                                                                                                          cTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACAt                                       >  4:285744/1‑149 (MQ=255)                                                                                                              aaaaaGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGCCGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGAGCGTGCACCGAAGGAACTGGCCAAGTACATCGcc                                   >  5:216253/1‑149 (MQ=255)                                                                                                                 aaGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCcac                                <  3:285744/149‑1 (MQ=255)                                                                                                                   gCCCTGAGCCCCACTACCCGCCTGGGCGGTCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACaa                              <  4:169791/149‑1 (MQ=255)                                                                                                                      cTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAggg                           <  3:237400/149‑1 (MQ=255)                                                                                                                       tGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGt                          >  5:88044/1‑149 (MQ=255)                                                                                                                           ccccACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGa                      >  7:8388/1‑149 (MQ=255)                                                                                                                           ccccACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGa                      >  6:181161/1‑149 (MQ=255)                                                                                                                            cccACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCTGCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGa                      <  1:241283/149‑1 (MQ=255)                                                                                                                            cccACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGAt                     <  5:210649/149‑1 (MQ=255)                                                                                                                             ccACTACCCGCCTGGGCGGGCCCCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGAt                     >  5:323798/1‑148 (MQ=255)                                                                                                                             ccACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGAt                     >  3:219042/1‑148 (MQ=255)                                                                                                                              cACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCa                   >  2:47912/1‑149 (MQ=255)                                                                                                                                  gcccGCCCAGGCGGGTACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGAt                     <  6:116209/142‑1 (MQ=255)                                                                                                                                  gcccGCCCAGGCGGGTACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGAt                     >  5:116209/2‑143 (MQ=255)                                                                                                                                   cccGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCACCg                <  5:21764/147‑1 (MQ=255)                                                                                                                                    ccGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCAcc                 >  6:302594/1‑145 (MQ=255)                                                                                                                                    ccGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCAcc                 <  5:302594/145‑1 (MQ=255)                                                                                                                                    ccGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCACCGTc              >  1:31562/1‑148 (MQ=255)                                                                                                                                    ccGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCACCGTc              >  1:249055/1‑148 (MQ=255)                                                                                                                                           ggcgggcACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCACCGTCGACGGCa       <  4:54302/148‑1 (MQ=255)                                                                                                                                            gcgggcACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCACCGTCGACGGCAcc     <  8:84316/149‑1 (MQ=255)                                                                                                                                              gggcACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCACCGTCGACGGCACCAg   >  2:319975/1‑149 (MQ=255)                                                                                                                                               ggcACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCACCGTCGACGGCACCAGc  <  4:243326/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CTCACGGCCGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTG‑GTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCC‑GCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCACCGTCGACGGCACCAGC  >  NC_002947/601973‑602263

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