Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,893,482	G→A	5.9%	E282E (GAG→GAA)	rarD →	chloramphenicol‑sensitive protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,893,482	0	G	A	5.9%	207.3 / 2.0	68	E282E (GAG→GAA)	rarD	chloramphenicol‑sensitive protein
Reads supporting (aligned to +/- strand):  ref base G (37/27);  new base A (2/2);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

TGCTGCCATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAG  >  NC_002947/3893334‑3893625                                                                                                                                                      |                                                                                                                                                 tgctgcCATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAg                                                                                                                                                  <  5:160556/149‑1 (MQ=255)  gctgcCATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAgg                                                                                                                                                 >  8:125867/1‑149 (MQ=255)    tgcCATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGc                                                                                                                                               >  7:341652/1‑149 (MQ=255)      ccATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCTGTGGCGTTGCTGCTGGGCGAAACAATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCt                                                                                                                                              >  1:283398/1‑149 (MQ=255)      ccATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGc                                                                                                                                               <  5:373473/147‑1 (MQ=255)      ccATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGc                                                                                                                                               <  2:34058/147‑1 (MQ=255)        aTTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTc                                                                                                                                            >  6:91234/1‑148 (MQ=255)        aTTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTc                                                                                                                                            <  5:91234/148‑1 (MQ=255)            aGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGc                                                                                                                                        <  7:125867/148‑1 (MQ=255)             gCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAc                                                                                                                                      >  1:9877/1‑149 (MQ=255)             gCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAAGCCCTCAAGCAc                                                                                                                                      <  1:371793/149‑1 (MQ=255)              ccTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAcc                                                                                                                                     >  7:193678/1‑149 (MQ=255)              ccTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAAGCCCTCAAGCAcc                                                                                                                                     <  4:213530/149‑1 (MQ=255)               cTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCt                                                                                                                                    <  2:283398/149‑1 (MQ=255)                tGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTg                                                                                                                                   >  1:236127/1‑149 (MQ=255)                  ttCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGc                                                                                                                                  >  5:163280/1‑148 (MQ=255)                      gCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGt                                                                                                                             >  7:206513/1‑149 (MQ=255)                      gCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGt                                                                                                                             <  7:253307/149‑1 (MQ=255)                        cTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCa                                                                                                                                           >  2:271638/1‑133 (MQ=255)                        cTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCa                                                                                                                                           <  1:271638/133‑1 (MQ=255)                         tGCTCAGCTACGTCGAACCCGTGCTTCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACATACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTAAGTccg                                                                                                                          >  2:260219/1‑147 (MQ=255)                           cTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCg                                                                                                                        >  8:28780/1‑149 (MQ=255)                                     tCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑g                                                                                                             <  1:348786/149‑1 (MQ=255)                                     tCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑g                                                                                                             <  2:9877/149‑1 (MQ=255)                                                   ctgGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGCTTACACCGAAcg                                                                                               <  2:236127/149‑2 (MQ=255)                                                       tAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCGGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGc                                                                                           <  4:42858/149‑1 (MQ=255)                                                       tAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGc                                                                                           <  7:28780/149‑1 (MQ=255)                                                       tAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGc                                                                                           <  8:193678/149‑1 (MQ=255)                                                              tggcgtTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTAccc                                                                                     >  7:32995/1‑149 (MQ=255)                                                                  gtTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAcc                                                                                                                                     <  8:178399/98‑1 (MQ=255)                                                                  gtTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAcc                                                                                                                                     >  7:178399/1‑98 (MQ=255)                                                                  gtTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCg                                                                                 >  2:336812/1‑149 (MQ=255)                                                                  gtTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCg                                                                                 >  4:345507/1‑149 (MQ=255)                                                                    tgctgcTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGct                                                                               >  6:252798/1‑149 (MQ=255)                                                                          tgGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCaa                                                                         >  2:329074/1‑149 (MQ=255)                                                                            ggCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAAc                                                                        <  8:32995/148‑1 (MQ=255)                                                                              cGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGc                                                                     >  1:173480/1‑149 (MQ=255)                                                                               gAAACCATCGGCCCGGATCAGTGGCAGACCTACCTGCCAAACTGGGCCGCAGTACTGGAGCTGGTGATTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGcc                                                                    >  8:239591/1‑149 (MQ=255)                                                                                   ccATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCa                                                                >  8:288801/1‑149 (MQ=255)                                                                                     aTCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACg                                                              >  7:15340/1‑149 (MQ=255)                                                                                         gCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCa                                                          <  6:163280/149‑1 (MQ=255)                                                                                         gCCCGGATCAGAGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCa                                                          <  8:341652/149‑1 (MQ=255)                                                                                          cccGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAcc                                                                                                                                     <  6:357688/74‑1 (MQ=255)                                                                                          cccGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAcc                                                                                                                                     >  5:357688/1‑74 (MQ=255)                                                                                          cccGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTgg                                                                                            <  5:212446/114‑1 (MQ=255)                                                                                          cccGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTgg                                                                                            >  6:212446/1‑114 (MQ=255)                                                                                          cccGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAg                                                         >  7:229813/1‑149 (MQ=255)                                                                                              gATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATc                                                     <  2:173480/149‑1 (MQ=255)                                                                                                tCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCa                                                   <  5:252798/149‑1 (MQ=255)                                                                                                          acctacctGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTcc                                           <  8:15340/147‑1 (MQ=255)                                                                                                            ctacctGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGc                                       >  7:351579/1‑149 (MQ=255)                                                                                                               cctGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCCTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAg                                    >  4:90711/1‑149 (MQ=255)                                                                                                                 tGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGcc                                  >  4:103469/1‑149 (MQ=255)                                                                                                                  gCCAATCTGGGCCGCAGTACTGGTACTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTTTAAGCCGAACCTGGCGTACCCGTCGATATTCAACCGCCTCCACGGTCAGAAACAAGCCATACCAGGCAAGCCg                                 >  7:235557/1‑149 (MQ=255)                                                                                                                   ccAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCg                                 <  1:336812/148‑1 (MQ=255)                                                                                                                    cAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGt                               >  8:319242/1‑149 (MQ=255)                                                                                                                         tGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCAcc                           >  6:104756/1‑148 (MQ=255)                                                                                                                          gggCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGa                         >  2:122884/1‑149 (MQ=255)                                                                                                                                cAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCg                   >  4:138780/1‑149 (MQ=255)                                                                                                                                     ctggtgctggtgCTTGAGGGCTTCAAGCACCTGCTACGTGAGCGTCGCCGTTC‑GGTGTTAGCCGTACCTGGCGTACGCGTCGCTGTTCAACCTCCTCCACGGTGAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAg               <  7:239591/148‑1 (MQ=255)                                                                                                                                            tggtgctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCAACGACCATCGGCAGGCGGTCCa       >  7:58326/1‑149 (MQ=255)                                                                                                                                            aggtgctTGAAGCCCTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTc         >  7:102770/2‑147 (MQ=255)                                                                                                                                            aggtgctTGAAGCCCTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCACCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTc         >  6:51066/2‑147 (MQ=255)                                                                                                                                              gtgctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGc     <  8:206513/149‑1 (MQ=255)                                                                                                                                               tgctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCa    <  8:229813/149‑1 (MQ=255)                                                                                                                                                gctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCAACCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGc     >  1:258558/1‑148 (MQ=255)                                                                                                                                                gctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTc         <  3:24926/143‑1 (MQ=255)                                                                                                                                                gctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTc         >  4:24926/1‑143 (MQ=255)                                                                                                                                                 ctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAg  >  5:370167/1‑149 (MQ=255)                                                                                                                                                  ttGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAg  >  4:357285/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                 TGCTGCCATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGC‑GTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTC‑GGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAG  >  NC_002947/3893334‑3893625

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