Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,466,056	A→T	6.6%	noncoding (104/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,056	0	A	T	6.6%	230.7 / 4.4	76	noncoding (104/119 nt)	PP_mr39	P15 ncRNA
Reads supporting (aligned to +/- strand):  ref base A (41/30);  new base T (2/3);  total (43/33)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

GATCGGAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAG  >  NC_002947/3465908‑3466203                                                                                                                                                       |                                                                                                                                                      gATCGGAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGa                                                                                                                                                       <  2:107108/149‑1 (MQ=255)      gAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAGTCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGaa                                                                                                                                                      <  3:342848/145‑1 (MQ=255)      gAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAGTCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGaa                                                                                                                                                      >  4:342848/1‑145 (MQ=255)       aGGCGGTAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGa                                                                                                                                                 >  2:55370/1‑149 (MQ=255)           ggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGCTTCCCGACCgg                                                                                                                                             <  1:420687/149‑1 (MQ=255)           ggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGCTTCCCGACCgg                                                                                                                                             <  1:427261/149‑1 (MQ=255)                 cTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGttt                                                                                                                                       <  2:186322/149‑1 (MQ=255)                  tGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCAGGAGAAAAAAAAAACCCCCGGTCGGTAAGCCGACCGGGGGtt                                                                                                                                        >  1:364209/1‑149 (MQ=255)                  tGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGCTTCCCGACCGGGGGttt                                                                                                                                       >  5:192526/1‑148 (MQ=255)                   ggCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCgggg                                                                                                                                           <  1:193953/143‑1 (MQ=255)                   ggCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCgggg                                                                                                                                           >  2:193953/1‑143 (MQ=255)                    gCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGa                                                                                                                                     <  8:184868/148‑1 (MQ=255)                    gCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGCTTCCCGACCGGGGGttt                                                                                                                                       <  6:192526/146‑1 (MQ=255)                     cccGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAg                                                                                                                                    >  2:309562/1‑148 (MQ=255)                      ccGAACGGGATTGTCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTa                                                                                                                                  <  4:370313/149‑1 (MQ=255)                      ccGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGTAAGCCGACCGGGGGTTTGAGTa                                                                                                                                  >  8:234560/1‑149 (MQ=255)                       cGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTAt                                                                                                                                 >  2:305999/1‑149 (MQ=255)                          aCGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATtct                                                                                                                              <  8:89416/149‑1 (MQ=255)                              gATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATg                                                                                                                          >  2:11900/1‑149 (MQ=255)                                   gCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTcgg                                                                                                                     <  6:311526/149‑1 (MQ=255)                                    cTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTcggc                                                                                                                    >  4:308455/1‑149 (MQ=255)                                    cTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAACGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTcggc                                                                                                                    >  3:9326/1‑149 (MQ=255)                                     tAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTcggcg                                                                                                                   >  6:395012/1‑149 (MQ=255)                                      aaaTCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTcggcgg                                                                                                                  >  5:263610/1‑149 (MQ=255)                                         tCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGccc                                                                                                               <  5:124295/149‑1 (MQ=255)                                         tCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAAACCACGGTCGGGAAGCCGACCGGGGGTTTGACTATTCTCATGTTCGGCGGccc                                                                                                               >  8:85270/1‑149 (MQ=255)                                           gCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCt                                                                                                             <  4:147012/149‑1 (MQ=255)                                             cAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑c                                                                                                           >  7:116251/1‑148 (MQ=255)                                                 cGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAgg                                                                                                      <  5:395012/149‑1 (MQ=255)                                                   tAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTg                                                                                                    <  6:285687/149‑1 (MQ=255)                                                      gtgtTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGCCCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGAc                                                                                                 >  6:72748/1‑149 (MQ=255)                                                        gtTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGc                                                                                               <  3:308455/149‑1 (MQ=255)                                                              gtagtaATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTcgg                                                                                                                     <  3:342356/122‑1 (MQ=255)                                                              gtagtaATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTcgg                                                                                                                     <  3:129274/122‑1 (MQ=255)                                                              gtagtaATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTcgg                                                                                                                     >  4:342356/1‑122 (MQ=255)                                                              gtagtaATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTcgg                                                                                                                     >  4:129274/1‑122 (MQ=255)                                                              gtagtaATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGCCGGGAAGCAGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCgg‑gtg                                                                                        >  8:227790/1‑149 (MQ=255)                                                               tagtaATACGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCgg‑gtgg                                                                                       <  7:379340/149‑1 (MQ=255)                                                                 gtaATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCgg‑gtggg                                                                                      >  7:73221/1‑148 (MQ=255)                                                                    ataaGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTAt                                                                                   <  6:263610/148‑1 (MQ=255)                                                                     taaGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCg                                                                                 >  3:66465/1‑149 (MQ=255)                                                                       aGTGGCGTAGTTCATGAATCTGTCGTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACGGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGgc                                                                               <  1:55370/149‑1 (MQ=255)                                                                           gCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGt                                                                                                     <  1:376014/124‑1 (MQ=255)                                                                           gCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGt                                                                                                     >  2:376014/1‑124 (MQ=255)                                                                             gTATTTCATTAATCTGTCCTCAAACTGAATCGGTGTTGCTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCTCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAg                                                                         <  2:196479/149‑1 (MQ=255)                                                                             gTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAg                                                                         >  8:118825/1‑149 (MQ=255)                                                                                      gAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAAc                                                                >  8:218058/1‑149 (MQ=255)                                                                                           tGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTa                                                                          >  3:223238/1‑134 (MQ=255)                                                                                           tGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTa                                                                          <  4:223238/134‑1 (MQ=255)                                                                                             tCCTCAAAGTTAATCGGTGTTGGTCAAAGGC‑‑CAGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGTCGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAacac                                                          <  6:35493/148‑1 (MQ=255)                                                                                             tCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCaa                                                             <  8:116251/145‑1 (MQ=255)                                                                                             tCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAAt                                                            <  1:305999/146‑1 (MQ=255)                                                                                               cTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATAccc                                                        >  3:154105/1‑148 (MQ=255)                                                                                                tCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCa                                                       >  5:247856/1‑148 (MQ=255)                                                                                                   aaGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAg                                                    >  2:5328/1‑148 (MQ=255)                                                                                                   aaGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGt                                                   >  2:42570/1‑149 (MQ=255)                                                                                                        aaTCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGGTGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATa                                               <  1:309562/149‑1 (MQ=255)                                                                                                         aTCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGc                                             >  7:103982/1‑149 (MQ=255)                                                                                                            ggTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTgg                                          <  2:364209/149‑1 (MQ=255)                                                                                                             gtgtTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGc                                         <  7:234560/149‑1 (MQ=255)                                                                                                                ttGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGACGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAgg                                      >  1:48603/1‑149 (MQ=255)                                                                                                                  ggTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGt                                     <  7:118825/148‑1 (MQ=255)                                                                                                                   gTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAAc                                                                <  8:365435/120‑1 (MQ=255)                                                                                                                   gTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAAc                                                                >  7:365435/1‑120 (MQ=255)                                                                                                                   gTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGTTGc                                   >  2:268847/1‑149 (MQ=255)                                                                                                                    tCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGc                                   >  3:63966/1‑148 (MQ=255)                                                                                                                     cAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCa                                 >  4:84386/1‑149 (MQ=255)                                                                                                                      aaaGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCAc                                >  6:412765/1‑149 (MQ=255)                                                                                                                      aaaGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCAc                                >  1:296037/1‑149 (MQ=255)                                                                                                                        aGGC‑‑CTGAAAAACAAAACCCCCGGCCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCCGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGc                              <  7:227790/149‑1 (MQ=255)                                                                                                                          gC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAg                            >  6:94852/1‑149 (MQ=255)                                                                                                                          gC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAg                            >  6:183290/1‑149 (MQ=255)                                                                                                                                gAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGAGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCgg                        >  6:416240/1‑149 (MQ=255)                                                                                                                                         aacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCa               <  7:218058/149‑1 (MQ=255)                                                                                                                                          acCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAg              >  3:195164/1‑149 (MQ=255)                                                                                                                                             cccGGTCGGCTTCCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTACGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGca            >  1:331948/1‑149 (MQ=255)                                                                                                                                                 gTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGc                                   >  2:252072/1‑121 (MQ=255)                                                                                                                                                 gTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGc                                   <  1:252072/121‑1 (MQ=255)                                                                                                                                                      gAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAg  >  7:107047/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                                      GATCGGAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGC‑‑CTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCT‑CGAAGGTGGACGCCGG‑GTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAG  >  NC_002947/3465908‑3466203

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