Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,549,906	A→T	14.2%	E102D (GAA→GAT)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,549,906	0	A	T	14.2%	167.1 / 15.1	69	E102D (GAA→GAT)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base A (32/27);  new base T (4/6);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCT  >  NC_002947/1549760‑1550053                                                                                                                                                    |                                                                                                                                                      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                                                                                                                                                    <  6:107755/140‑1 (MQ=255)           ccAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGc                                                                                                                                                    >  5:107755/1‑140 (MQ=255)            cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCg                                                                                                                                        >  1:49127/1‑149 (MQ=255)            cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCg                                                                                                                                        >  4:283507/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                                                                                                              <  6:105338/148‑1 (MQ=255)                                       cGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGAc                                                                                                              <  4:134400/148‑1 (MQ=255)                                       cGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGAc                                                                                                              <  5:78013/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                                                                                       >  1:164/1‑133 (MQ=255)                                                                             ccAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGc                                                                                       <  2:164/133‑1 (MQ=255)                                                                              cAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAg                                                                                         <  1:201048/130‑1 (MQ=255)                                                                              cAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAg                                                                                         >  2:201048/1‑130 (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                             >  7:161480/1‑144 (MQ=255)                                                                                                                             cTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAg                             <  8:161480/144‑1 (MQ=255)                                                                                                                                tttCAGGCCTTCGTTGACGATGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGc                     >  5:233010/3‑149 (MQ=255)                                                                                                                                     ggCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCagg                >  8:39495/1‑149 (MQ=255)                                                                                                                                             tCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTAcc        <  6:6739/149‑1 (MQ=255)                                                                                                                                              cAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCa       >  3:235182/1‑149 (MQ=255)                                                                                                                                              cAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCa       >  5:98279/1‑149 (MQ=255)                                                                                                                                                   aaGGCGTGAAA‑G‑CCGTCGCTGGCGGCATTAAACCGATGGACCTGAAGCGCGTTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAATAAGCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATctct  <  4:242171/149‑1 (MQ=255)                                                                                                                                                   aaGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATctct  <  8:171183/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCTCT  >  NC_002947/1549760‑1550053

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