Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,466,057	A→T	9.1%	noncoding (103/119 nt)	PP_mr39 ←	P15 ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,466,057	0	A	T	9.1%	177.2 / 7.9	66	noncoding (103/119 nt)	PP_mr39	P15 ncRNA
Reads supporting (aligned to +/- strand):  ref base A (30/30);  new base T (4/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

CGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCA  >  NC_002947/3465917‑3466193                                                                                                                                             |                                                                                                                                             cggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCg                                                                                                                                   <  2:288295/149‑1 (MQ=255)  ggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCgg                                                                                                                                  <  5:65214/149‑1 (MQ=255)  ggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCgg                                                                                                                                  >  2:216275/1‑149 (MQ=255)   gcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCATGAGTAGGTGGTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCTAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCggg                                                                                                                                 <  1:36957/149‑1 (MQ=255)   gcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGCTTCCC‑‑GA‑CCggg                                                                                                                                 <  8:85935/149‑1 (MQ=255)    cAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCgggg                                                                                                                                <  1:214788/149‑1 (MQ=255)     aGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCggggg                                                                                                                               <  1:203336/149‑1 (MQ=255)      ggCTGGCCCGAACGGGAATGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGt                                                                                                                              <  7:47481/149‑1 (MQ=255)       gCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGCTTCCC‑‑GA‑CCGGGGGtt                                                                                                                             >  6:99900/1‑149 (MQ=255)           gCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTTAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGCTTCCC‑‑GA‑CCGGGGGttttg                                                                                                                          <  2:369033/148‑3 (MQ=255)             ccGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTa                                                                                                                       <  7:221767/149‑1 (MQ=255)               gAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTATTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTAtt                                                                                                                     <  8:11738/149‑1 (MQ=255)                aaCGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTAtt                                                                                                                     >  2:281137/1‑148 (MQ=255)                aaCGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATtc                                                                                                                    >  1:150455/1‑149 (MQ=255)                aaCGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATtc                                                                                                                    <  2:150455/149‑1 (MQ=255)                           cTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTcggc                                                                                                         <  8:52664/149‑1 (MQ=255)                            tAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTcggcg                                                                                                        >  5:341525/1‑149 (MQ=255)                            tAAATCGCCAGGCGTAGGAGTTTGCGTAGTAATAAGTGGCGTAGATAATGAATCTGTCCTCAAAGTGAATCGGTGATCGTCAAAGGCCTGAAAAACAAAACCCCCGGACGGCAAGCC‑‑GA‑CCGTGGGTTTGAGTATTCTCATGTTcggcg                                                                                                        >  2:324070/1‑149 (MQ=255)                             aaaTCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTcggcgg                                                                                                       >  2:341525/1‑149 (MQ=255)                                       gcgtttgtgtTTGCGTTGTAATTTGTGGGGTTGTTCATTTATCTGTTCTCAAAGTGATTCGGTGTTGGTCAAATGTCTGAAACACAATACGCCCTGTCGTGAAGCC‑‑TA‑CGGGGGGTTTGAGTATTCTCATGTTCTGCGG‑CCCCTCGaa                                                                                             <  8:137833/142‑1 (MQ=255)                                        ggttggTTTTTGTGTTTTATTAAGTGGTGTATTTCATTATCTTGTCCTCAAAGTGAATCTGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTGGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAgg                                                                                           <  1:324070/145‑1 (MQ=255)                                          tAGGTTTTTGCGTAGTAATAAGTGGCGTAGTTCATTAATCTGTCCTCAAAGTGAATCGGTGGTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTTGAAGGTg                                                                                         <  7:308928/149‑1 (MQ=255)                                             gtgtTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGAc                                                                                      >  3:64291/1‑149 (MQ=255)                                                   gCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCgg                                                                                 <  6:56639/148‑1 (MQ=255)                                                             aaGTGGCGTAGTTCATGATTCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTAtt                                                                                                                     <  8:137548/103‑1 (MQ=255)                                                             aaGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTAtt                                                                                                                     >  7:137548/1‑103 (MQ=255)                                                              aGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGgc                                                                     <  6:341525/149‑1 (MQ=255)                                                                  gCGTAGTTTATGGATATGTACTAAAAGTGAATCGGTGTTGGTCAAAGGCCAGAAAAAAAAAAACCCCCGTGGGGGAGCC‑‑GA‑CCGGGGGTTTTAGTATTATCATGTTAGGCGG‑GCCCTCGAAGGTGGGCGCCGGGTGGGTATCGGCGcc                                                                  >  2:201035/1‑148 (MQ=255)                                                                   cGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGcc                                                                  >  6:240891/1‑147 (MQ=255)                                                                    gTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAg                                                                                            >  7:351200/1‑120 (MQ=255)                                                                    gTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAg                                                                                            <  8:351200/120‑1 (MQ=255)                                                                     tAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGt                                                              >  5:265217/1‑149 (MQ=255)                                                                           aTGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTa                                                         >  8:252342/1‑148 (MQ=255)                                                                             gAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCACGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAAc                                                      >  5:92852/1‑149 (MQ=255)                                                                              aaTCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCAGGAGCCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTa                                                         >  1:315536/1‑148 (MQ=255)                                                                              aaTCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTaaa                                                       <  1:281137/148‑1 (MQ=255)                                                                                tCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTCGTGGGTAAGCCac                                                   >  8:293017/1‑148 (MQ=255)                                                                                tCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCaa                                                   >  5:276683/1‑149 (MQ=255)                                                                                    tCCTCAAAGTGAATCGGTGTTGGTCAAAGTCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATAc                                                <  1:232966/148‑1 (MQ=255)                                                                                    tCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATAcc                                               >  5:363168/1‑149 (MQ=255)                                                                                     ccTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTgg                                                            >  7:63181/1‑135 (MQ=255)                                                                                     ccTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTgg                                                            <  8:63181/135‑1 (MQ=255)                                                                                      cTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCgggtgg                                                                             >  6:202951/1‑117 (MQ=255)                                                                                      cTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCgggtgg                                                                             <  5:202951/117‑1 (MQ=255)                                                                                         aaaGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAg                                          <  6:363168/149‑1 (MQ=255)                                                                                              gAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTaa                                       >  6:366888/1‑147 (MQ=255)                                                                                                 tCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCt                                  >  5:208861/1‑149 (MQ=255)                                                                                                 tCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCt                                  <  1:201035/149‑1 (MQ=255)                                                                                                     tgtTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCa                              <  1:159376/149‑1 (MQ=255)                                                                                                            cAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGCGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCa                       >  7:287012/1‑149 (MQ=255)                                                                                                             aaaGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCAc                      <  7:252342/149‑1 (MQ=255)                                                                                                               aGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACg                     >  5:136251/1‑148 (MQ=255)                                                                                                                 accTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAg                  >  5:60754/2‑149 (MQ=255)                                                                                                                  ccTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGc                 >  4:40890/1‑149 (MQ=255)                                                                                                                     gAAAAACAAAACCCCCAGTCGGGAAGCT‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGTCACGCAGCCgg              >  1:134519/1‑149 (MQ=255)                                                                                                                      aaaaaCAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGc             <  1:341525/149‑1 (MQ=255)                                                                                                                       aaaacaaaacCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGc             <  6:208861/148‑1 (MQ=255)                                                                                                                              aacCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCAGCAACGGCCa     <  6:265217/149‑1 (MQ=255)                                                                                                                              aacCCCCGGTCGGCTTCCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCa     >  6:228045/1‑149 (MQ=255)                                                                                                                               acCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAg    >  5:45210/1‑149 (MQ=255)                                                                                                                                cccccGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAg    >  7:353125/1‑148 (MQ=255)                                                                                                                                cccccGGTCGGCTTCCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATAACCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGc   >  7:213460/1‑149 (MQ=255)                                                                                                                                 ccccGGTCGGGAAGCC‑‑GT‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCa                       <  6:109872/128‑1 (MQ=255)                                                                                                                                 ccccGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGca  >  8:196166/1‑149 (MQ=255)                                                                                                                                 ccccGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGca  <  5:366888/149‑1 (MQ=255)                                                                                                                                  cccGGTCGGCTTCCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGca  >  8:250255/1‑148 (MQ=255)                                                                                                                                             |                                                                                                                                             CGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCA  >  NC_002947/3465917‑3466193

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑