Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,549,895	A→T	10.9%	I99F (ATC→TTC)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,549,895	0	A	T	10.9%	164.5 / 9.3	62	I99F (ATC→TTC)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base A (29/26);  new base T (3/4);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGG  >  NC_002947/1549749‑1550039                                                                                                                                                     |                                                                                                                                                 ccAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGCCGCTGCCGGTGACGGCACCACCACCGCAACCGTTCTGGC‑‑TCAGGCCATc                                                                                                                                                >  1:285042/1‑149 (MQ=255)     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                                                                         >  5:177690/1‑146 (MQ=255)                                                                           aaGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGc                                                                         <  6:177690/146‑1 (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                                                          <  5:321962/148‑1 (MQ=255)                                                                                        ccAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcgtcgt                                                          >  6:321962/1‑148 (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                                >  4:183728/1‑148 (MQ=255)                                                                                                                  ggCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGc                                >  8:115786/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)                                                                                                                                                     |                                                                                                                                                 CCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGC‑‑TCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGG  >  NC_002947/1549749‑1550039

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