Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,466,055	G→C	10.0%	noncoding (105/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,055	0	G	C	10.0%	176.6 / 8.5	66	noncoding (105/119 nt)	PP_mr39	P15 ncRNA
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base C (5/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGGAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCA  >  NC_002947/3465909‑3466193                                                                                                                                                   |                                                                                                                                               aTCGGAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCggg                                                                                                                                                >  7:52664/1‑147 (MQ=255) aTCGGAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCggg                                                                                                                                                >  5:56639/1‑147 (MQ=255)         cggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCg                                                                                                                                   <  2:288295/149‑1 (MQ=255)          ggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCgg                                                                                                                                  >  2:216275/1‑149 (MQ=255)          ggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCgg                                                                                                                                  <  5:65214/149‑1 (MQ=255)           gcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCATGAGTAGGTGGTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCTAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCggg                                                                                                                                 <  1:36957/149‑1 (MQ=255)           gcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGCTTCCC‑‑GA‑CCggg                                                                                                                                 <  8:85935/149‑1 (MQ=255)            cAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCggg                                                                                                                                                <  8:92139/136‑1 (MQ=255)            cAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCggg                                                                                                                                                >  7:92139/1‑136 (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                                                            <  8:63181/135‑1 (MQ=255)                                                                                             ccTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTgg                                                            >  7:63181/1‑135 (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)                                                                                                                                                   |                                                                                                                                               ATCGGAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCC‑‑GA‑CCGGGGGTTTGAGTATTCTCATGTTCGGCGG‑CCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCA  >  NC_002947/3465909‑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: ‑