Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,390,663	T→A	6.2%	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	6.2%	187.8 / 2.2	64	noncoding (31/87 nt)	PP_mr56	C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base T (27/33);  new base A (2/2);  total (29/35)
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.45e-01

CTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGACCCGCCTCACCACC  >  NC_002947/5390521‑5390802                                                                                                                                                |                                                                                                                                               ctgcctgcGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTtgcgt                                                                                                                                           <  5:238352/148‑1 (MQ=255)     ctgcGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCAAGCGTGcccc                                                                                                                                      <  2:23194/149‑1 (MQ=255)      tgcGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCg                                                                                                                                     >  5:137608/1‑149 (MQ=255)        cgcCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGCGCAC‑GCAAGCGTGccccatg                                                                                                                                   <  5:228260/149‑4 (MQ=255)         gcCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGtt                                                                                                                                  >  2:230607/1‑149 (MQ=255)              ttttCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCtgt                                                                                                                             <  1:230607/149‑1 (MQ=255)                 tCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGtt                                                                                                                          >  3:132591/1‑149 (MQ=255)                    gTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGc                                                                                                                       <  2:313070/149‑1 (MQ=255)                    gTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGc                                                                                                                       <  2:150163/149‑1 (MQ=255)                    gTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGc                                                                                                                       <  2:154988/149‑1 (MQ=255)                     ttGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATAAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGcc                                                                                                                      >  2:253302/1‑149 (MQ=255)                           tCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATAAACGACCGGCTTTGACAGGCCGCAACGATACACAGAGCTCATAACCGTTCATGGCGGTTGTGAATGGGGCAC‑GCTTCAGTGCTCCTGTTTATGTTTTAGCCGGTCTg                                                                                                                >  5:94937/1‑149 (MQ=255)                                     aGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCAt                                                                                                      <  6:137608/148‑1 (MQ=255)                                        gTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGc‑ac                                                                                                 <  3:9360/149‑1 (MQ=255)                                        gTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGc‑a                                                                                                  >  7:85324/1‑148 (MQ=255)                                           ttAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAg                                                                                               <  8:273684/148‑1 (MQ=255)                                            tAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGc                                                                                              <  4:132591/148‑1 (MQ=255)                                              gTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGc                                                                                           >  8:274019/1‑149 (MQ=255)                                                 tcGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGcca‑c                                                                                       >  2:138243/1‑149 (MQ=255)                                                   ggTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGt                                                                                                                   >  2:12753/1‑122 (MQ=255)                                                   ggTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGt                                                                                                                   <  1:12753/122‑1 (MQ=255)                                                        gTCGTTGATGAGCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAAc                                                                                                          >  7:189301/1‑125 (MQ=255)                                                        gTCGTTGATGAGCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAAc                                                                                                          <  8:189301/125‑1 (MQ=255)                                                         tCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGt                                                                               <  3:102426/149‑1 (MQ=255)                                                         gcGTTGCTCATCAACGACCGGCTTTGACAGGCCGCCACGATACACAGCGCCCCTAACCTTTCATGGCTGTTGTGCATGTGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGt                                                                               <  6:94937/148‑1 (MQ=255)                                                          cGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGtt                                                                              <  4:201513/149‑1 (MQ=255)                                                           gTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGttt                                                                             >  4:43178/1‑149 (MQ=255)                                                           gTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGttt                                                                             <  1:138243/149‑1 (MQ=255)                                                                  atcaACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGt                                                                                                               >  2:234357/1‑111 (MQ=255)                                                                  atcaACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGt                                                                                                               <  1:234357/111‑1 (MQ=255)                                                                  atcaACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCg                                                                      >  5:128660/1‑149 (MQ=255)                                                                     aaCGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGA‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGt                                                                   >  1:58069/1‑149 (MQ=255)                                                                           cGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGTCCAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCg                                                              <  1:225077/149‑1 (MQ=255)                                                                                 tGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGTCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCAGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGc                                                         <  1:246696/147‑1 (MQ=255)                                                                                 tGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGc                                                         <  2:120752/147‑1 (MQ=255)                                                                                  gACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGccc                                                       <  8:85324/148‑1 (MQ=255)                                                                                   aCAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGcccc                                                      >  5:124524/1‑148 (MQ=255)                                                                                        ccGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCga                                                <  2:105087/149‑1 (MQ=255)                                                                                         cGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCgag                                               <  7:274019/149‑1 (MQ=255)                                                                                          gcaaataTACACCGCGCGCATAACCGTTCAAGGCGGTTGTGCATGTGGCAC‑GCTTGCGAGCCCAGGTTTCACAGATCCCCGCACTGT‑CAACCCAAGC‑AAAGCTGCCA‑CCAATCCGTTTGACAGCGGGTTGCCGGTGGCACCTAACgag                                               >  2:253380/7‑148 (MQ=255)                                                                                          gCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCgg                                                                     >  1:281653/1‑126 (MQ=255)                                                                                          gCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCgg                                                                     <  2:281653/126‑1 (MQ=255)                                                                                          gCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCgag                                               >  5:297358/1‑148 (MQ=255)                                                                                             aCGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAaca                                           <  6:297358/149‑1 (MQ=255)                                                                                               gATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTgg                                                          >  8:123284/1‑132 (MQ=255)                                                                                               gATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTgg                                                          <  7:123284/132‑1 (MQ=255)                                                                                               gATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAaca                                           >  5:177461/1‑147 (MQ=255)                                                                                                    acaGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTaa                                     >  4:55648/1‑148 (MQ=255)                                                                                                    acaGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTaa                                     >  3:84139/1‑148 (MQ=255)                                                                                                      aGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAAt                                  <  3:55648/149‑1 (MQ=255)                                                                                                        cgcgcATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGCTTCAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATg                                 <  1:253380/149‑1 (MQ=255)                                                                                                         gcgcATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGc                                >  1:312238/1‑148 (MQ=255)                                                                                                              tAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGa                          <  6:128660/149‑1 (MQ=255)                                                                                                                 ccGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCg                       >  6:196631/1‑149 (MQ=255)                                                                                                                   gTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGtt                     <  5:293994/149‑1 (MQ=255)                                                                                                                   gTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGtt                     >  8:38904/1‑149 (MQ=255)                                                                                                                    ttCTTGGGGTTTGTGAAT‑GGGCACGGATTGCTTGCGCCGGTTACTGTGTTCGCGGGCGTGT‑CACCCATTGA‑ACAGCTGCAATACATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGtt                     <  1:253302/149‑1 (MQ=255)                                                                                                                      cATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTcc                   <  6:177461/148‑1 (MQ=255)                                                                                                                          gCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGAc               <  3:43178/148‑1 (MQ=255)                                                                                                                                   cATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGGCCCGGCTcgcg     >  2:262541/1‑146 (MQ=255)                                                                                                                                   accggcgCAC‑GCAAGCGTGCCCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGACCCGCCTcac      >  6:209822/7‑148 (MQ=255)                                                                                                                                     tGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGACCCGCCTcaccac   <  6:124524/149‑1 (MQ=255)                                                                                                                                      ggggCAC‑GATTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAAACCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACAAAGTAAAATGCTAAAGATCTTTCCTGACCCGCCTCcccacc  >  8:101044/1‑149 (MQ=255)                                                                                                                                      ggcgCAC‑GCAAGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGACCCGCCTCaccacc  >  1:293601/4‑149 (MQ=255)                                                                                                                                                |                                                                                                                                               CTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTTGCTCATCAACGACCGGCTTTGACAGGCCGCAACGATACACAGCGCGCATAACCGTTCATGGCGGTTGTGCATGGGGCAC‑GCTTGCGTGCGCCGGTTTCTGTGTTCGCCGGTCTGT‑CAACCCATGC‑ACAGCTGCCA‑CCATTCCGTTTGACAGCGGGTTGGCGGTGGCCCCTATCGAGAACACAGTAAAATGCTAAAGATCGTTCCTGACCCGCCTCACCACC  >  NC_002947/5390521‑5390802

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