Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,390,656	G→C	5.9%	noncoding (24/87 nt)	PP_mr56 →	C4 ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,390,656	0	G	C	5.9%	210.6 / 3.6	68	noncoding (24/87 nt)	PP_mr56	C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base C (2/2);  total (35/33)
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.58e-01

TCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGACCCGCCTCACCACCTC  >  NC_002947/5390510‑5390804                                                                                                                                                    |                                                                                                                                                     tCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATgggg                                                                                                                                                     >  1:329833/1‑148 (MQ=255)   ccAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCat                                                                                                                                                  >  6:112639/1‑148 (MQ=255)   ccAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCa                                                                                                                                                   <  1:308435/148‑1 (MQ=255)    cAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCAc                                                                                                                                                  >  3:210003/1‑148 (MQ=255)          tACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTtgc                                                                                                                                            <  1:290767/148‑1 (MQ=255)               cctgcGCCTGTTTTCCGTTGCTGTCATTTTTTCAGGGTTTTAGTCTCGGTGTTTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTgcgc                                                                                                                                      <  1:125607/149‑1 (MQ=255)               cctgcGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGCGCACGCAagcgtgccc                                                                                                                                      >  7:230306/1‑147 (MQ=255)               cctgcGCCTG‑TTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTAATGGC‑GGTTTTGCATGGGGCACGCTTGCGTGCGcc                                                                                                                                     >  5:302814/1‑149 (MQ=255)                   cgcCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGCGCACGCAAGCGTGccccat                                                                                                                                   <  8:230306/148‑3 (MQ=255)                    gcCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGtt                                                                                                                                 >  5:362192/1‑149 (MQ=255)                    gcCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCAagcgtgccc                                                                                                                                      >  5:133336/1‑142 (MQ=255)                    gcCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCAagcgtgccc                                                                                                                                      <  6:133336/144‑3 (MQ=255)                       tGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGCGCACGCAAGCGTGccccat                                                                                                                                   <  4:115217/144‑3 (MQ=255)                       tGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGCGCACGCAAGCGTGccccat                                                                                                                                   >  3:115217/1‑142 (MQ=255)                         ttttCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCt                                                                                                                              <  7:75226/147‑1 (MQ=255)                              cGTTGCGGTCTTTTTTTCACGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGCCCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGA‑GGTTGTGCATGGTGCACGCTTGCGTGCGACGGTTTCTGTGTTCg                                                                                                                       >  8:223958/1‑149 (MQ=255)                               gTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGCACGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGc                                                                                                                      >  8:259873/1‑149 (MQ=255)                                  gCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCg                                                                                                                    <  2:329833/148‑1 (MQ=255)                                   cGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGt                                                                                                                  <  4:210003/149‑1 (MQ=255)                                     gTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCtgt                                                                                                                            >  5:232065/1‑137 (MQ=255)                                     gTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCtgt                                                                                                                            <  6:232065/137‑1 (MQ=255)                                     gTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCt                                                                                                                <  5:347310/149‑1 (MQ=255)                                      tCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTg                                                                                                               <  3:308110/149‑1 (MQ=255)                                         ttttttCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTc                                                                                                             >  8:231396/1‑148 (MQ=255)                                               cAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCAt                                                                                                      >  2:149027/1‑149 (MQ=255)                                                 gggTTTTAGTCTCGGTTGGTCGTTGCTCATCAAGGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGc                                                                                                    <  7:133942/149‑1 (MQ=255)                                                     tttAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAg                                                                                                >  7:202443/1‑149 (MQ=255)                                                      ttAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCa                                                                                                       <  2:259143/141‑1 (MQ=255)                                                      ttAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCa                                                                                                       >  1:259143/1‑141 (MQ=255)                                                        aGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCTGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTg                                                                                             >  2:8477/1‑149 (MQ=255)                                                         gTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGc                                                                                            <  1:219187/149‑1 (MQ=255)                                                          tctcGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGc                                                                                            >  3:308117/1‑148 (MQ=255)                                                              ggTGGGTCGTTGCTCATAAAAGAACGGATTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGGATGGGGCACGCTTGCGTGCGACGGTTTCTGTGTTAGCCGGTCTGTCAACCCATGCACAGCTGacacca                                                                                       >  2:64694/1‑149 (MQ=255)                                                                 cggTCGTTGATGAGCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTc                                                                                    >  6:191929/2‑149 (MQ=255)                                                                      gTTTCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTTCGTGCGCCGGTTTCTGTGTTCTCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGttt                                                                               <  4:173143/149‑1 (MQ=255)                                                                      gTTGCTCATCAACGACCGGCTTTGGCAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCTACCCATGCACAGCTGCCACCATTCCGttt                                                                               <  7:157567/149‑1 (MQ=255)                                                                           tcatcaACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAg                                                                          >  4:98170/1‑149 (MQ=255)                                                                           tcatcaACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAg                                                                          >  4:70907/1‑149 (MQ=255)                                                                              tcaACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCg                                                                        >  1:309839/1‑149 (MQ=255)                                                                                aaCGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCggg                                                                      <  1:215097/148‑1 (MQ=255)                                                                                aaCGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGt                                                                     >  4:172974/1‑149 (MQ=255)                                                                                   gtccGTCTTTGACAGGCAGCAACGATACACAGCGCGAATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGGGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGATTCCACCATTCCGTTTGACAGCTGGTTgg                                                                  <  7:223958/147‑1 (MQ=255)                                                                                      cGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCgg                                                               >  1:370406/1‑149 (MQ=255)                                                                                        gCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTg                                                             >  6:140805/1‑149 (MQ=255)                                                                                             gACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGccc                                                         <  8:202443/148‑1 (MQ=255)                                                                                             gACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGccc                                                         <  1:149027/148‑1 (MQ=255)                                                                                              aCAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGcccc                                                        >  7:356791/1‑148 (MQ=255)                                                                                              aCAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCt                                                       <  6:362192/149‑1 (MQ=255)                                                                                                 ggCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATc                                                    <  4:308117/149‑1 (MQ=255)                                                                                                    cGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCgag                                                 <  5:112235/149‑1 (MQ=255)                                                                                                      cAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGaa                                               <  3:155659/149‑1 (MQ=255)                                                                                                          gATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAAcaca                                           <  3:172974/149‑1 (MQ=255)                                                                                                              cacaGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTa                                        <  7:231396/148‑1 (MQ=255)                                                                                                               acaGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTaa                                       >  1:164088/1‑148 (MQ=255)                                                                                                                caGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTaaaa                                     <  5:191929/149‑1 (MQ=255)                                                                                                                 aGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGttt                                                                               <  2:70089/106‑1 (MQ=255)                                                                                                                 aGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGttt                                                                               >  1:70089/1‑106 (MQ=255)                                                                                                                          aaCCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTa                                >  6:234028/1‑144 (MQ=255)                                                                                                                          aaCCGTTCATGGC‑GGTAGTGCATGGGGCACGCTTGCGAGCGCCGGTTTCTGGGTACGCCGGACGGTCAACCCATGCACCGCTGCCACCAGTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTAGCGAGAACGCAGTAAAATGCTa                                >  6:257458/1‑144 (MQ=255)                                                                                                                           aCCGTTCATGGC‑GGTTGTGCATGGGGCATGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATc                          <  5:112639/149‑1 (MQ=255)                                                                                                                             cGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGt                        <  8:356791/149‑1 (MQ=255)                                                                                                                              gTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGtt                       >  1:184517/1‑149 (MQ=255)                                                                                                                               ttCATGGC‑GGTTGTGGATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCATGTGTGTCAACCCATGCACATCTGCCACCATTCCGTTCGACAGCGGGTTGGCGGTGGCTCCTATCTAGAACACAGTAAAATGCTAAAGATCGTTc                      <  1:64694/149‑1 (MQ=255)                                                                                                                               ttCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTc                      >  8:112690/1‑149 (MQ=255)                                                                                                                                   tGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGa                  <  7:112690/149‑1 (MQ=255)                                                                                                                                          ttGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGACCCGcc            >  5:351901/1‑149 (MQ=255)                                                                                                                                               cATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGACCCGCCTcac        >  5:86802/1‑148 (MQ=255)                                                                                                                                                aTGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGACCCGCCTcacca      >  6:246611/1‑149 (MQ=255)                                                                                                                                                 tGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGACCCGCCTcaccac     >  6:354805/1‑149 (MQ=255)                                                                                                                                                    ggCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTACAGATCGTTCCTGACCCGCCTCACCACCTc  <  2:370406/149‑1 (MQ=255)                                                                                                                                                    ggCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGACCCGCCTCACCACCTc  <  6:302814/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGACCCGCCTCACCACCTC  >  NC_002947/5390510‑5390804

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