Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,287,685	A→G	8.8%	Q271R (CAG→CGG) ‡	PP_2016 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,287,685	0	A	G	8.8%	185.9 / 5.6	69	Q271R (CAG→CGG) ‡	PP_2016	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/36);  new base G (4/2);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTACCCGGTCGAGGCGTACCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCACACGCTGGTTCAC  >  NC_002947/2287538‑2287831                                                                                                                                                       |                                                                                                                                                   aCTACCCGGTCGAGGCGTACCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAg                                                                                                                                                   <  4:352080/149‑1 (MQ=255)     cccGGTCGAGGCGTACCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCAACGACACCGGTGTGAGCAATATTCCCAAG‑‑T‑ACACGGGCTACCAGtggc                                                                                                                                               >  8:170182/1‑146 (MQ=255)               gCGTACCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCa                                                                                                                                     <  6:24951/149‑1 (MQ=255)               gCGTACCCGAAGAAGATTCGCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCa                                                                                                                                     <  6:265682/149‑1 (MQ=255)                 gTACCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCAACTGGTAGGCCGTgtacttgc                                                                                                                                   <  1:187721/149‑8 (MQ=255)                   aCCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTa                                                                                                                                  >  1:55255/1‑148 (MQ=255)                     ccGAAGTTGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTaaaa                                                                                                                               <  6:83439/149‑1 (MQ=255)                       gaagaaGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAg                                                                                                                             <  3:78292/149‑1 (MQ=255)                       gaagaaGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAg                                                                                                                             >  8:135475/1‑149 (MQ=255)                       gaagaaGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATATTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAg                                                                                                                             >  3:228450/1‑149 (MQ=255)                             gaTTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATTTTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTg                                                                                                                       <  1:93231/149‑1 (MQ=255)                               ttGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGca                                                                                                                     <  4:59878/149‑1 (MQ=255)                                     gCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCg                                                                                                               >  6:52924/1‑149 (MQ=255)                                        tCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTAc                                                                                                            >  2:304084/1‑149 (MQ=255)                                           aCCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCg                                                                                                          >  1:39/1‑148 (MQ=255)                                             cTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGaa                                                                                                        >  8:21398/1‑148 (MQ=255)                                             cTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGaa                                                                                                        <  7:21398/148‑1 (MQ=255)                                                 aCAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGGCCTGCACATGCGTACCGAACAAcc                                                                                                   <  6:70983/149‑1 (MQ=255)                                                   aaGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTaaaa                                                                                                                               <  7:63673/119‑1 (MQ=255)                                                   aaGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTaaaa                                                                                                                               >  8:63673/1‑119 (MQ=255)                                                     ggACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCaa                                                                                               <  7:135475/149‑1 (MQ=255)                                                                    gggCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAg                                                                                <  6:178627/149‑1 (MQ=255)                                                                    gggCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAg                                                                                <  5:179086/149‑1 (MQ=255)                                                                     ggCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAgg                                                                               <  8:110926/149‑1 (MQ=255)                                                                      gCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGa                                                                              >  1:22474/1‑149 (MQ=255)                                                                       cACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGAt                                                                             >  5:28537/1‑149 (MQ=255)                                                                        aCGCTGCTGGGCATCAAGGGCCATTAACTTATCTTCGACACTGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCATTTGGCCGTGCTCCAGTAAAAATGACCTGCACATTTGTACCGAACAACCGCAAGTGATc                                                                                         <  7:56088/136‑1 (MQ=255)                                                                        aCGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACTGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATc                                                                                         >  8:56088/1‑136 (MQ=255)                                                                        aCGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGAt                                                                             >  7:175791/1‑148 (MQ=255)                                                                                 ggCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGgc                                                                     <  6:365267/147‑1 (MQ=255)                                                                                   cATCAAGGGCCAGTACCTGATCTTCGACTCCTGTGTGATCAATATTAGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACCTGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCg                                                                 <  4:228450/149‑1 (MQ=255)                                                                                   cATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCTCAAGCCTCTCACGTAGGGCGAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGgcg                                                                    <  1:312112/149‑1 (MQ=255)                                                                                    aTCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCg                                                                 >  4:368142/1‑148 (MQ=255)                                                                                         gggCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTAc                                                           <  3:226510/149‑1 (MQ=255)                                                                                             cAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGa                                                       >  5:368626/1‑149 (MQ=255)                                                                                               gTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATc                                                     <  6:368626/149‑1 (MQ=255)                                                                                                tACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATc                                                     >  4:82699/1‑148 (MQ=255)                                                                                                   cTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATcgac                                                  >  2:363634/1‑148 (MQ=255)                                                                                                    tGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGaacaa                                                                                                     >  5:276356/1‑96 (MQ=255)                                                                                                    tGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGaacaa                                                                                                     <  6:276356/96‑1 (MQ=255)                                                                                                      aTCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATcgacga                                                <  5:52924/147‑1 (MQ=255)                                                                                                      aTCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATcgacga                                                <  2:55255/147‑1 (MQ=255)                                                                                                            gACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCAcc                                        >  2:79817/1‑149 (MQ=255)                                                                                                              cacCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTCCACATCCGTACCGAACAACCGCAAGTGATCTACCTAAAGGATTACCAGGCGCCGGAGAACCTGATCAACGACACGCACCt                                       >  7:245953/1‑148 (MQ=255)                                                                                                               acCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGa                                     <  1:362525/149‑1 (MQ=255)                                                                                                                ccGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGAc                                    <  3:182952/149‑1 (MQ=255)                                                                                                                 cGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGa                                     >  3:325095/1‑147 (MQ=255)                                                                                                                 cGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGa                                     <  4:325095/147‑1 (MQ=255)                                                                                                                 cGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGAcc                                   >  5:274340/1‑149 (MQ=255)                                                                                                                      tgATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCg                               <  3:297461/148‑1 (MQ=255)                                                                                                                      tgATCAATATACGCAAG‑‑T‑AGACGGCCTACCTTTTGGCCGTGCTCCAGTAATAAGGACCTGCACATGCGTACCGAACAACCGCTAGTGATCTACCTCATGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGa                              <  7:170182/149‑1 (MQ=255)                                                                                                                        aTCATTATTCGGAAG‑‑T‑ACTCGGTCTAGCAGTTGGCCGTGCTCCATTAAAAAGGAGCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGTATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAg                             <  8:245953/148‑1 (MQ=255)                                                                                                                           aaTATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGt                         <  3:368142/149‑1 (MQ=255)                                                                                                                             tatTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTc                       <  6:28537/149‑1 (MQ=255)                                                                                                                                 cGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAg                    <  8:175791/148‑1 (MQ=255)                                                                                                                                 cGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAg                    <  2:22474/148‑1 (MQ=255)                                                                                                                                   cAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGAc                 >  5:84803/1‑149 (MQ=255)                                                                                                                                     aG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACgaga                                              >  1:179588/1‑118 (MQ=255)                                                                                                                                     aG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACgaga                                              <  2:179588/118‑1 (MQ=255)                                                                                                                                      ctggagcacGGCCAACTGGTAGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGAc                 >  2:30662/7‑149 (MQ=255)                                                                                                                                       tggagcacGGCCAACTGGTAGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCAAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGATGTTCGAGGAc                 >  1:109682/6‑148 (MQ=255)                                                                                                                                           gcacGGCCAACTGGTAGGCCGTGCTCCAGTAAAAAGGACCTGCACATTCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCacac            >  8:143455/2‑149 (MQ=255)                                                                                                                                            cacGGCCAACTGGTAGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGCCGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCACACg           >  6:242463/1‑149 (MQ=255)                                                                                                                                             acGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGa                     >  5:240056/1‑138 (MQ=255)                                                                                                                                             acGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGa                     <  6:240056/138‑1 (MQ=255)                                                                                                                                             acGGCCAACTGGTAGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCACACGc          <  1:30662/149‑1 (MQ=255)                                                                                                                                                  cTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCACACGCTGGtt     <  1:304084/149‑1 (MQ=255)                                                                                                                                                    aCCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCACACGCTGGTTCa   >  6:327807/1‑149 (MQ=255)                                                                                                                                                     ccAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCACACGCTGGTTCAc  <  2:39/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                   ACTACCCGGTCGAGGCGTACCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCACACGCTGGTTCAC  >  NC_002947/2287538‑2287831

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