Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	873,421	G→A	8.2%	G18S (GGC→AGC)	PP_5457 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	873,421	0	G	A	8.2%	173.8 / 5.3	62	G18S (GGC→AGC)	PP_5457	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/28);  new base A (2/3);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CCGAGGATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTACAA  >  NC_002947/873273‑873569                                                                                                                                                      |                                                                                                                                                     ccGAGGATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑Tggcagg                                                                                                                                                     >  1:680349/1‑150 (MQ=255) ccGAGGATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑Tggcagg                                                                                                                                                     >  2:932665/1‑150 (MQ=255)  cGAGGATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑Tggcaagc                                                                                                                                                    >  2:566301/1‑150 (MQ=255)    aGGATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACTAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAAGCAt                                                                                                                                                  <  1:1419799/150‑1 (MQ=255)       aTTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCt                                                                                                                                               >  2:1030147/1‑150 (MQ=255)        ttGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCtt                                                                                                                                              <  1:565622/150‑1 (MQ=255)        ttGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCtt                                                                                                                                              <  2:680349/150‑1 (MQ=255)          gCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAAGCATGCCTGc                                                                                                                                            <  1:566301/150‑1 (MQ=255)           ccAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGcc                                                                                                                                           >  1:650432/1‑150 (MQ=255)             aGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAg                                                                                                                                         >  2:1294561/1‑150 (MQ=255)                tgtgCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGCAAAGCGATGTAGCATATCTACGAACCGGGAAGCCTGCTTCCGGC‑CAAATTCTTTTCAGGGATGCTTGCCAGCGa                                                                                                                                      >  1:1095991/1‑150 (MQ=255)                tgtgCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGa                                                                                                                                      <  1:736771/150‑1 (MQ=255)                  tgCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAgg                                                                                                                                    <  1:1079119/150‑1 (MQ=255)                     tGTCAGCCTGTGGTATGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGa                                                                                                                                 <  1:648789/150‑1 (MQ=255)                          gCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGAt                                                                                                                            <  1:1188244/150‑1 (MQ=255)                             tgtgGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTa                                                                                                                         <  2:1120323/150‑1 (MQ=255)                                   gtgCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAAGCATGCCTGCCAGcat                                                                                                                                      >  2:1145600/1‑129 (MQ=255)                                   gtgCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAAGCATGCCTGCCAGcat                                                                                                                                      <  1:1145600/131‑3 (MQ=255)                                    tgCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGGGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTgg                                                                                                                  >  1:1205257/1‑150 (MQ=255)                                     gCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGt                                                                                                                 >  1:1114488/1‑150 (MQ=255)                                            gCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGt                                                                                                          >  2:1212195/1‑150 (MQ=255)                                            gCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGt                                                                                                          <  1:1212195/150‑1 (MQ=255)                                              aGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGt                                                                                                          <  2:771462/148‑1 (MQ=255)                                              aGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGt                                                                                                          >  1:771462/1‑148 (MQ=255)                                               ggTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGAc                                                                                                       <  1:932665/150‑1 (MQ=255)                                                gTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGAcc                                                                                                      >  1:492600/1‑150 (MQ=255)                                                   cATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAgg                                                                                                                                    <  1:937696/117‑1 (MQ=255)                                                   cATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAgg                                                                                                                                    >  2:937696/1‑117 (MQ=255)                                                        tGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGt                                                                                              <  2:492600/150‑1 (MQ=255)                                                        gtaCGAAGCGATGCTTGCAGTCAGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGt                                                                                              <  1:878249/148‑1 (MQ=255)                                                            gAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTgccg                                                                                          >  2:645732/1‑150 (MQ=255)                                                             aaTGGATGCTTGAAGTACCCCAGGTGGACGACACCTGATGCAGCGTATCTACGAACAGGAACGCCTTCTCGAGGCGTCTTGAC‑CGGGACTTCTTCTTGCCAGCGAGGGGATCGATTTACACGTGGTCGGCCGTGACCTTGTCGGTgccgc                                                                                         <  1:1142296/150‑1 (MQ=255)                                                                      ttGAAGTCCGCCGCGAGGACTACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCt                                                                                >  1:489551/1‑150 (MQ=255)                                                                      ttGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCt                                                                                >  2:913137/1‑150 (MQ=255)                                                                      ttGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCt                                                                                >  2:1213842/1‑150 (MQ=255)                                                                               gccgcGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCa                                                                       <  1:1294561/150‑1 (MQ=255)                                                                                     aGGACGGCACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCct                                                                 >  1:857254/1‑150 (MQ=255)                                                                                       gacgacACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCctgc                                                               >  2:123962/1‑150 (MQ=255)                                                                                       gacgacACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCctgc                                                               <  2:1114488/150‑1 (MQ=255)                                                                                          gacACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGc                                                                   <  1:69677/143‑1 (MQ=255)                                                                                          gacACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGc                                                                   >  2:69677/1‑143 (MQ=255)                                                                                            cacCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTgggg                                                          <  1:123962/150‑1 (MQ=255)                                                                                            cacCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTgggg                                                          <  2:1205257/150‑1 (MQ=255)                                                                                              ccTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCt                                                        >  2:446860/1‑150 (MQ=255)                                                                                                  aTGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGAt                                                                                                                            >  2:1205983/1‑78 (MQ=255)                                                                                                  aTGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGAt                                                                                                                            <  1:1205983/78‑1 (MQ=255)                                                                                                    gCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGAGGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTcgtcgt                                                  <  2:489551/150‑1 (MQ=255)                                                                                                             ccgatctccGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCAAGCGTTGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCtgatga                                         <  2:1046897/143‑1 (MQ=255)                                                                                                             cTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCtgatga                                         >  1:250611/1‑150 (MQ=255)                                                                                                                cGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCa                                      <  1:446860/150‑1 (MQ=255)                                                                                                                 gAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAg                                     >  2:178865/1‑150 (MQ=255)                                                                                                                 gAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAg                                     <  1:178865/150‑1 (MQ=255)                                                                                                                   aCCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGAAGGATATACACGTGGCTTCCTGTGAGCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCGCGTCGTGGCTGATGAGCAGGc                                   >  1:530897/1‑150 (MQ=255)                                                                                                                     cGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAgagatcgg                                 >  1:1046897/1‑143 (MQ=255)                                                                                                                      ggAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTc                                                                                                 >  1:751110/1‑85 (MQ=255)                                                                                                                      ggAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTc                                                                                                 <  2:751110/85‑1 (MQ=255)                                                                                                                       gAAAGCCGGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGaa                               >  2:240216/1‑150 (MQ=255)                                                                                                                        aaaGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAAt                              <  1:645732/150‑1 (MQ=255)                                                                                                                        aaaGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAAt                              >  1:92054/1‑150 (MQ=255)                                                                                                                            ccTGCTTGAGGCGCAAATGC‑TGGCAGGCATGATAGCAAACGACCGGATAGATGTACATATGGGCGGCCGTGTGGTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATgcgc                          >  1:387734/1‑150 (MQ=255)                                                                                                                                cTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACgg                      <  2:250611/150‑1 (MQ=255)                                                                                                                                      gcgcAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTcgtcgt                                                  >  2:492909/1‑116 (MQ=255)                                                                                                                                      gcgcAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTcgtcgt                                                  <  1:492909/116‑1 (MQ=255)                                                                                                                                        gcAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGAt              <  1:1030147/150‑1 (MQ=255)                                                                                                                                           ctcgC‑TGGCAAGCATGCCTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGa           >  2:1209505/4‑150 (MQ=255)                                                                                                                                                   gcaggcaTGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTAc    >  2:1297398/1‑150 (MQ=255)                                                                                                                                                     aggcaTGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTACaa  >  2:917987/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CCGAGGATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTACAA  >  NC_002947/873273‑873569

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