Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,530,913	G→A	20.3%	G37S (GGC→AGC) ‡	PP_2219 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,530,913	0	G	A	20.3%	84.8 / 18.9	49	G37S (GGC→AGC) ‡	PP_2219	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (14/25);  new base A (6/4);  total (20/29)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-02

CTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACC  >  NC_002947/2530777‑2531059                                                                                                                                               |                                                                                                                                                   cTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaagatcaacggcaacgg                                                                                                                                       >  1:1567375/1‑150 (MQ=9)  tGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaagatcaacggcaacggc                                                                                                                                      <  2:1488254/150‑1 (MQ=9)      gCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaagatcaacggcaacggcaacg                                                                                                                                  <  1:1361471/150‑1 (MQ=255)               ggtGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaagatcaacggcaacggcaacggcaacggca                                                                                                                         >  1:390523/3‑150 (MQ=255)                    cGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaagatcaacggcaacggcaacggcaacggcaacggc                                                                                                                    <  2:293880/150‑1 (MQ=255)                     gggCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATTAAGATCAAGATCAAGATCAAGATCAAGCACGCGTTcaagagcaacggcaacggcaacggcaacggcaacggca                                                                                                                   >  1:1739125/1‑150 (MQ=14)                     gggCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaagatcaacggcaacggcaagggcaacggcagcagca                                                                                                                   >  2:750111/1‑150 (MQ=14)                     gggCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaagatcaacggcaacggcaacggcaacggcaacggca                                                                                                                   >  2:1296203/1‑150 (MQ=255)                      ggCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaacggcaacggcaacggcaacggcaacggcaacggcaa                                                                                                                  >  2:1443053/1‑150 (MQ=255)                       gCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaagatcaacggcaacggcaacggcaacggcaacggcaac                                                                                                                 >  1:1915876/1‑150 (MQ=255)                       gCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaagatcaacggcaacggcaacggcaacggcaacggcaac                                                                                                                 >  2:1458268/1‑150 (MQ=255)                             aTGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaagatcaacggcaacggcaacggcaacggcaacggcaacggcaac                                                                                                           >  1:1905005/1‑150 (MQ=255)                             aTGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaagatcaacggcaacggcaacggcaacggcaacggcaacggcaac                                                                                                           >  2:2348534/1‑150 (MQ=255)                             aTGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGAACGGAAGcaagatcaacggcaacggcaacggcaacggcaacggcaacggcaac                                                                                                           >  1:148497/1‑150 (MQ=14)                                cccGTCACCTTCCTCCGCAAGGCCTCCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCaagatcaacggcaacggcaacggcaacggcaacggcaacggcaacggc                                                                                                        <  2:1078267/150‑1 (MQ=17)                                 ccGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaagatcaacggcaacggcaac                                                                                                                                   <  2:2256058/122‑1 (MQ=255)                                 ccGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaagatcaacggcaacggcaac                                                                                                                                   >  1:2256058/1‑122 (MQ=255)                                 ccGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAagatcaacggcaacggcaacggcaacggcaacggcaacggcaacggca                                                                                                       >  1:2058622/1‑150 (MQ=255)                                      aCCTGCCTCTGCACGGCCTGCGTCACGCTTCCTGAACTCGCG‑‑‑‑‑‑ACGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACGGCAACGGCaacggcaacggcaacggcaacggcaacggcaacggc                                                                                                  <  2:2309936/150‑1 (MQ=17)                                       ccTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACGGCAACGGCAacggcaacggcaacggcaacggcaacggcaacggca                                                                                                 >  1:2309936/1‑150 (MQ=17)                                        cTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACGGCAACGGCAacggcaacggcaacggcaacggcaacggcaacggcaa                                                                                                >  2:612776/1‑149 (MQ=17)                                           ccTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACGGCAACGGCAacggcaacggcaacggcaacggcaacggcaacggcaacag                                                                                             >  2:1768202/1‑146 (MQ=17)                                            cTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACGGCAACGGCAacggcaacggcaacggcaacggcaacggcaacggcaacagc                                                                                            >  1:816481/1‑145 (MQ=17)                                            cTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACGGCAACGGCAacggcaacggcaacggcaacggcaacggcaacggcaacagc                                                                                            <  2:1915876/150‑6 (MQ=17)                                            cTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACGGCAACGGCAacggcaacggcaacggcaacggcaacggcaacggcaacagc                                                                                            >  1:1906218/1‑145 (MQ=17)                                            cTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACGGCAACGGCAacggcaacggcaacggcaacggcaacggcaacggcaacagc                                                                                            <  1:1296203/150‑6 (MQ=17)                                                  aaGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaagatcaacggcaacggcaacggcaacggcaacg                                                                                                                      <  1:528959/118‑1 (MQ=14)                                                  aaGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaagatcaacggcaacggcaacggcaacggcaacg                                                                                                                      >  2:528959/1‑118 (MQ=255)                                                     gCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACGGCAACGTCAACGGCAACGGCAACGGCAACCGCAACGGCAACGGCAACAGCAGGTGGATg                                                                                         <  1:1261629/150‑1 (MQ=17)                                                      ccTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATgg                                                                                        <  1:1458268/150‑1 (MQ=18)                                                        tGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGc                                                                                      <  1:1768202/150‑1 (MQ=18)                                                             tCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCa                                                                           <  1:1443053/150‑1 (MQ=255)                                                                               gaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTgg                                                               >  2:1555320/8‑150 (MQ=255)                                                                                   atcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGcaca                                                           <  2:2058622/147‑1 (MQ=255)                                                                                     caagatcaagatcaagatcaagatcaagatcaagatcaagatcacgatcaagatcAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCAcagc                                                         <  2:1739125/149‑1 (MQ=255)                                                                                     caagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcAACGGCAACGGCAACGGTAACGGCAACAGCAACGGCAACAGCAACGGCAACAGGAGGGAGATGGGCACAGAATCTCAGGTTACTGCTGGcagccg                                                         >  2:983019/2‑146 (MQ=18)                                                                                      aagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaacggcaacggcaacggcaacggc                                                                                                                                >  2:1858577/1‑78 (MQ=39)                                                                                      aagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagagcaacggcaacggcaacggcaacggc                                                                                                                                <  1:1858577/78‑1 (MQ=21)                                                                                       agatcaagatcaagatcaagatcaagatcaagatcaagatcaagacaaagagcaACGGCAACGGCAACGGGTACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCAcagcag                                                       >  1:1624312/1‑150 (MQ=25)                                                                                               atcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCagctgagctg                                               >  1:1108413/1‑150 (MQ=255)                                                                                                 caagatcaagatcaagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTggg                                             <  2:1624312/150‑1 (MQ=255)                                                                                                   agatcaagatcaagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACCCCagctcagctgctgt                                           <  2:663162/150‑5 (MQ=255)                                                                                                   agatcaagatcaagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGACAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACCCCagctcagctgctgt                                           <  1:1959693/150‑5 (MQ=21)                                                                                                     atcaagatcaagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGc                                         <  2:828346/150‑1 (MQ=255)                                                                                                      tcaagatcaagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGcc                                        >  1:904637/1‑150 (MQ=255)                                                                                                      tcaagatcaagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGcc                                        >  2:1381236/1‑150 (MQ=255)                                                                                                      tcaagatcaagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGcc                                        >  1:1639726/1‑150 (MQ=255)                                                                                                       caagatcaagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCg                                       >  1:1791349/1‑150 (MQ=255)                                                                                                        aagatcaagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGa                                      <  2:1567375/150‑1 (MQ=255)                                                                                                          gatcaagatcaagatcaagatcaagatcaagatcaaCAGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACCCCAGCTCAGCTGCTGTGccagta                                    <  1:983019/150‑5 (MQ=21)                                                                                                           atcaagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACcg                                   >  2:1875228/1‑150 (MQ=255)                                                                                                            tcaagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACcgg                                  >  2:2163620/1‑150 (MQ=255)                                                                                                             caagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACcggc                                 <  1:1651138/150‑1 (MQ=255)                                                                                                                gatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGa                              <  2:1143393/150‑1 (MQ=255)                                                                                                                  tcaagctcaaaatcaagctcaataccatCGTCTACACCCACGCGAAAGTCAACGGCAACGGGAACTGTAACGGCAACAGCAGGTGGATGTGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTTAGCTGGGGTGCCGACCGGCGGAtt                            <  2:148497/150‑1 (MQ=255)                                                                                                                    aagatcaagatcaagatcaagatcaaGATCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATtgt                          >  1:1864398/1‑150 (MQ=255)                                                                                                                     agatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATtgtg                         >  1:905806/1‑150 (MQ=255)                                                                                                                       atcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCATTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAg                       <  2:1791349/150‑1 (MQ=255)                                                                                                                       atcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTgg                                                               >  2:2280878/1‑110 (MQ=255)                                                                                                                       atcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTgg                                                               <  1:2280878/110‑1 (MQ=255)                                                                                                                       atcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAg                       <  2:1639726/150‑1 (MQ=255)                                                                                                                       ataaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCACCAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGCGCTGGGGTGCCGACCGGCGGATTGTTAg                       <  2:761035/147‑1 (MQ=255)                                                                                                                              tcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTtactac                                                                  <  2:1695667/100‑1 (MQ=255)                                                                                                                              tcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTtactac                                                                  >  1:1695667/1‑100 (MQ=255)                                                                                                                                aagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAg              <  1:750111/150‑1 (MQ=255)                                                                                                                                aagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAg              <  2:904637/150‑1 (MQ=255)                                                                                                                                       agatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCGACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTc       <  2:1644322/150‑1 (MQ=255)                                                                                                                                          tcaacggcaacggcaacggcaacggcaacggcaacggcaacg                                                                                                                >  1:1798237/1‑42 (MQ=32)                                                                                                                                          tcaacggcaacggcaacggcaacggcaacggcaacggcaacg                                                                                                                <  2:1798237/42‑1 (MQ=32)                                                                                                                                          tcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCa    >  1:1568678/1‑150 (MQ=255)                                                                                                                                          tcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCa    <  1:2348534/150‑1 (MQ=255)                                                                                                                                          tcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCa    <  1:1555320/150‑1 (MQ=255)                                                                                                                                          tcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCa    <  2:1144406/150‑1 (MQ=255)                                                                                                                                            aacggcaacggcaacggcaacggcaacggcaacggcaacggcaacggcAACAGCAGGTGGATGGACACAGAATCTCAGATTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGTCGGATTGGTAGCAGTGGAAGCAAATTCGCAcc  >  2:831317/1‑150 (MQ=38)                                                                                                                                               |                                                                                                                                                   CTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACC  >  NC_002947/2530777‑2531059

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