Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,549,900	C→T	9.2%	V100V (GTC→GTT)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,549,900	0	C	T	9.2%	179.4 / 6.9	65	V100V (GTC→GTT)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base C (28/31);  new base T (3/3);  total (31/34)
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.60e-02

GGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACC  >  NC_002947/1549753‑1550047                                                                                                                                                      |                                                                                                                                                    ggACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCa                                                                                                                                                    <  2:38494/149‑1 (MQ=255)     ggCGTTTCCGTCGCTAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACg                                                                                                                                                 >  4:237645/1‑148 (MQ=255)       cGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAg                                                                                                                                              <  3:147349/149‑1 (MQ=255)            ccGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCt                                                                                                                                          <  4:11933/148‑1 (MQ=255)             cGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCTTCGTTGACGATGGCCt                                                                                                                                          >  4:187894/1‑147 (MQ=255)             cGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCTTCGTTGACGATGGCCTGa                                                                                                                                        <  3:187894/149‑1 (MQ=255)               tCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGTC‑‑TCAGGCCATCGTCAACGAAGGCCTGaaa                                                                                                                                      <  2:305800/149‑1 (MQ=255)               tCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGaaa                                                                                                                                      <  7:71129/149‑1 (MQ=255)                 gCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTTTCAGGCCTTCGTTGACGATGGCCTGagc                                                                                                                                      <  5:146825/149‑3 (MQ=255)                              aGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCgg                                                                                                                        >  5:203643/1‑148 (MQ=255)                                cTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCAt                                                                                                                     <  1:171543/149‑1 (MQ=255)                                cTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCAt                                                                                                                     >  3:160190/1‑149 (MQ=255)                                   aaaGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCt                                                                                                                                          >  6:93007/1‑125 (MQ=255)                                   aaaGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGaa                                                                                                                  >  7:51030/1‑149 (MQ=255)                                   aaaGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGa                                                                                                                   <  6:167896/148‑1 (MQ=255)                                   aaaGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCCCTGCCGGCATGaa                                                                                                                  >  5:295174/1‑149 (MQ=255)                                    aaGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAAc                                                                                                                 <  8:186728/149‑1 (MQ=255)                                           tttCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGa                                                                                                                   <  2:71218/140‑1 (MQ=255)                                           tttCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGa                                                                                                                   >  1:71218/1‑140 (MQ=255)                                             tCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGAc                                                                                                        >  2:304363/1‑149 (MQ=255)                                                   aCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCTTCGTTGACGATGGCCTg                                                                                                                                         >  7:7815/1‑110 (MQ=255)                                                   aCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCTTCGTTGACGATGGCCTg                                                                                                                                         <  8:7815/110‑1 (MQ=255)                                                        ggCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGcgcg                                                                                              <  1:304363/148‑1 (MQ=255)                                                             ccAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATc                                                                                         >  3:78187/1‑148 (MQ=255)                                                                 cTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACaa                                                                                    <  4:78187/149‑1 (MQ=255)                                                                     tCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGcc                                                                                <  4:74644/149‑1 (MQ=255)                                                                       aaGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGc                                                                                 <  6:177690/146‑1 (MQ=255)                                                                       aaGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGc                                                                                 >  5:177690/1‑146 (MQ=255)                                                                        aGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCAct                                                                             <  6:203643/149‑1 (MQ=255)                                                                         ggAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCAct                                                                             >  7:77681/1‑148 (MQ=255)                                                                               tGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATTGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGcc                                                                       <  4:160190/148‑1 (MQ=255)                                                                                gCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTc                                                                                                                                <  8:210632/90‑1 (MQ=255)                                                                                gCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTAACGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCCCGGTATCCACAAGGCCACTGCTGCcgt                                                                     >  5:305894/1‑149 (MQ=255)                                                                                  ttCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcgtcg                                                                   >  4:221345/1‑149 (MQ=255)                                                                                  ttCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcgtcg                                                                   >  8:268087/1‑149 (MQ=255)                                                                                  ttCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCAGAAGCGAGGTATCGACAAGGCCACTGCTGCcgtcg                                                                   >  3:197289/1‑149 (MQ=255)                                                                                    ccAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcgtcgt                                                                  >  6:321962/1‑148 (MQ=255)                                                                                    ccAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcgtcgt                                                                  <  5:321962/148‑1 (MQ=255)                                                                                      aaGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGc                                                               >  3:79755/1‑149 (MQ=255)                                                                                              cgacgCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTga                                                       <  8:77681/149‑1 (MQ=255)                                                                                                acgCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTgaag                                                     <  8:51030/149‑1 (MQ=255)                                                                                                  gCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTgaagaa                                                   >  1:236011/1‑149 (MQ=255)                                                                                                   cgtccGCTGACCGCCCCCCCCCCCCTCCCGGTCTGTC‑‑TCCGGCCTTCGCCCACGCAGGCCTGATAGCCGGCCGTGCCGGCTTGCACCCGCTGGACCTGATGCGCGGTTTCGTCAATGCCTCTGCTGCCGTCGTTGCCGAACTGAAGAAc                                                  <  2:285042/146‑1 (MQ=255)                                                                                                   atgacgttGACGGCACAACCACAGCTACAGTTCTTGC‑‑TCACGCCATCGTCAACTAACGCCTTTAAGCCGTCGCTGCCGGCATGAACCCGATTGACCTTAAGCGCGTTATCTACAAGGCCACTGCTGCCGTCGTTGCCGTACTGAAGAAc                                                  <  2:238661/142‑1 (MQ=255)                                                                                                      ccGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCt                                                <  3:221345/148‑1 (MQ=255)                                                                                                      ccGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTg                                               >  6:93764/1‑149 (MQ=255)                                                                                                           gACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCaa                                          >  2:5392/1‑149 (MQ=255)                                                                                                              ggCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGc                                        >  8:115786/1‑148 (MQ=255)                                                                                                              ggCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGc                                        >  4:183728/1‑148 (MQ=255)                                                                                                               gcacacccaccGCTACCGTTCTGGC‑‑TCAGGCCATAGTCAACGAAGGCCTGAAAGCCGTAGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGACACTGCTGCCGTCGTTTCAGAACTGAAGAACCTGTCCAAGACa                                      >  8:264729/7‑149 (MQ=255)                                                                                                                 accaccaccGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATg                                    >  1:243640/1‑149 (MQ=255)                                                                                                                  ccaccaccGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATgc                                   <  3:183728/149‑1 (MQ=255)                                                                                                                     caaccGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCAGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACGGAAGAACCTGTCCAAGCCATGCGcc                                >  7:146865/3‑149 (MQ=255)                                                                                                                       accGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGa                              <  1:5392/149‑1 (MQ=255)                                                                                                                         cGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGcc                                                              <  7:57853/115‑1 (MQ=255)                                                                                                                         cGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGcc                                                              >  8:57853/1‑115 (MQ=255)                                                                                                                            tACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTcc                          >  5:21452/1‑148 (MQ=255)                                                                                                                              ccGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAg                       <  4:79755/149‑1 (MQ=255)                                                                                                                                      gC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATc                 <  5:63362/147‑1 (MQ=255)                                                                                                                                          tCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGcc              >  6:79862/1‑148 (MQ=255)                                                                                                                                           ctgggaaTCTTCAACGACGCCCTGAAAGCGGTCGGTGCCGGCTTGTACCCGATGGACCTGATGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGATCTGAAGATCCTGTCCTAGCCATGCGCCGACTCCAAGGCCATCGccc             <  4:197289/142‑1 (MQ=255)                                                                                                                                             ggCCTTCGTTGACGATGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCagg          >  2:215563/1‑149 (MQ=255)                                                                                                                                                  tCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCaggtaggt     <  2:243640/149‑1 (MQ=255)                                                                                                                                                    gTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTAc   <  5:79862/149‑1 (MQ=255)                                                                                                                                                     tCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTAcc  <  6:21452/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    GGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACC  >  NC_002947/1549753‑1550047

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