Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,466,056	A→T	8.0%	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	8.0%	254.4 / 10.0	87	noncoding (104/119 nt)	PP_mr39	P15 ncRNA
Reads supporting (aligned to +/- strand):  ref base A (43/37);  new base T (4/3);  total (47/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

CGGAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCG  >  NC_002947/3465911‑3466200                                                                                                                                                   |                                                                                                                                                 cGGAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGc                                                                                                                                               >  6:333730/1‑149 (MQ=255)  ggAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGcc                                                                                                                                              <  3:145865/149‑1 (MQ=255)   gAGGCGGCAGGCTGGCCCGACCGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTGCATGAATCTGTCCTCAAAGTGACTCGGTGTTGGTACAAGGCC‑TGACAAACAAAACCCCCGGTCGGGAAccc                                                                                                                                              >  3:81207/1‑148 (MQ=255)   gAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGcc                                                                                                                                              >  5:197326/1‑148 (MQ=255)        ggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCgg                                                                                                                                        >  5:268081/1‑149 (MQ=255)         gcAGGCTGGCCCGAACGGGATGGGCTAAAGAGCCAGGCGTAGGTGGTTGCGTAGTATTAATTGTCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTTGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCggg                                                                                                                                       <  4:81207/149‑1 (MQ=255)          cAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGCTTGCGTAGTAATGAGAGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCgggg                                                                                                                                      <  5:357253/149‑1 (MQ=255)           aGGCTGGCCCGAACGGGATTGGTTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGCTTCCCGACCggggg                                                                                                                                     <  8:117525/149‑1 (MQ=255)              cTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGttt                                                                                                                                  <  1:264131/149‑1 (MQ=255)               tGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGttt                                                                                                                                  >  4:116899/1‑148 (MQ=255)                ggCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGGGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGa                                                                                                                                <  3:198974/149‑1 (MQ=255)                   ccGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTa                                                                                                                             >  2:367476/1‑149 (MQ=255)                        cGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATtctc                                                                                                                        <  2:47509/149‑1 (MQ=255)                           gATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCAt                                                                                                                      >  3:85558/1‑148 (MQ=255)                           gATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCAt                                                                                                                      <  4:85558/148‑1 (MQ=255)                            aTTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAA‑CCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGtt                                                                                                                   >  2:264850/1‑149 (MQ=255)                               ggCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATCCTCATGTTc                                                                                                                  <  3:116899/148‑1 (MQ=255)                                gCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTcgg                                                                                                                >  6:152584/1‑149 (MQ=255)                                   aaaTCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCggg                                                                                                                                       >  2:233585/1‑123 (MQ=255)                                   aaaTCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCggg                                                                                                                                       <  1:233585/123‑1 (MQ=255)                                    aaTCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCTGACCGGGGGTTTGAGTATTCTCATGTTCggcggc                                                                                                            >  5:296868/1‑149 (MQ=255)                                    aaTCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCggcggc                                                                                                            >  4:69733/1‑149 (MQ=255)                                                tAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTg                                                                                                >  3:79540/1‑149 (MQ=255)                                                tAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTg                                                                                                >  1:26306/1‑149 (MQ=255)                                                      tttGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGc                                                                                           <  8:1608/148‑1 (MQ=255)                                                        tGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCg                                                                                         <  4:79540/148‑1 (MQ=255)                                                        gtcGTAGTAATATGGGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCg                                                                                         <  4:251972/146‑1 (MQ=255)                                                         gCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGa                                                                                                                                >  2:333634/1‑108 (MQ=255)                                                         gCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGa                                                                                                                                <  1:333634/108‑1 (MQ=255)                                                           gtagtaATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggtg                                                                                     >  7:127760/1‑149 (MQ=255)                                                            tagtaATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggtgg                                                                                    <  6:306024/149‑1 (MQ=255)                                                             agtaATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggtggg                                                                                   >  6:87146/1‑149 (MQ=255)                                                             agtaATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggtggg                                                                                   >  7:315945/1‑149 (MQ=255)                                                               taataaGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGACGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTa                                                                                 >  6:223201/1‑149 (MQ=255)                                                                aataaGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTAt                                                                                >  3:241015/1‑149 (MQ=255)                                                                    aGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCTGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGgc                                                                            <  6:296868/149‑1 (MQ=255)                                                                    aGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGgc                                                                            >  7:372826/1‑149 (MQ=255)                                                                        gCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggtgg                                                                                    >  7:35951/1‑137 (MQ=255)                                                                        gCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggtgg                                                                                    <  8:35951/137‑1 (MQ=255)                                                                             gTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTgg                                                                   >  2:389372/1‑149 (MQ=255)                                                                             gTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTgg                                                                   <  5:87146/149‑1 (MQ=255)                                                                                cATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCt                                                                 <  3:69733/148‑1 (MQ=255)                                                                                 aTGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTaa                                                               <  5:333730/149‑1 (MQ=255)                                                                                 aTGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTa                                                                >  4:403856/1‑148 (MQ=255)                                                                                 aTGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTa                                                                >  4:93116/1‑148 (MQ=255)                                                                                   gAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAAc                                                             >  7:109376/1‑149 (MQ=255)                                                                                    aaTCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTg                                                                    >  6:394659/1‑141 (MQ=255)                                                                                    aaTCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTg                                                                    <  5:394659/141‑1 (MQ=255)                                                                                       cTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggt                                                                                      <  3:79647/120‑1 (MQ=255)                                                                                       cTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggt                                                                                      >  4:79647/1‑120 (MQ=255)                                                                                          tCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATa                                                        <  8:315945/147‑1 (MQ=255)                                                                                            cTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATAccc                                                     >  7:321664/1‑148 (MQ=255)                                                                                            cTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCa                                                    <  3:286338/149‑1 (MQ=255)                                                                                             tCAAAGTGAATCGGTGTTGGTCAAAGGCCGTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTATCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTATTGGCTAAACCAAAAccc                                                     >  1:31324/1‑148 (MQ=255)                                                                                              cAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCaa                                                          >  6:38364/1‑141 (MQ=255)                                                                                              cAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCaa                                                          <  5:38364/141‑1 (MQ=255)                                                                                                  gTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTaa                                              <  1:367476/149‑1 (MQ=255)                                                                                                             ttGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggtgggt                                                                                  <  6:234799/102‑1 (MQ=255)                                                                                                             ttGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggtgggt                                                                                  >  5:234799/1‑102 (MQ=255)                                                                                                             ttGGTCAAAGGCC‑TGAAAAACAAAA‑CCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTTGCTAAACCAATACCCAAAGTAATAGCTGGCAGGt                                  <  1:264850/149‑1 (MQ=255)                                                                                                               ggTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAATGTGGACGGCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGt                                  <  6:268081/148‑1 (MQ=255)                                                                                                                gTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGt                                  >  6:153096/1‑147 (MQ=255)                                                                                                                gTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGt                                  <  5:153096/147‑1 (MQ=255)                                                                                                                gTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGc                                <  2:31324/149‑1 (MQ=255)                                                                                                                gTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGc                                >  5:337541/1‑149 (MQ=255)                                                                                                                gTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGc                                <  2:310046/149‑1 (MQ=255)                                                                                                                   aaaGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCa                              >  6:163409/1‑148 (MQ=255)                                                                                                                   aaaGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCa                              <  5:163409/148‑1 (MQ=255)                                                                                                                      ggCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCa                          <  5:152584/149‑1 (MQ=255)                                                                                                                         c‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGTCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGcc                       >  2:400982/1‑149 (MQ=255)                                                                                                                             aaaaaCAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCg                      <  2:26306/147‑1 (MQ=255)                                                                                                                               aaacaaaacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCt                                         >  6:127762/1‑126 (MQ=255)                                                                                                                               aaacaaaacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCt                                         <  5:127762/126‑1 (MQ=255)                                                                                                                                 acaaaacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACg                >  1:51156/1‑149 (MQ=255)                                                                                                                                   aaaacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGc              <  1:389372/149‑1 (MQ=255)                                                                                                                                     aacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCa            >  5:232132/1‑149 (MQ=255)                                                                                                                                        ccccGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGca         >  7:6144/1‑149 (MQ=255)                                                                                                                                        ccccGGTCGGCTTCCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATa                                                        <  7:350238/102‑1 (MQ=39)                                                                                                                                        ccccGGTCGGCTTCCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATa                                                        >  8:350238/1‑102 (MQ=39)                                                                                                                                        ccccGGTCGGCTTCCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTgg                                       <  3:196092/119‑1 (MQ=255)                                                                                                                                        ccccGGTCGGCTTCCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTgg                                       >  4:196092/1‑119 (MQ=255)                                                                                                                                        ccccGGTCGGCTTCCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGca         >  2:45855/1‑149 (MQ=255)                                                                                                                                         cccGGTCGGCTTCCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGca         >  6:332099/1‑148 (MQ=255)                                                                                                                                          ccGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAACCCAATACGCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGca         >  5:319951/1‑147 (MQ=255)                                                                                                                                                 ggAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGgcgc                                                                          >  1:393781/1‑75 (MQ=255)                                                                                                                                                 ggAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGgcgc                                                                          <  2:393781/75‑1 (MQ=255)                                                                                                                                                 ggAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACACAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGgcg  >  7:260334/1‑147 (MQ=255)                                                                                                                                                   |                                                                                                                                                 CGGAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCC‑TGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCG  >  NC_002947/3465911‑3466200

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