Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,221,235	A→C	9.1%	intergenic (+50/+212)	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,235	0	A	C	9.1%	210.6 / 11.7	77	intergenic (+50/+212)	PP_3699/PP_3700	hypothetical protein/chromosome partitioning ATPase
Reads supporting (aligned to +/- strand):  ref base A (35/35);  new base C (3/4);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

TTGCGTGCAACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACATTTCCGCGGCCAA  >  NC_002947/4221089‑4221383                                                                                                                                                   |                                                                                                                                                      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                                                                                                                                      <  2:934823/131‑1 (MQ=255)                                  tGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctcagc                                                                                                                                      >  1:934823/1‑131 (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                                                                                                                   <  1:1005220/106‑1 (MQ=255)                                                                             tttGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑aaa                                                                                                                   >  2:1005220/1‑106 (MQ=255)                                                                             tttGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGaa                                                                                         <  1:1699928/132‑1 (MQ=255)                                                                             tttGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGaa                                                                                         >  2:1699928/1‑132 (MQ=255)                                                                             tttGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAAc                                                                       <  2:1369420/150‑1 (MQ=255)                                                                             tttGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAAc                                                                       <  2:1314647/150‑1 (MQ=255)                                                                                 aaCAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCTGTCCGATGccccc                                                                                                                                              <  2:1295211/76‑2 (MQ=38)                                                                                 aaCAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCTGTCCGATGccccc                                                                                                                                              >  1:1295211/1‑75 (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                                                                                 <  2:1636651/133‑1 (MQ=255)                                                                                    aaGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAAt                                                                                 >  1:1636651/1‑133 (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)                                                                                                                                                  cAGATACCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACATTTCCGCGGCCaa  <  1:19689/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      TTGCGTGCAACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTC‑AAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACATTTCCGCGGCCAA  >  NC_002947/4221089‑4221383

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