Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,661,483	T→G	5.2%	V52G (GTC→GGC)	PP_1456 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,661,483	0	T	G	5.2%	235.3 / 4.0	78	V52G (GTC→GGC)	PP_1456	membrane protein
Reads supporting (aligned to +/- strand):  ref base T (38/36);  new base G (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.11e-01

GATTGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGG  >  NC_002947/1661340‑1661628                                                                                                                                                |                                                                                                                                                    gATTGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCa                                                                                                                                                >  5:95529/1‑149 (MQ=255)  aTTGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCAt                                                                                                                                               >  4:422929/1‑149 (MQ=255)   ttGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATg                                                                                                                                              >  8:238615/1‑149 (MQ=255)      ccTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCtg                                                                                                                                           >  2:309186/1‑149 (MQ=255)       cTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCtgg                                                                                                                                          <  6:65250/149‑1 (MQ=255)       cTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCtgg                                                                                                                                          >  4:112461/1‑149 (MQ=255)            tCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCa                                                                                                                                     >  5:106625/1‑149 (MQ=255)            tCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCa                                                                                                                                     >  3:5209/1‑149 (MQ=255)                 agcGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGt                                                                                                                                <  7:252587/149‑1 (MQ=255)                     tAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTg                                                                                                                            >  2:62005/1‑149 (MQ=255)                       gTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTc                                                                                                                          <  3:422929/149‑1 (MQ=255)                         gCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCtg                                                                                                                                           <  1:172340/130‑1 (MQ=255)                          cTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTa                                                                                                                       <  6:106625/149‑1 (MQ=255)                          cTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTa                                                                                                                       <  5:240690/149‑1 (MQ=255)                            gATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATg                                                                                                                     >  4:306043/1‑149 (MQ=255)                            gATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGACCAGGGCGCCATgac                                                                                                                                            >  5:31864/1‑124 (MQ=255)                            gATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGACCAGGGCGCCATgac                                                                                                                                            <  6:31864/126‑3 (MQ=255)                              tCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGt                                                                                                                    >  4:223612/1‑148 (MQ=255)                                   tACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGt                                                                                                                                                     >  6:96217/1‑110 (MQ=255)                                     cGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCg                                                                                                             <  1:346358/148‑1 (MQ=255)                                     cGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCGGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGGCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCg                                                                                                             <  2:285257/148‑1 (MQ=255)                                              cAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCa                                                                                                   >  5:283253/1‑149 (MQ=255)                                                   ccccGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGcc                                                                                              >  8:334602/1‑149 (MQ=255)                                                     ccGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGCCCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGCCACTGGTGCTGTCCTATGTCACCTCGGAGCCTGCCATCGCCCt                                                                                            >  8:351113/1‑149 (MQ=255)                                                       gCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTgg                                                                                          <  6:95529/149‑1 (MQ=255)                                                             ggACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAgcg                                                                                    >  7:76381/1‑149 (MQ=255)                                                              gACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTc                                                                                                                   <  8:346726/117‑1 (MQ=255)                                                              gACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTc                                                                                                                   >  7:346726/1‑117 (MQ=255)                                                               aCATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGcc                                                                                  >  2:221018/1‑149 (MQ=255)                                                                 atcatcTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGTCCTGGCCAGCGCCAc                                                                                <  1:2354/149‑1 (MQ=255)                                                                   catcTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCACTGGCCAGCTCCACTg                                                                              >  8:319632/1‑149 (MQ=255)                                                                    atcTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTGCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGc                                                                             >  4:127619/1‑149 (MQ=255)                                                                     tcTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCg                                                                            >  3:9906/1‑149 (MQ=255)                                                                         gTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGttct                                                                        >  7:85437/1‑149 (MQ=255)                                                                         gTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGttct                                                                        >  5:307453/1‑149 (MQ=255)                                                                         gTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGttc                                                                         <  1:255225/148‑1 (MQ=255)                                                                         gTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGttc                                                                         <  8:76381/148‑1 (MQ=255)                                                                         gTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGtt                                                                          <  3:306043/147‑1 (MQ=255)                                                                              ccTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGttct‑tc                                                                     <  5:103330/146‑1 (MQ=255)                                                                                  gggCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCGTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTc                                                               <  1:243795/147‑1 (MQ=255)                                                                                  gggCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTcc                                                              <  1:221018/148‑1 (MQ=255)                                                                                     cggcCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGc                                                          >  4:248125/1‑149 (MQ=255)                                                                                     cggcCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGc                                                          <  6:307453/149‑1 (MQ=255)                                                                                     cggcCTGGGGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGGCGCCATGCTGGCTGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGc                                                          <  8:308519/149‑1 (MQ=255)                                                                                              gtTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGAct                                                 >  1:212656/1‑149 (MQ=255)                                                                                                   tCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGActggctg                                            <  8:85437/149‑1 (MQ=255)                                                                                                   tCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGActggctg                                            >  7:94363/1‑149 (MQ=255)                                                                                                      tGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGgtg                                         >  5:178440/1‑149 (MQ=255)                                                                                                      tGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGgtg                                         >  7:340813/1‑149 (MQ=255)                                                                                                        cgcgACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGgtg                                         <  8:55818/147‑1 (MQ=255)                                                                                                           gACATGGTGCATACCCGTTTTGGTCATGTCGCCCTGGTCGCCATGCTGGTGGCTCTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAg                                    <  7:351113/149‑1 (MQ=255)                                                                                                            aCATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGc                                   <  2:332783/149‑1 (MQ=255)                                                                                                            aCATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGc                                   <  3:223612/149‑1 (MQ=255)                                                                                                            aCATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGc                                   <  7:334602/149‑1 (MQ=255)                                                                                                              aTGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCa                                     <  5:264726/145‑1 (MQ=255)                                                                                                              aTGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCa                                     >  6:264726/1‑145 (MQ=255)                                                                                                                 gTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCa                               <  2:212656/148‑1 (MQ=255)                                                                                                                   gCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCAccc                            >  3:402890/1‑149 (MQ=255)                                                                                                                   gCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCAccc                            >  8:204262/1‑149 (MQ=255)                                                                                                                     aGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGt                          >  6:419409/1‑149 (MQ=255)                                                                                                                      gACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTc                         <  1:62005/149‑1 (MQ=255)                                                                                                                       aCCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTCTCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTc                         <  1:309186/149‑1 (MQ=255)                                                                                                                        cccGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCACTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCgg                       <  7:319632/149‑1 (MQ=255)                                                                                                                         ccGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGc                      >  1:43433/1‑149 (MQ=255)                                                                                                                            ttttGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGc                   >  5:170830/1‑149 (MQ=255)                                                                                                                                  tCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGTCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGcgcg             >  6:325704/1‑149 (MQ=255)                                                                                                                                  tCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTgcgc              >  6:243344/1‑148 (MQ=255)                                                                                                                                  gcATGGCGACCAGGGCGCCATGCTGGTGGCACTGGTGCTGGCCTATGTCACGTCGGAGCCGGCCATAGCCCTGGCCAGCGCCACTGCGTGCT‑ACC‑TGTCCGAGCTGAGCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTCcgcg             >  5:357219/2‑149 (MQ=255)                                                                                                                                   cATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGa            >  1:417536/1‑149 (MQ=255)                                                                                                                                   cATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGa            >  4:423264/1‑149 (MQ=255)                                                                                                                                    aTGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCATTTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGa            >  8:66412/1‑149 (MQ=255)                                                                                                                                    aTGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGAc           <  3:127619/149‑1 (MQ=255)                                                                                                                                         gcCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGtt                                                                          >  8:273362/1‑83 (MQ=255)                                                                                                                                         gcCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCCCCCTGGCCAGCGCCACTGCGtt                                                                          <  7:273362/83‑1 (MQ=255)                                                                                                                                          cccTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCtg     <  4:355078/149‑1 (MQ=255)                                                                                                                                          cccTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCtg     <  6:170830/149‑1 (MQ=255)                                                                                                                                           ccTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCtgt    <  4:402890/149‑1 (MQ=255)                                                                                                                                               gTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTgg  <  3:112461/147‑1 (MQ=255)                                                                                                                                               gTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTgg  <  5:195478/147‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                    GATTGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCT‑TCA‑TTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGG  >  NC_002947/1661340‑1661628

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