Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,549,901	A→G	14.8%	N101D (AAC→GAC)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,549,901	0	A	G	14.8%	154.4 / 19.2	69	N101D (AAC→GAC)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base A (32/26);  new base G (4/6);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

CGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCA  >  NC_002947/1549759‑1550048                                                                                                                                                |                                                                                                                                                      cGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCaa                                                                                                                                                      >  2:61127/1‑144 (MQ=255) cGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCaa                                                                                                                                                      <  1:61127/144‑1 (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                                                                                                                                   >  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                                                                                                                                           <  7:153394/136‑1 (MQ=255)                    aTCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGa                                                                                                                                           >  8:149914/1‑136 (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                                                                                                         <  6:105338/148‑1 (MQ=255)                                        cGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGAc                                                                                                         <  4:134400/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  >  5:98279/1‑149 (MQ=255)                                                                                                                                               cAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCa  >  3:235182/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                      CGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCG‑GTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAA‑G‑CCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCA  >  NC_002947/1549759‑1550048

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