Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,221,233	A→T	9.2%	intergenic (+48/+214)	PP_3699 → / ← PP_3700	hypothetical protein/chromosome partitioning ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,221,233	0	A	T	9.2%	224.7 / 11.8	76	intergenic (+48/+214)	PP_3699/PP_3700	hypothetical protein/chromosome partitioning ATPase
Reads supporting (aligned to +/- strand):  ref base A (35/34);  new base T (3/4);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-01

GACACTTGCGTGCAACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACATTTCCGCGGCC  >  NC_002947/4221084‑4221381                                                                                                                                                      |                                                                                                                                                      gACACTTGCGTGCAACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGa                                                                                                                                                       <  1:1236824/150‑1 (MQ=255)      ttGCGTGCAACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTAAAGTGGAAAATGCTGACCGTCGTGCATAGCGCGCTCAGGCCGCATCTCACAGAt                                                                                                                                                  <  2:1900034/150‑1 (MQ=255)         cGTGCAACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGcc                                                                                                                                               >  2:2284617/1‑150 (MQ=255)           tGCAACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCTGTCCGATGcccc                                                                                                                                             <  1:1908030/150‑1 (MQ=255)           tGCAACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGcccc                                                                                                                                             >  1:90216/1‑150 (MQ=255)              aaCGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGCAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGCGGGCATCGGACAGATACCCCAAc                                                                                                                                          >  2:19689/1‑150 (MQ=255)                      gTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGcc                                                                                                                                               <  1:288012/137‑1 (MQ=255)                      gTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGcc                                                                                                                                               >  2:288012/1‑137 (MQ=255)                          aCGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctcagcctcagc                                                                                                                              >  2:1533159/1‑150 (MQ=255)                            gCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctcagcctcagcct                                                                                                                            >  1:711811/1‑150 (MQ=255)                                   aCGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCagcctcagcctcagcctc                                                                                                                     <  1:2290516/150‑1 (MQ=255)                                      cTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCTGTCCGATGccccct                                                                                                                                           <  2:2168613/125‑3 (MQ=255)                                      cTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCTGTCCGATGccccct                                                                                                                                           <  2:2166176/125‑3 (MQ=255)                                      cTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCTGTCCGATGccccct                                                                                                                                           >  1:2168613/1‑123 (MQ=255)                                      cTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCTGTCCGATGccccct                                                                                                                                           >  1:2166176/1‑123 (MQ=255)                                       tGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctcagc                                                                                                                                    >  1:934823/1‑131 (MQ=255)                                       tGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctcagc                                                                                                                                    <  2:934823/131‑1 (MQ=255)                                           gATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTa                                                                                                            <  1:2335850/150‑1 (MQ=255)                                            aTCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAGCCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAg                                                                                                           <  2:1403159/150‑1 (MQ=255)                                               cTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTa                                                                                                        >  1:186439/1‑150 (MQ=255)                                                tGAAGCTTTGGTCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTAt                                                                                                       <  2:1625145/150‑1 (MQ=255)                                                tGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctcagc                                                                                                                                    >  2:1179059/1‑122 (MQ=255)                                                tGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctcagc                                                                                                                                    <  1:1179059/122‑1 (MQ=255)                                                  aaGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATgg                                                                                                     <  1:1533159/150‑1 (MQ=255)                                                  aaGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATgg                                                                                                     >  1:1369420/1‑150 (MQ=255)                                                     cTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCg                                                                                                  >  1:1582982/1‑150 (MQ=255)                                                       ttGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAg                                                                                                >  1:2322309/1‑150 (MQ=255)                                                            cAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAgcgcgc                                                                                           <  1:768229/150‑1 (MQ=255)                                                                        gCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAAt                                                                               >  1:731905/1‑150 (MQ=255)                                                                         cTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAAtc                                                                              <  1:2284617/150‑1 (MQ=255)                                                                            gaggagTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCATCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCa                                                                           <  1:897613/150‑1 (MQ=255)                                                                            gaggagTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCa                                                                           <  2:731905/150‑1 (MQ=255)                                                                                 gTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCaa                                                                      >  2:1415009/1‑150 (MQ=255)                                                                                  tttGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑aaa                                                                                                                 >  2:1005220/1‑106 (MQ=255)                                                                                  tttGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑aaa                                                                                                                 <  1:1005220/106‑1 (MQ=255)                                                                                  tttGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGaa                                                                                       >  2:1699928/1‑132 (MQ=255)                                                                                  tttGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGaa                                                                                       <  1:1699928/132‑1 (MQ=255)                                                                                  tttGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAAc                                                                     <  2:1369420/150‑1 (MQ=255)                                                                                  tttGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAAc                                                                     <  2:1314647/150‑1 (MQ=255)                                                                                      aaCAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCTGTCCGATGccccc                                                                                                                                            >  1:1295211/1‑75 (MQ=38)                                                                                      aaCAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCTGTCCGATGccccc                                                                                                                                            <  2:1295211/76‑2 (MQ=38)                                                                                        cAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctcagcctca                                                                                                                                <  2:389667/86‑1 (MQ=255)                                                                                        cAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctcagcctca                                                                                                                                >  1:389667/1‑86 (MQ=255)                                                                                        cAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGc                                                               <  2:1582982/150‑1 (MQ=255)                                                                                         aaGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACaa                                                                                >  1:1952956/1‑132 (MQ=255)                                                                                         aaGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACaa                                                                                <  2:1952956/132‑1 (MQ=255)                                                                                         aaGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAAt                                                                               >  1:1636651/1‑133 (MQ=255)                                                                                         aaGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAAt                                                                               <  2:1636651/133‑1 (MQ=255)                                                                                         aaGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGcc                                                              >  1:2376739/1‑150 (MQ=255)                                                                                         aaGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCCTCGCGTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGcc                                                              >  1:1748411/1‑150 (MQ=255)                                                                                              tGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAACGTAAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGGCACAATCTCATGCAACCGCTGCCGTACg                                                         >  1:41851/1‑150 (MQ=255)                                                                                               gCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGCGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGa                                                        <  2:90216/150‑1 (MQ=255)                                                                                                cAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAg                                                       >  1:441462/1‑150 (MQ=255)                                                                                                         tgGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAg                                              >  2:1823864/1‑150 (MQ=255)                                                                                                           gCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGACAGACGACGCAGTAATGGCAGTCGACCGCTGCCGTACGAGTCGGCACAGTa                                            >  1:2340040/1‑150 (MQ=255)                                                                                                            cAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGTCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTc                                                     <  2:503781/140‑1 (MQ=255)                                                                                                            cAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGTCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTc                                                     >  1:503781/1‑140 (MQ=255)                                                                                                             aagtGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCg                                                    >  1:2405856/4‑140 (MQ=255)                                                                                                             aaaTGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctcagcctcag                                                                                                                               <  2:1587212/66‑1 (MQ=255)                                                                                                             aaaTGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctcagcctcag                                                                                                                               >  1:1587212/1‑66 (MQ=255)                                                                                                             aaaTGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCg                                                    <  2:2405856/140‑1 (MQ=255)                                                                                                                  cTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCAtc                                     <  1:523399/150‑1 (MQ=255)                                                                                                                  cTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGAACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCAtc                                     <  2:1350877/150‑1 (MQ=255)                                                                                                                    gACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCAtctc                                   <  2:186439/150‑1 (MQ=255)                                                                                                                       cgtcgtGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAg                                <  2:441462/150‑1 (MQ=255)                                                                                                                                tAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCa                       >  1:2262646/1‑150 (MQ=255)                                                                                                                                  gcgcgcTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATg                     >  1:316787/1‑150 (MQ=255)                                                                                                                                      gcTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCt                 <  1:1415009/150‑1 (MQ=255)                                                                                                                                          aGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATg                     <  2:2010515/142‑1 (MQ=255)                                                                                                                                          aGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATg                     >  1:2010515/1‑142 (MQ=255)                                                                                                                                               gCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACATTTcc        <  1:1823864/150‑1 (MQ=255)                                                                                                                                                   cGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACATTTCCGCgg    >  2:2179130/1‑150 (MQ=255)                                                                                                                                                    ggACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAAACCTAGCTATGGCGGAGCGCGCCGAAGCCACAATCTCATGCAAc                                                                     >  1:1079985/1‑85 (MQ=255)                                                                                                                                                    ggACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCCCATGCAAc                                                                     <  2:1079985/85‑1 (MQ=255)                                                                                                                                                     gACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCAt                                      <  2:598019/114‑1 (MQ=255)                                                                                                                                                     gACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCAt                                      >  1:598019/1‑114 (MQ=255)                                                                                                                                                     gACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACATTTCCGCGGcc  >  2:734793/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GACACTTGCGTGCAACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACATTTCCGCGGCC  >  NC_002947/4221084‑4221381

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