Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,175,005	G→C	7.1%	W27S (TGG→TCG)	PP_1030 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,175,005	0	G	C	7.1%	175.1 / 4.0	56	W27S (TGG→TCG)	PP_1030	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/22);  new base C (2/2);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

GTCGCTACCTGAGCGCCTGAGCGGTAGAATTGCGTCATTACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCC  >  NC_002947/1174859‑1175148                                                                                                                                                     |                                                                                                                                                ttCGCTAGCTGAGCGCCTGAGCGGTAGATTTGCGTCATTACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTATACCGTCCCGGTGCTTGAAGCCCGTGCACG‑CGCAATCCTGCgc‑tgg                                                                                                                                                <  8:255507/147‑1 (MQ=255) gTCGCTACCTGAGCGCCTGAGCGGTAGAATTGCGTCATTACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCgc‑tggc                                                                                                                                               >  6:47640/1‑149 (MQ=255)  tCGCTACCTGAGCGCCTGAGCGGTAGAATTGCGTCATTACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCgc‑tggct                                                                                                                                              >  1:257181/1‑149 (MQ=255)  tCGCTACCTGAGCGCCTGAGCGGTAGAATTGCGTCATTACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCgc‑tggct                                                                                                                                              >  3:3876/1‑149 (MQ=255)  tCGCTACCTGAGCGCCTGAGCGGTAGAATTGCGTCATTACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCgc‑tggc                                                                                                                                               <  5:101476/148‑1 (MQ=255)        cgtgagcgcctgagggGTAGAATTGCGTCATTACGGGCGTTTCTCAGGCTGGTCCCGTCCATTACCACCCATAGCATCATCGTTCCAACGCTGGCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTgc                                                                                                                                         <  1:93119/146‑1 (MQ=255)         ctgagcgcctgagcgGTAGAATTGCGTCATTACCGGAGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTgcgc                                                                                                                                       <  2:154784/149‑1 (MQ=255)          tgagcgcctgagcgGTAGAATTGCGTCATTACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTgcgc                                                                                                                                       >  1:295203/1‑148 (MQ=255)              agcctgagcgGTAGAATTGCGTCATTACCGGCGTTTATCAGGCTTGACCCGTTCATGAACAACCAGAGCATCATCGTTCCAAAGCAGTCTACCGACACGGTGCATGAAGCCCGTGCCAG‑CGCCATCCTGCGC‑TGGATCGTGCGCGaa                                                                                                                                    >  6:90563/2‑147 (MQ=255)               gcctgagcgGTAGAATTGCGTCATTACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCgc‑tggc                                                                                                                                               <  2:290903/135‑1 (MQ=255)               gcctgagcgGTAGAATTGCGTCATTACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAg                                                                                                                                 >  4:173609/1‑149 (MQ=255)               gcctgagcgGTAGAATTGCGTCATTACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTATACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAg                                                                                                                                 <  1:339246/149‑1 (MQ=255)                 ctgagcgGTAGAATTGCGTCATTACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGgt                                                                                                                               >  1:12965/1‑149 (MQ=255)                       ggTAGAATTGCGTCATTACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGgtcgtc                                                                                                                           >  3:74126/1‑147 (MQ=255)                             aTTGCGTCATTACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGagc                                                                                                                    <  4:3876/148‑1 (MQ=255)                                  gTCATTACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGa                                                                                                              >  1:176212/1‑149 (MQ=255)                                  gTCATTACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGa                                                                                                              >  3:263436/1‑149 (MQ=255)                                  gTCATTACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGa                                                                                                              <  3:173609/149‑1 (MQ=255)                                      ttACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGaccac                                                                                                          >  1:284457/1‑149 (MQ=255)                                       tACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TCGCTGGTGCGCTAAAAGGTCGTCGAGGAGCAGCTGaccacc                                                                                                         >  6:163807/1‑149 (MQ=255)                                       tACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TCGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGaccacc                                                                                                         >  6:168624/1‑149 (MQ=255)                                         ccGGAGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCAGGTGCATGAAGCCCGTGCCAG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCt                                                                                                        <  5:90563/148‑1 (MQ=255)                                          cGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATAATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCtgt                                                                                                      >  7:300402/1‑149 (MQ=255)                                             cGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGc                                                                                                   <  5:47640/149‑1 (MQ=255)                                                   tcAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAACGGTCGTCGAGGAGCAGCTGACCCCCTGTGGCCGCAcc                                                                                             >  8:166802/1‑149 (MQ=255)                                                     aGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCgg                                                                                           >  6:39190/1‑149 (MQ=255)                                                        cTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTaa                                                                                        <  2:12965/149‑1 (MQ=255)                                                           gTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAAcc                                                                                      <  4:263436/148‑1 (MQ=255)                                                             cccGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTAACGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCTC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTCTGGCCGCACCGGTAACCGCa                                                                                   >  5:72492/1‑149 (MQ=255)                                                              ccGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCg                                                                                            <  4:270720/139‑1 (MQ=255)                                                              ccGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCg                                                                                            >  3:270720/1‑139 (MQ=255)                                                                gTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATgg                                                                                >  8:160537/1‑149 (MQ=255)                                                                 ttCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATggg                                                                               >  7:212917/1‑149 (MQ=255)                                                                    aTGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGc                                                                                                   >  2:269634/1‑126 (MQ=255)                                                                    aTGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGc                                                                                                   <  1:269634/126‑1 (MQ=255)                                                                     tGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCAt                                                                           >  7:119019/1‑149 (MQ=255)                                                                      gAACACCCAGAGCATCATCCTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGAGCCCG‑CGCCATCCTGCGC‑TGGCTCGTGAGCGAAAAGGTCGTCGAGGAGCAGCAGACAACCTCTGGCCGCACCGGTAACCGCATGGGCCATg                                                                          >  6:273067/1‑149 (MQ=255)                                                                       aaCACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGc                                                                         >  5:120190/1‑149 (MQ=255)                                                                       aaCACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TCGCTGGTGCGCGAAAAGGTCGTCGAGTAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATg                                                                          <  5:163807/148‑1 (MQ=255)                                                                       aaCACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TCGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATg                                                                          <  5:168624/148‑1 (MQ=255)                                                                            ccAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCG‑AAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCAt                                                                           >  8:219938/1‑141 (MQ=255)                                                                            ccAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCG‑AAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCAt                                                                           <  7:219938/141‑1 (MQ=255)                                                                                  catcatCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAg                                                              <  2:77242/149‑1 (MQ=255)                                                                                   atcatcGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAgg                                                             >  5:172457/1‑149 (MQ=255)                                                                                    tcatcGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAggg                                                            <  2:295203/149‑1 (MQ=255)                                                                                              aaaGCTGTCTACCGTCCCGGTGCATGAAGCCCGAGCCCC‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAg                                                  <  7:306624/149‑1 (MQ=255)                                                                                                     tCTACCGTCCCGGTGCATGAAGCCCGTGCCCGCCGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGccc                                            >  1:41352/1‑149 (MQ=255)                                                                                                     tCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGcc                                             <  8:300402/147‑1 (MQ=255)                                                                                                       tACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGCGTGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTg                                                >  1:65289/1‑143 (MQ=255)                                                                                                       tACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTg                                                <  2:65289/142‑1 (MQ=255)                                                                                                               cGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGa                                  >  5:213369/1‑148 (MQ=255)                                                                                                                ggTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGaa                                 >  2:45000/1‑148 (MQ=255)                                                                                                                 gTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGTGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGaaa                                >  4:317096/1‑148 (MQ=255)                                                                                                                  tGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGaaaaa                              <  2:41352/149‑1 (MQ=255)                                                                                                                       gAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGcc                         >  1:23997/1‑149 (MQ=255)                                                                                                                            ccGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCg                    <  7:160537/149‑1 (MQ=255)                                                                                                                                 cccG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGaa               >  7:91127/1‑149 (MQ=255)                                                                                                                                  ccG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAAc              <  6:120190/149‑1 (MQ=255)                                                                                                                                           tCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAAc      <  4:74126/149‑1 (MQ=255)                                                                                                                                            ccTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACg     >  2:291136/1‑149 (MQ=255)                                                                                                                                             cTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTTACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACg     <  6:72492/148‑1 (MQ=255)                                                                                                                                               gcgc‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGcc  >  1:52412/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                GTCGCTACCTGAGCGCCTGAGCGGTAGAATTGCGTCATTACCGGCGTTTCTCAGGCTTGTCCCGTTCATGAACACCCAGAGCATCATCGTTCCAAAGCTGTCTACCGTCCCGGTGCATGAAGCCCGTGCCCG‑CGCCATCCTGCGC‑TGGCTGGTGCGCGAAAAGGTCGTCGAGGAGCAGCTGACCACCTGTGGCCGCACCGGTAACCGCATGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCC  >  NC_002947/1174859‑1175148

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