Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,390,663	T→A	8.7%	noncoding (31/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,663	0	T	A	8.7%	185.5 / 7.7	69	noncoding (31/87 nt)	PP_mr56	C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base T (32/31);  new base A (3/3);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

TACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGACCCGCCTCACCACCTCCACTCCC  >  NC_002947/5390519‑5390811                                                                                                                                                  |                                                                                                                                                     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                                                                                                                                             <  6:133336/144‑3 (MQ=255)           gcCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCAagcgtgccc                                                                                                                                             >  5:133336/1‑142 (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:70907/1‑149 (MQ=255)                                                                  tcatcaACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAg                                                                                 >  4:98170/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                                                                                      >  1:70089/1‑106 (MQ=255)                                                                                                        aGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGttt                                                                                      <  2:70089/106‑1 (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)                                                                                                                                              aCGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGACCCGCCTCACCACCTCCAc      >  5:350255/1‑149 (MQ=255)                                                                                                                                                gCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGACCCGCCTCACCACCTCCACTc    >  4:15111/1‑149 (MQ=255)                                                                                                                                                  ttGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTTCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGACCCGCCTCACCACCTCCACTccc  >  7:213207/1‑149 (MQ=255)                                                                                                                                                  ttGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTc                             <  7:254124/122‑1 (MQ=255)                                                                                                                                                  ttGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTc                             >  8:254124/1‑122 (MQ=255)                                                                                                                                                  |                                                                                                                                                     TACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGC‑GGTTGTGCATGGGGCACGCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCAACCCATGCACAGCTGCCACCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGACCCGCCTCACCACCTCCACTCCC  >  NC_002947/5390519‑5390811

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