Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,633,556	C→G	8.1%	P169A (CCC→GCC)	lepB →	signal peptidase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,633,556	0	C	G	8.1%	190.8 / 6.5	62	P169A (CCC→GCC)	lepB	signal peptidase I
Reads supporting (aligned to +/- strand):  ref base C (25/32);  new base G (2/3);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

AGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTGCCTGC  >  NC_002947/1633413‑1633704                                                                                                                                                  |                                                                                                                                                     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                                                                                                                                                    <  6:257418/120‑1 (MQ=255)                           gTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGccc                                                                                                                                                    >  5:257418/1‑120 (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                                                                                                                      <  8:179525/142‑1 (MQ=255)                                   gtTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCg                                                                                                                      >  7:179525/1‑142 (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                                               >  6:284683/1‑145 (MQ=255)                                                                                                        tCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAAt                                               <  5:284683/145‑1 (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)                                                                                                                                            ggacAGCACGTTGCAGAACAACTGGTGGGCAGCGAGCCGGGTACCCTTGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGTCGCGGCCGAACACCTGATCCGCAATGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGt        <  4:245023/146‑1 (MQ=255)                                                                                                                                            ggCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGt        <  4:274003/149‑1 (MQ=255)                                                                                                                                             gCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGAGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTg       <  3:6729/149‑1 (MQ=255)                                                                                                                                                aGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTCTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGCAAATGAGCAGTTATCGCATGCCGCCGCACCAGCAGTGGACCGTCCCt    >  2:123343/1‑149 (MQ=255)                                                                                                                                                  cccGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTGCCTGc  <  4:185825/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     AGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCT‑ACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTGCCTGC  >  NC_002947/1633413‑1633704

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