Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,549,895	A→T	14.3%	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	14.3%	171.6 / 17.5	70	I99F (ATC→TTC)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base A (32/28);  new base T (4/6);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

ACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGG  >  NC_002947/1549748‑1550039                                                                                                                                                     |                                                                                                                                                   aCCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCCCCCCCGCTACCGTTCTGGCTCCGGCCa                                                                                                                                                    >  4:68143/1‑148 (MQ=255) aCCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCa                                                                                                                                                    <  3:68143/148‑1 (MQ=255)  ccAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAATGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATc                                                                                                                                                  <  1:200525/149‑1 (MQ=255)   cAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCg                                                                                                                                                 >  5:6739/1‑149 (MQ=255)            cGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCaa                                                                                                                                             <  1:61127/144‑1 (MQ=255)            cGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCaa                                                                                                                                             >  2:61127/1‑144 (MQ=255)             gTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGCGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAg                                                                                                                                        <  1:243772/149‑1 (MQ=255)              tttCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGc                                                                                                                                      <  4:137396/149‑1 (MQ=255)                 ccGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCt                                                                                                                                    >  6:268964/1‑148 (MQ=255)                  cGTCTCCATAGAAATCGAGCTGAAAGTCGCTTTCGAAAACATGGGCGCGCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTg                                                                                                                                   <  5:22508/148‑1 (MQ=255)                  cGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGa                                                                                                                                  <  1:225112/149‑1 (MQ=255)                  cGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGa                                                                                                                                  <  8:103988/149‑1 (MQ=255)                  cGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGa                                                                                                                                  <  2:248817/149‑1 (MQ=255)                    tCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGaaa                                                                                                                                <  3:267509/149‑1 (MQ=255)                      gCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑g                                                                                                                              <  4:84130/148‑1 (MQ=255)                       ccAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGc                                                                                                                                      >  5:107755/1‑140 (MQ=255)                       ccAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGc                                                                                                                                      <  6:107755/140‑1 (MQ=255)                        cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCg                                                                                                                          >  4:283507/1‑149 (MQ=255)                        cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCg                                                                                                                          >  1:49127/1‑149 (MQ=255)                         aaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGt                                                                                                                         >  3:209907/1‑149 (MQ=255)                              aaTCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTg                                                                                                                    >  8:185194/1‑149 (MQ=255)                               aTCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGCCGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGa                                                                                                                                  <  7:149914/136‑1 (MQ=255)                               aTCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGa                                                                                                                                  >  8:153394/1‑136 (MQ=255)                               aTCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGa                                                                                                                                  >  8:149914/1‑136 (MQ=255)                               aTCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGa                                                                                                                                  <  7:153394/136‑1 (MQ=255)                                  gAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCg                                                                                                                 <  2:7949/148‑1 (MQ=255)                                  gAGCGGAAACACGCATTCGAAAACAAGGGCGCCCAGCCGGGCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCgg                                                                                                                <  1:45543/149‑1 (MQ=255)                                      tGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGa                                                                                                                                  >  2:308511/1‑129 (MQ=255)                                      tGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGa                                                                                                                                  <  1:308511/129‑1 (MQ=255)                                      tGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCAt                                                                                                             >  6:78013/1‑148 (MQ=255)                                       gAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATg                                                                                                            >  8:145403/1‑148 (MQ=255)                                        aaaGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGa                                                                                                           <  5:265865/148‑1 (MQ=255)                                                   cGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGAcc                                                                                               >  3:212517/1‑149 (MQ=255)                                                   cGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGAc                                                                                                <  5:78013/148‑1 (MQ=255)                                                   cGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGAc                                                                                                <  4:134400/148‑1 (MQ=255)                                                   cGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGAc                                                                                                <  6:105338/148‑1 (MQ=255)                                                     aaaaCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCt                                                                                              >  8:96599/1‑148 (MQ=255)                                                      aaaCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTg                                                                                             >  1:293669/1‑148 (MQ=255)                                                       aaCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGAccgg                                                                                             >  2:266681/1‑145 (MQ=255)                                                              gcgcCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAATGAGCGGCGCGGCCGGCATGAACCCGATGGACCTGAAGcgcg                                                                                      >  1:109494/1‑149 (MQ=255)                                                              gcgcCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCgg                                                                                     >  1:238408/1‑148 (MQ=255)                                                                 cccAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATc                                                                                 <  2:49127/149‑1 (MQ=255)                                                                  ccAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGACGCTGCCGGCATCAACCCGATGGACCTGAAGCGCGGTATc                                                                                 >  2:160947/1‑148 (MQ=255)                                                                      cTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACaa                                                                            <  4:209907/149‑1 (MQ=255)                                                                         gTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGc                                                                         >  3:309246/1‑149 (MQ=255)                                                                            aaGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCAc                                                                      >  3:242171/1‑149 (MQ=255)                                                                              ggAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCActg                                                                    >  7:210496/1‑149 (MQ=255)                                                                                 aGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCActgctg                                                                 >  2:106441/1‑149 (MQ=255)                                                                                     gCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcgt                                                             >  7:171183/1‑149 (MQ=255)                                                                                       ttCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCt                                                                                              >  6:147362/1‑114 (MQ=255)                                                                                       ttCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCt                                                                                              <  5:147362/114‑1 (MQ=255)                                                                                        tCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcgtcgt                                                          >  4:33375/1‑149 (MQ=255)                                                                                         ccAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGc                                                                         <  2:164/133‑1 (MQ=255)                                                                                         ccAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGc                                                                         >  1:164/1‑133 (MQ=255)                                                                                          cAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAg                                                                           >  2:201048/1‑130 (MQ=255)                                                                                          cAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAg                                                                           <  1:201048/130‑1 (MQ=255)                                                                                                 aaCGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAACGCCACTCCTGCAGTCGTTGCCGacat                                                 >  8:275401/1‑146 (MQ=255)                                                                                                   cgacgCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTg                                                >  3:91731/1‑148 (MQ=255)                                                                                                       gCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTgaagaa                                           <  2:293669/149‑1 (MQ=255)                                                                                                        cTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAAc                                          <  7:145403/149‑1 (MQ=255)                                                                                                           ccG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAAcc                                         <  4:309246/147‑1 (MQ=255)                                                                                                            cG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCAGGCATGAACCCGATGGACCTGAAGCGCGATATCGACAAGGCCCCTGCTGACGTCGTTGCCGAACTGAAGAACCTga                                      >  3:263527/1‑148 (MQ=255)                                                                                                                 gACGGCCCCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGCACCCGATGGCCCTGAAGCGCGGTATCGGCAAGGCCGCTGCTGCcg                                                              >  4:93382/1‑121 (MQ=255)                                                                                                                 gACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCTACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGATCACGATGGACCTGTAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCaa                                  <  5:229635/149‑1 (MQ=255)                                                                                                                 gACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcg                                                              <  3:93382/121‑1 (MQ=255)                                                                                                                  aCGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCaa                                  <  4:212517/148‑1 (MQ=255)                                                                                                                           ccaccGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGcc                        <  5:268964/149‑1 (MQ=255)                                                                                                                                    ccGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAg               <  4:91731/149‑1 (MQ=255)                                                                                                                                         cTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAg               <  8:161480/144‑1 (MQ=255)                                                                                                                                         cTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAg               >  7:161480/1‑144 (MQ=255)                                                                                                                                            tttCAGGCCTTCGTTGACGATGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGc       >  5:233010/3‑149 (MQ=255)                                                                                                                                                 ggCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCagg  >  8:39495/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                   ACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGG  >  NC_002947/1549748‑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: ‑