Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	873,428	T→C	8.9%	L20P (CTT→CCT)	PP_5457 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	873,428	0	T	C	8.9%	178.5 / 4.2	64	L20P (CTT→CCT)	PP_5457	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/28);  new base C (3/3);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTACAATGATGCCC  >  NC_002947/873279‑873577                                                                                                                                                       |                                                                                                                                                      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                                                                                                                  >  1:771462/1‑148 (MQ=255)                                        aGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGt                                                                                                                  <  2:771462/148‑1 (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                                                                  <  2:1205257/150‑1 (MQ=255)                                                                                      cacCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTgggg                                                                  <  1:123962/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                                             <  1:178865/150‑1 (MQ=255)                                                                                                           gAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAg                                             >  2:178865/1‑150 (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:92054/1‑150 (MQ=255)                                                                                                                  aaaGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAAt                                      <  1:645732/150‑1 (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)                                                                                                                                                 gcaTGCCTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCg                                                                                                        >  2:399947/1‑54 (MQ=255)                                                                                                                                                 gcaTGCCTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCg                                                                                                        <  1:399947/54‑1 (MQ=255)                                                                                                                                                   aTGCTTGCCAGCGAGGTGATCGATGTACACGTGGTCGGCCGTGACCTGGTCGGTGCCGCCGGCACGCTGCCCCGTCAGGGCCTGCATAGGCTCGTCGTGACTGATGTTCAGGGCGACTCACCACGGCAGCTGATCGATGCGTACAatgat      >  2:94802/1‑150 (MQ=255)                                                                                                                                                   aTGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTACAatgat      <  1:913137/150‑1 (MQ=255)                                                                                                                                                       ttGCCGGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTACAATGATGccc  >  2:717428/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                      ATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGC‑TGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTACAATGATGCCC  >  NC_002947/873279‑873577

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