Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	713,884	A→G	6.5%	D117G (GAC→GGC)	fkpB →	peptidylprolyl isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	713,884	0	A	G	6.5%	182.6 / 2.5	63	D117G (GAC→GGC)	fkpB	peptidylprolyl isomerase
Reads supporting (aligned to +/- strand):  ref base A (24/34);  new base G (2/2);  total (26/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

AGAACGCGTTCGGCCAGCCTAACCCGCAGAACGTGCAGGTCATGCCGAGGTCCAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGCCGAGCATGCAAATCAAACTCGCCAACCCTCGCGGCTTCTGCGCCGGGGT  >  NC_002947/713737‑714030                                                                                                                                                    |                                                                                                                                                   agaACGCGTTCGGCCAGCCTAACCCGCAGAACGTGCAGGTCATGCCGAGGTCCAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATgac                                                                                                                                                   <  2:256758/149‑1 (MQ=255)  gaACGCGTTCGGCCAGCCTAACCCGCAGAACGTGCAGGTCATGCCGAGGTCCAACTACGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCACCACCGAGCTGCCGGGTGAGGTCAAAGCGTTCGATgacg                                                                                                                                                  >  5:210625/1‑149 (MQ=255)     cgcgTTCGGCCAGCCTAACCCGCAGAACGTGCAGGTCATGCCGAGGTCCAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATgacgac                                                                                                                                                >  5:149763/1‑148 (MQ=255)      gcgTTCGGCCAGCCTAACCCGCAGAACGTGCAGGTCATGCCGAGGTCCAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCACCACCGAGCTGCCGGGTGTGGTCAACGCGTTCGATGACGACCa                                                                                                                                              <  5:282864/149‑1 (MQ=255)       cgTTCGGCCAGCCTAACCCGCAGAACGTGCAGGTCATGCCGAGGTCCAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAg                                                                                                                                             <  2:2710/149‑1 (MQ=255)                  ccTAACCCGCAGAACGTGCAGGTCATGCCGAGGTCCAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGa                                                                                                                                  <  5:349797/149‑1 (MQ=255)                  ccTAACCCGCAGAACGTGCAGGTCATGCCGAGGTCCAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGa                                                                                                                                  >  3:76680/1‑149 (MQ=255)                   cTAACCCGCAGAACGTGCAGGTCATGCCGAGGTCCAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGAc                                                                                                                                 <  2:20752/149‑1 (MQ=255)                      aCCCGCAGAACGTGCAGGTCATGCCGAGTTCCAACTTCGAAGGCATTGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACtt                                                                                                                               <  5:81385/148‑1 (MQ=255)                       cccGCAGAACGTGCAGGTCATGCCGAGGTCAAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTc                                                                                                                              >  8:316814/1‑148 (MQ=255)                        ccGCAGAACGTGCAGGTCATGCCGAGGTCCAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCaa                                                                                                                            >  6:24576/1‑149 (MQ=255)                           cAGAACGTGCAGGTCATGCCGAGGTCCAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGGCGACCAGGTCACCATCGACTTCAATCa                                                                                                                         <  3:218822/149‑1 (MQ=255)                            agcaCGTGCAGGTCATGCCGAGGGCCAACTTCGAAGGCATGGAGTTGGCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCAc                                                                                                                        <  6:210625/146‑1 (MQ=255)                            aGAACGTGCAGGTCATGCCGAGGTCCAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGGGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCAc                                                                                                                        <  2:76302/149‑1 (MQ=255)                                      ggTCATGCCGAGGTCCAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCAGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCActggctg                                                                                                              >  2:91267/1‑149 (MQ=255)                                             ccGAGGTCCAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCGAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGAc                                                                                                       <  7:321764/149‑1 (MQ=255)                                             ccGAGGTCCAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGa                                                                                                        <  5:24576/148‑1 (MQ=255)                                               gAGGTCCAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGAccc                                                                                                     >  2:228170/1‑149 (MQ=255)                                                    ccAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGAcc                                                                                                >  3:334960/1‑149 (MQ=255)                                                     cAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTc                                                                                                                                          >  2:73545/1‑106 (MQ=255)                                                     cAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCGCCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCt                                                                                               >  7:209245/1‑149 (MQ=255)                                                        cTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTc                                                                                             >  6:115770/1‑148 (MQ=255)                                                            gAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGt                                                                                        >  8:120499/1‑149 (MQ=255)                                                             aaGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTg                                                                                       <  6:6082/149‑1 (MQ=255)                                                              aGGCTTGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCCCAAACACCGCGCTGCCGGGTGTGGTCAAAGCGTTCTATGTCGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTgg                                                                                      <  8:58255/149‑1 (MQ=255)                                                              aGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTgg                                                                                      <  3:176026/149‑1 (MQ=255)                                                                  aTGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTg                                                                                       >  4:9757/1‑144 (MQ=255)                                                                  aTGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGATg                                                                                       >  3:198060/1‑144 (MQ=255)                                                                      aGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTg                                                                                                   >  5:328684/1‑128 (MQ=255)                                                                      aGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTg                                                                                                   <  6:328684/128‑1 (MQ=255)                                                                        ttGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGa                                                                            <  6:149763/149‑1 (MQ=255)                                                                        ttGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGa                                                                            <  4:76680/149‑1 (MQ=255)                                                                        ttGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGa                                                                            <  8:269070/149‑1 (MQ=255)                                                                           tCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAgg                                                                          <  7:120499/148‑1 (MQ=255)                                                                           tCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGt                                                                         <  2:330475/149‑1 (MQ=255)                                                                            ccGACGGGCTGCTGATCATCTTCAACGCCGCCGCCCACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTg                                                                        <  4:198060/149‑1 (MQ=255)                                                                               aaGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAg                                                                     <  4:334960/149‑1 (MQ=255)                                                                                  ggCTGCTGATCATCTTCAACTACGCCGCCAACACCGAGCTCCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATAGACTTCAAACACCCACTGGCTGGCAAGACCCTGACATTCGAGTTGGAGATCCTCGAGGTGAAGGcc                                                                  >  8:35211/1‑149 (MQ=255)                                                                                       ctgATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGa                                                                            >  7:90503/1‑134 (MQ=255)                                                                                       ctgATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGa                                                                            <  8:90503/134‑1 (MQ=255)                                                                                       ctgATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTC‑AAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTaa                                                            <  7:247657/149‑1 (MQ=255)                                                                                           tcatcTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGcc                                                         <  7:316814/149‑1 (MQ=255)                                                                                             atcTTCAACGACGCCGCCAACACCGCGCTGCCGGGTGTGGTCAACGCTTTCGATGACGACCAGGGCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCt                                                        <  6:6615/148‑1 (MQ=255)                                                                                             atcTTCAACGACGCCGCCAACACCGAGCTGCCAGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCTAGACCCTGAGCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCt                                                        <  1:91267/148‑1 (MQ=255)                                                                                               cTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTg                                                       >  6:82790/1‑147 (MQ=255)                                                                                                  cAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGa                                                     <  7:39501/146‑1 (MQ=255)                                                                                                  cAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGcc                                                  >  7:122499/1‑149 (MQ=255)                                                                                                  cAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGc                                                   <  4:86474/148‑1 (MQ=255)                                                                                                    acgacgCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGCCGa                                                <  5:82790/149‑1 (MQ=255)                                                                                                     cgacgCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGCCGAg                                               >  7:140684/1‑149 (MQ=255)                                                                                                         gccgccAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGCCGAGCATg                                           <  8:140684/149‑1 (MQ=255)                                                                                                               aaCACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGcc                                                                  >  7:20922/1‑120 (MQ=255)                                                                                                               aaCACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGcc                                                                  <  8:20922/120‑1 (MQ=255)                                                                                                                  acCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGCCGAGCATGCAAATCaaa                                  <  3:9757/149‑1 (MQ=255)                                                                                                                   ccgttcTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGCCGAGCATGCAAATCAAAc                                 <  3:260382/144‑1 (MQ=255)                                                                                                                           ccGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGCCGAGCATGCAAATCAAACTCGCCAAc                         >  4:112368/1‑149 (MQ=255)                                                                                                                               gtgtGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGCCGAGCATGCAAATCaaa                                  <  3:282531/136‑1 (MQ=255)                                                                                                                               gtgtGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGCCGAGCATGCAAATCaaa                                  >  4:282531/1‑136 (MQ=255)                                                                                                                               gtgtGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGCCGAGCATGCAAATCAAACTCGCCAACCCTc                     <  1:228170/149‑1 (MQ=255)                                                                                                                                  tgGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGCCGAGCATGCAAATCAAACTCGCCAACCCTcgcg                  >  5:310652/1‑149 (MQ=255)                                                                                                                                             gtCGATGGTGACCTGGTCGTCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGCCGAGCATGCAAATCAAACTCGCCAACCCTCGCGGCTTc             <  4:281396/142‑1 (MQ=255)                                                                                                                                             gtCGATGGTGACCTGGTCGTCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGCCGAGCATGCAAATCAAACTCGCCAACCCTCGCGGCTTc             >  3:281396/2‑143 (MQ=255)                                                                                                                                             gtCGATGGTGACCTGGTCGTCATCGACTTCAATCACCCACTGGCTGGCAACACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGCCGAGCATGCAAATCAAACTCGCCAACCCTCGCGGCTTCTGAGcc       >  5:55569/2‑149 (MQ=255)                                                                                                                                                gATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGCCGAGCATGCAAATCAAACTCGCCAACCCTCGCGGCTTCTGCGCCggg    >  7:226993/1‑149 (MQ=255)                                                                                                                                                gATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGCCGAGCATGCAAATCAAACTCGCCAACCCTCGCGGCTTCTGCGCCggg    <  3:112368/149‑1 (MQ=255)                                                                                                                                                  tGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGCCGAGCATGCAAATCAAACTCGCCAACCCTCGCGGCTTCTGCGCCGGGGt  >  3:303601/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   AGAACGCGTTCGGCCAGCCTAACCCGCAGAACGTGCAGGTCATGCCGAGGTCCAACTTCGAAGGCATGGAGTTGTCCGAAGGGCTGCTGATCATCTTCAACGACGCCGCCAACACCGAGCTGCCGGGTGTGGTCAAAGCGTTCGATGACGACCAGGTCACCATCGACTTCAATCACCCACTGGCTGGCAAGACCCTGACCTTCGAGGTGGAGATCCTCGAGGTGAAGGCCCTGTAAGCCTGGAGCCGAGCATGCAAATCAAACTCGCCAACCCTCGCGGCTTCTGCGCCGGGGT  >  NC_002947/713737‑714030

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