Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,633,552	C→G	7.8%	G167G (GGC→GGG)	lepB →	signal peptidase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,633,552	0	C	G	7.8%	203.3 / 6.4	64	G167G (GGC→GGG)	lepB	signal peptidase I
Reads supporting (aligned to +/- strand):  ref base C (26/33);  new base G (2/3);  total (28/36)
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.36e-01

AGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTG  >  NC_002947/1633413‑1633699                                                                                                                                              |                                                                                                                                                    aGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGtt                                                                                                                                            <  1:307524/149‑1 (MQ=255)  ggTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTg                                                                                                                                           <  7:33254/149‑1 (MQ=255)    tCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTgac                                                                                                                                         <  5:112434/149‑3 (MQ=255)      ggTGACCCGCAACGCGGGGATGTAATGGGGTTCCGCTACCCGAGTGACCCGA‑ATTTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGa                                                                                                                                       <  6:61025/149‑1 (MQ=255)       gTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGaa                                                                                                                                      >  7:92404/1‑149 (MQ=255)        tGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGaa                                                                                                                                      >  7:241890/1‑148 (MQ=255)          aCCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTCACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTGACGaaca                                                                                                                                    <  1:240428/149‑1 (MQ=255)            ccGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACAGCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACt                                                                                                                                 >  5:261825/1‑149 (MQ=255)                aaagcgTGGATGTAAAGGTGTTCCGCTACCCGAGTGAACCTA‑ATCTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtggtg                                                                                                                             <  8:125930/146‑1 (MQ=255)                 aCGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtggtgg                                                                                                                            <  2:139649/149‑1 (MQ=255)                   gcgGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGc                                                                                                                          <  6:196612/149‑1 (MQ=255)                   gcgGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGGGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGc                                                                                                                          <  7:53322/149‑1 (MQ=255)                          tGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGc                                                                                                                   <  6:228959/149‑1 (MQ=255)                           gTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGccc                                                                                                                                               >  5:257418/1‑120 (MQ=255)                           gTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGccc                                                                                                                                               <  6:257418/120‑1 (MQ=255)                             aaTGGTGTTCCGCTACCCGAGTGACCCGAGATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCg                                                                                                                 >  1:24657/1‑149 (MQ=255)                              aTGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCgg                                                                                                                <  2:24657/148‑1 (MQ=255)                               tGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGccc                                                                                                                                               <  2:242083/116‑1 (MQ=255)                               tGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGccc                                                                                                                                               >  1:242083/1‑116 (MQ=255)                                   gtTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCg                                                                                                                 >  7:179525/1‑142 (MQ=255)                                   gtTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCg                                                                                                                 <  8:179525/142‑1 (MQ=255)                                    ttCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTAccc                                                                                                          <  7:169902/148‑1 (MQ=255)                                    ttCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCt                                                                                                         >  3:274003/1‑149 (MQ=255)                                      ccGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCAGGGTACCCTgg                                                                                                       >  5:77665/1‑149 (MQ=255)                                         cTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCa                                                                                                    <  7:247662/149‑1 (MQ=255)                                                    aCCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTg                                                                                         >  2:281240/1‑149 (MQ=255)                                                      ccGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTa                                                                                       >  1:190218/1‑149 (MQ=255)                                                                aaCTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCt                                                                                               <  1:302910/132‑1 (MQ=255)                                                                aaCTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCt                                                                                               >  2:302910/1‑132 (MQ=255)                                                                 aCTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACTGCCAGCCCGTTGCCTAACAACTGGTTGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAg                                                                             <  5:157235/149‑1 (MQ=255)                                                                 aCTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAg                                                                             <  2:158514/149‑1 (MQ=255)                                                                    aCATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTAccc                                                                                                          <  7:112609/117‑1 (MQ=255)                                                                    aCATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTAccc                                                                                                          >  8:112609/1‑117 (MQ=255)                                                                    aCATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTc                                                                          >  4:6729/1‑149 (MQ=255)                                                                           gCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCACCGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCgg                                                                   <  2:299774/149‑1 (MQ=255)                                                                               gtgGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCg                                                                >  5:207187/1‑148 (MQ=255)                                                                                tgGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGaa                                                              <  1:281240/149‑1 (MQ=255)                                                                                  gTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAaca                                                            >  7:241687/1‑149 (MQ=255)                                                                                   ttGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAacac                                                           <  8:241890/149‑1 (MQ=255)                                                                                     ggCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACAcc                                                          <  8:241687/148‑1 (MQ=255)                                                                                            cGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGc                                                  <  6:207187/149‑1 (MQ=255)                                                                                                    gTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAAt                                          <  2:190218/149‑1 (MQ=255)                                                                                                     tGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAAt                                          >  2:10343/1‑148 (MQ=255)                                                                                                        tCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAAt                                          <  5:284683/145‑1 (MQ=255)                                                                                                        tCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAAt                                          >  6:284683/1‑145 (MQ=255)                                                                                                          cGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGcc                                     <  8:92404/148‑1 (MQ=255)                                                                                                          cGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCg                                    >  1:285573/1‑149 (MQ=255)                                                                                                                cacCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAac                                                             <  8:76901/119‑1 (MQ=255)                                                                                                                   cAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGATCACCTGATCCGCAAGGAAATGAGACGTTATCGCa                            >  3:245023/1‑149 (MQ=255)                                                                                                                      cGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGaaaa                                                                              >  6:242232/1‑96 (MQ=255)                                                                                                                      cGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGaaaa                                                                              <  5:242232/96‑1 (MQ=255)                                                                                                                       gACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATgcc                        >  5:45337/1‑149 (MQ=255)                                                                                                                              ggCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACAcc                                                          >  3:214439/1‑108 (MQ=255)                                                                                                                                   ttCGTCAACGGCCAGCCCGTTTCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACcagc             <  2:285573/148‑1 (MQ=255)                                                                                                                                     cGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACcagcag           <  1:10343/148‑1 (MQ=255)                                                                                                                                      ggcAACGGGCTGGCCGTTGACGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTg         >  7:131846/3‑149 (MQ=255)                                                                                                                                      ggcAACGGGCTGGCCGTTGACGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTg         <  8:131846/147‑1 (MQ=255)                                                                                                                                      ggcAACGGGCTGGCCGTTGACGAACAACTGGTGGGCAGCGAGCCGGGTACCCT‑GGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGCTCCGCAAGGAAATGAGCAGTTATCGCATGCCGCCGGACCAGCAGTg         >  2:272472/3‑148 (MQ=255)                                                                                                                                      gTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCGGCAGTg         >  8:61947/1‑149 (MQ=255)                                                                                                                                      gTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTg         <  6:261825/149‑1 (MQ=255)                                                                                                                                       tCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCACAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTg         >  5:329427/1‑148 (MQ=255)                                                                                                                                         aaCGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGAc      >  3:185825/1‑149 (MQ=255)                                                                                                                                            ggCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGt   <  4:274003/149‑1 (MQ=255)                                                                                                                                             gCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGAGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTg  <  3:6729/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                    AGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTG  >  NC_002947/1633413‑1633699

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