Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,786,986	A→G	8.5%	T259A (ACC→GCC)	tsf →	elongation factor Ts

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,786,986	0	A	G	8.5%	192.6 / 4.9	68	T259A (ACC→GCC)	tsf	elongation factor Ts
Reads supporting (aligned to +/- strand):  ref base A (25/37);  new base G (4/2);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGGAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGC‑TCACGCGCGCAGCCTCTTTGT  >  NC_002947/1786838‑1787132                                                                                                                                                     |                                                                                                                                                    gCCGGAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCAc                                                                                                                                                    <  1:329890/150‑1 (MQ=255) gCCGGAAAACATCGTTGAGAACATGATCAACGGTCCTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCCTCCTCCGGCACACCAGTTTCAAGGCTGGCGGGGTTATCGAGAAGGCCGGTGCAGAAATCGTTTCCTTCAc                                                                                                                                                    <  1:739694/150‑1 (MQ=255)     gAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTAc                                                                                                                                                <  2:1277665/150‑1 (MQ=255)           aTCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                                                          <  1:234418/150‑1 (MQ=255)           aTCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                                                          <  2:288326/150‑1 (MQ=255)           aTCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                                                          <  2:645524/150‑1 (MQ=255)           aTCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                                                          <  2:598108/150‑1 (MQ=255)            tCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAg                                                                                                                                         >  2:125254/1‑150 (MQ=255)              gTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCCAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTa                                                                                                                                       <  2:1364933/150‑1 (MQ=255)                   gaACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGa                                                                                                                                  <  1:629346/150‑1 (MQ=255)                     aCATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCGCCTACTTTGAAGTAGGTGaa                                                                                                                                 >  2:359790/1‑149 (MQ=255)                     aCATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCGCCTACTTTGAAGTAGGTGaa                                                                                                                                 <  1:359790/149‑1 (MQ=255)                       aTGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTGGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGc                                                                                                                              >  2:343165/1‑150 (MQ=255)                         gATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCAt                                                                                                                            >  1:1100055/1‑150 (MQ=255)                          aTCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAGGTTGGGGTGGAGTGCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATc                                                                                                                           <  1:1316313/150‑1 (MQ=255)                          aTCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATc                                                                                                                           <  1:437417/150‑1 (MQ=255)                            cAACGGTCGTATCACCAATTGCAAGGGCGAAGCCTCGCCGAAAGAGCAAGCCTTCGTCATGACCCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCga                                                                                                                         <  1:102/150‑1 (MQ=255)                            cAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCga                                                                                                                         <  1:125254/150‑1 (MQ=255)                              aCGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGGCTGCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                                                          <  1:1095650/131‑1 (MQ=255)                              aCGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                                                          >  2:1095650/1‑131 (MQ=255)                                ggTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAg                                                                                                                     >  2:1316648/1‑150 (MQ=255)                                  tCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGcc                                                                                                                   <  2:846661/150‑1 (MQ=255)                                   cGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCa                                                                                                                  <  1:867533/150‑1 (MQ=255)                                        cACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTcga                                                                                                             >  1:782878/1‑150 (MQ=255)                                         aCCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTcgac                                                                                                            >  1:309479/1‑150 (MQ=255)                                         aCCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTcgac                                                                                                            >  2:1205574/1‑150 (MQ=255)                                         aCCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTcgac                                                                                                            <  1:1062603/150‑1 (MQ=255)                                          ccAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTcgacg                                                                                                           >  2:42100/1‑150 (MQ=255)                                          ccAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTcgacg                                                                                                           >  2:925831/1‑150 (MQ=255)                                           cAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTcgacga                                                                                                          >  1:362924/1‑150 (MQ=255)                                           cAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTcgacga                                                                                                          <  1:64468/150‑1 (MQ=255)                                                        gAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGaa                                                                                             <  2:1383919/150‑1 (MQ=255)                                                         aaGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAg                                                                                            <  1:1205574/150‑1 (MQ=255)                                                           gCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGTCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGtt                                                                                          <  1:1316648/150‑1 (MQ=255)                                                              tCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGcc                                                                                       <  2:1243860/150‑1 (MQ=255)                                                                gCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAg                                                                                                                                >  1:1061175/1‑107 (MQ=255)                                                                gCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAg                                                                                                                                <  2:1061175/107‑1 (MQ=255)                                                                         gCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCgc                                                                            <  2:362924/150‑1 (MQ=255)                                                                          cAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGaa                                                                                                                                 >  2:359715/1‑96 (MQ=255)                                                                          cAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGaa                                                                                                                                 <  1:359715/96‑1 (MQ=255)                                                                           aaGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                                                          <  1:42516/86‑1 (MQ=255)                                                                           aaGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                                                          >  2:42516/1‑86 (MQ=255)                                                                           aaGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAgg                                                                                                                                     >  2:227893/1‑91 (MQ=255)                                                                           aaGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAgg                                                                                                                                     <  1:227893/91‑1 (MQ=255)                                                                           aaGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCgctg                                                                          >  2:1235301/1‑150 (MQ=255)                                                                             gCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCgctgct                                                                        >  2:863802/1‑150 (MQ=255)                                                                               cTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCAc                                                                                                                                                    >  2:1235109/1‑72 (MQ=255)                                                                               cTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCAc                                                                                                                                                    <  1:1235109/72‑1 (MQ=255)                                                                                ttCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGaa                                                                                                                      >  2:1160850/1‑101 (MQ=255)                                                                                ttCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGaa                                                                                                                      <  1:1160850/101‑1 (MQ=255)                                                                                  cGTCATGACCCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCACAAGTCGTGACTGGAGTTTAGGTCAACGTAGGCGAAGGCATCGACACGCCAGTCGACGCCTTCTCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCaa                                                                   <  1:324526/150‑1 (MQ=255)                                                                                      aTGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAg                                                               <  1:1235301/150‑1 (MQ=255)                                                                                        gAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTa                                                             >  2:575519/1‑150 (MQ=255)                                                                                            ccGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCgct                                                                           <  2:1053059/132‑1 (MQ=255)                                                                                            ccGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCgct                                                                           >  1:1053059/1‑132 (MQ=255)                                                                                                      aGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCt                                               <  2:782878/150‑1 (MQ=255)                                                                                                                 cTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCaaa                                    <  1:42100/150‑1 (MQ=255)                                                                                                                    gCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAagag                                 >  1:1139366/1‑150 (MQ=255)                                                                                                                      cAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAagag                                 >  2:266924/1‑148 (MQ=255)                                                                                                                      cAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAagag                                 <  1:266924/148‑1 (MQ=255)                                                                                                                               cGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGC‑Tc                     >  1:1025518/1‑150 (MQ=255)                                                                                                                                          tCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGC‑TCACGCGCGCAGc          <  2:1025518/150‑1 (MQ=255)                                                                                                                                          tCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGC‑TCACGCGCGCAGc          <  1:925831/150‑1 (MQ=255)                                                                                                                                                  ttCGCCTACTTTGAAGTAGGTGAAGGCATCGAGAAGCCAGTCGACGCATTGTCGGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGC‑TCACGCGCGCAGCCTCTTTGt  >  1:968689/1‑150 (MQ=255)                                                                                                                                                  ttCGCCTACTTTGAAGTAGGTGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGC‑Tc                     >  2:519198/1‑131 (MQ=255)                                                                                                                                                  ttCGCCTACTTTGAAGTAGGTGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGC‑TCAcgc                 >  2:902228/1‑135 (MQ=255)                                                                                                                                                  ttCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGGGCCTGCGCGC‑GCACGCGCGCAGCCTCTTTGt  >  2:373385/1‑150 (MQ=255)                                                                                                                                                  ttCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCGCGCGTGA‑GCGGGCAGCCTCTTTgg  <  1:592199/150‑2 (MQ=38)                                                                                                                                                     |                                                                                                                                                    GCCGGAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGC‑TCACGCGCGCAGCCTCTTTGT  >  NC_002947/1786838‑1787132

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