Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,590,016	G→T	6.1%	L179L (CTG→CTT)	nicP‑I →	porin‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,590,016	0	G	T	6.1%	209.0 / 3.1	66	L179L (CTG→CTT)	nicP‑I	porin‑like protein
Reads supporting (aligned to +/- strand):  ref base G (31/31);  new base T (2/2);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

GTGCAGGTGAACTGCAACCGAACCTGCCGGTGCTGACCTTCAGTGATATT‑CGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACG  >  NC_002947/3589871‑3590155                                                                                                                                                   |                                                                                                                                            gTGCAGGTGAACTGCAACCGAACCTGCCGGTGCTGACCTTCAGTGATATTGCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCAACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCg                                                                                                                                           >  1:230876/1‑149 (MQ=255)         gAACTGCAACCGAACCTGCCGGTGCTGACCTTCAGTGATATT‑CGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAgg                                                                                                                                  >  1:308559/1‑149 (MQ=255)           aCTGCAACCGAACCTGCCGGTGCTGACCTTCAGTGATATT‑CGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGGCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTg                                                                                                                                <  5:6870/149‑1 (MQ=255)             tGCAACCGAACCTGCCGGTGCTGACCTTCAGTGATATT‑CGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTTGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGc                                                                                                                              >  7:87496/1‑149 (MQ=255)               cAACCGAACCTGCCGGTGCTGACCTTCAGTGATATT‑CGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTTCGTTACGCAGGTgggtg                                                                                                                            <  4:256152/149‑4 (MQ=255)               cAACCGAACCTGCCGGTGCTGACCTTCAGTGATATT‑CGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCAACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGa                                                                                                                            <  2:230876/149‑1 (MQ=255)                aaCCGAACCTGCCGGTGCTGACCTTCAGTGATATT‑CGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTTCGTTACGCAGGTg                                                                                                                                <  8:68241/144‑1 (MQ=255)                aaCCGAACCTGCCGGTGCTGACCTTCAGTGATATT‑CGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTTCGTTACGCAGGTg                                                                                                                                >  7:68241/1‑144 (MQ=255)                    gAACCTGCCGGTGCTGACCTTCAGTGATATT‑CACCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGc                                                                                                                        >  1:188355/1‑148 (MQ=255)                       ccTGCCGGTGCTGACCTTCAGTGATATT‑CGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGc                                                                                                                        >  1:85957/1‑145 (MQ=255)                       ccTGCCGGTGCTGACCTTCAGTGATATT‑CGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGc                                                                                                                        >  4:69652/1‑145 (MQ=255)                            cGGTGCTGACCTTCAGTGATATT‑CGCCTGTTGCCACCCACCTATCAGGGCGTGAGTATGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAAc                                                                                                               >  3:232367/1‑149 (MQ=255)                            cGGTGCCGACCTTCAGTGATATT‑CGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTCAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAAc                                                                                                               <  1:217524/149‑1 (MQ=255)                                   gACCTTCAGTGATATT‑CGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGc                                                                                                        >  5:259835/1‑149 (MQ=255)                                          aGTGATATT‑CGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCa                                                                                                 >  4:221602/1‑149 (MQ=255)                                                atT‑CGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGcc                                                                                           >  1:78051/1‑149 (MQ=255)                                                atT‑CGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGcc                                                                                           >  4:133781/1‑149 (MQ=255)                                                  t‑CGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGa                                                                                         >  4:51115/1‑149 (MQ=255)                                                  t‑CGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGa                                                                                         >  2:231498/1‑149 (MQ=255)                                                         gTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCaa                                                                                                                      >  5:255907/1‑114 (MQ=255)                                                         gTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCaa                                                                                                                      <  6:255907/114‑1 (MQ=255)                                                          ttGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGTAAATCGACTGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCTATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATTCGa                                                                                    <  6:102464/147‑1 (MQ=255)                                                             ccacccacCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAgg                                                                                <  6:288869/148‑1 (MQ=255)                                                              cacccaccTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCAGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCAGATGCGACAGGcc                                                                              >  8:140934/1‑149 (MQ=255)                                                                  caccTATCAGGGCGTGAGTTTGAGCTTGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGaaa                                                                          <  8:87496/149‑1 (MQ=255)                                                                   accTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAg                                                                         >  3:279616/1‑149 (MQ=255)                                                                          aGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGAcgc                                                                   <  6:259835/148‑1 (MQ=255)                                                                               gTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGc                                                                                            <  7:31785/118‑1 (MQ=255)                                                                               gTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGc                                                                                            >  8:31785/1‑118 (MQ=255)                                                                               gTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCaa                                                             <  3:51115/149‑1 (MQ=255)                                                                               gTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCaa                                                             <  2:308559/149‑1 (MQ=255)                                                                                   gTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTAt                                                         <  1:231498/149‑1 (MQ=255)                                                                                    tttGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTAtg                                                        >  7:255962/1‑149 (MQ=255)                                                                                      tGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTAtgt                                                       <  2:85957/148‑1 (MQ=255)                                                                                       gAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATgtgg                                                     >  5:230539/1‑149 (MQ=255)                                                                                        aGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTAt                                                         >  2:282951/1‑144 (MQ=255)                                                                                        aGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTAt                                                         <  1:282951/144‑1 (MQ=255)                                                                                        aGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATgtgg                                                     <  1:10771/148‑1 (MQ=255)                                                                                        aGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATgtgg                                                     <  3:221602/148‑1 (MQ=255)                                                                                         gCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGAGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATgtgggt                                                   >  6:81746/1‑149 (MQ=255)                                                                                             gtcgGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATgtgggtgg                                                 <  3:133781/147‑1 (MQ=255)                                                                                              tcgGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGa                                              >  4:123402/1‑149 (MQ=255)                                                                                                 gAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTa                                           >  4:96924/1‑149 (MQ=255)                                                                                                           ggCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGAAAGATGAACGCGATGCGCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAAt                                 <  4:232367/149‑1 (MQ=255)                                                                                                              tGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGAGAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGc                               >  1:252123/1‑148 (MQ=255)                                                                                                                aCTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGccaa                            <  2:252123/149‑1 (MQ=255)                                                                                                                 cTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGccaa                            <  4:279616/148‑1 (MQ=255)                                                                                                                   aTGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGAt                                                                                                                           >  4:50697/1‑51 (MQ=255)                                                                                                                   aTGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGAt                                                                                                                           <  3:50697/51‑1 (MQ=255)                                                                                                                    tGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGccaaccaa                        <  3:123402/149‑1 (MQ=255)                                                                                                                     gCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGACGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGccaagatcg                       >  3:160016/1‑144 (MQ=255)                                                                                                                     gCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTATGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGcc                              <  4:160016/142‑1 (MQ=255)                                                                                                                      cAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCaaa                       <  6:230539/148‑1 (MQ=255)                                                                                                                      cAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCaaa                       <  7:140934/148‑1 (MQ=255)                                                                                                                             ggCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGt                >  5:225678/1‑148 (MQ=255)                                                                                                                             ggCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTc               <  3:96924/149‑1 (MQ=255)                                                                                                                              gCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGt                >  5:301838/1‑147 (MQ=255)                                                                                                                                aCTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGc            >  7:251543/1‑149 (MQ=255)                                                                                                                                  ttGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCt          >  3:97661/1‑149 (MQ=255)                                                                                                                                     aaaaGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGtt                                     >  6:30330/1‑119 (MQ=255)                                                                                                                                     aaaaGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGtt                                     <  5:30330/119‑1 (MQ=255)                                                                                                                                      aaaGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTgg      <  3:69652/149‑1 (MQ=255)                                                                                                                                      aaaGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTgg      <  2:78051/149‑1 (MQ=255)                                                                                                                                      aaaGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTg       <  8:251543/148‑1 (MQ=255)                                                                                                                                          tcgccaccTTCGTTACGCAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACg  >  8:20639/4‑149 (MQ=255)                                                                                                                                          cacccaccTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACg  <  4:97661/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                            GTGCAGGTGAACTGCAACCGAACCTGCCGGTGCTGACCTTCAGTGATATT‑CGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACG  >  NC_002947/3589871‑3590155

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