Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,786,910	Δ1 bp	5.6%	coding (741/864 nt)	tsf →	elongation factor Ts

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,786,910	0	G	.	5.6%	231.4 / 6.6	71	coding (741/864 nt)	tsf	elongation factor Ts
Reads supporting (aligned to +/- strand):  ref base G (35/32);  new base . (2/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

CAGCTGAACGCCGACAAGATCGCCGGCAAGCCGGAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCA  >  NC_002947/1786767‑1787055                                                                                                                                                 |                                                                                                                                                  cAGCTGAACGCCGACAAGATCGCCGGCAAGCCGGAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCa                                                                                                                                               <  8:152630/148‑1 (MQ=255)             gACAAGATCGCCGGCAAGCCGGAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTg                                                                                                                                   >  1:298019/1‑148 (MQ=255)                   aTCGCCGGCAAGCCGGAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAAt                                                                                                                            <  1:296444/149‑1 (MQ=255)                          gCAAGCCGGAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCg                                                                                                                          <  8:37667/144‑1 (MQ=255)                          gCAAGCCGGAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCg                                                                                                                          >  7:37667/1‑144 (MQ=255)                            aaGCCGGAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCtt                                                                                                                   <  6:184540/149‑1 (MQ=255)                            aaGCCGGAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCtt                                                                                                                   >  1:295238/1‑149 (MQ=255)                            aaGCCGGAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCt                                                                                                                    >  3:24914/1‑148 (MQ=255)                              gCCGGAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCa                                                                                                                 >  3:62414/1‑149 (MQ=255)                                  gAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCAcc                                                                                                               <  1:123589/147‑1 (MQ=255)                                         tCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                      <  3:54609/149‑1 (MQ=255)                                         tCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                      >  3:295840/1‑149 (MQ=255)                                           gTTTAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGt                                                                                                    <  5:109358/149‑1 (MQ=255)                                            ttGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGaaa                                                                                                                             <  1:289609/123‑1 (MQ=255)                                            ttGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGaaa                                                                                                                             >  2:289609/1‑123 (MQ=255)                                             tGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                      >  1:305653/1‑145 (MQ=255)                                             tGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                      <  2:305653/145‑1 (MQ=255)                                               agaACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAgg                                                                                                 <  8:271106/148‑1 (MQ=255)                                                  aCATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCg                                                                                                                                      >  3:148617/1‑108 (MQ=255)                                                  aCATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCg                                                                                                                                      <  4:148617/108‑1 (MQ=255)                                                    aTGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAAt                                                                                                                            <  7:16100/116‑1 (MQ=255)                                                    aTGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAAt                                                                                                                            >  8:16100/1‑116 (MQ=255)                                                       aTCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCAt                                                                                        >  3:171271/1‑149 (MQ=255)                                                        tCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATc                                                                                       <  4:295840/149‑1 (MQ=255)                                                         cAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGTCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATc                                                                                       >  1:237893/1‑149 (MQ=255)                                                           aCGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCa                                                                                                                                 <  1:103140/104‑1 (MQ=255)                                                           aCGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCa                                                                                                                                 >  2:103140/1‑104 (MQ=255)                                                           aCGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTa                                                                                                   <  5:163596/134‑1 (MQ=255)                                                           aCGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCAACTACTTCAAAGTa                                                                                                   >  6:163596/1‑134 (MQ=255)                                                               tCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAAGCAAGGTTCGCGAGCTGGCCAAGAAAGCCGGGGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGc                                                                                <  2:24965/149‑1 (MQ=255)                                                                              cAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTc                                                                  <  4:62414/148‑1 (MQ=255)                                                                               aaGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGACAGCCGGTGCAGAACTCGTTTCCTGCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGCCGACGTCgc                                                                >  7:139696/1‑149 (MQ=255)                                                                               aaGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCgc                                                                >  1:104673/1‑149 (MQ=255)                                                                               aaGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCgc                                                                >  7:143789/1‑149 (MQ=255)                                                                                   ccGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTcgac                                                                        <  2:73353/137‑1 (MQ=255)                                                                                   ccGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTcgac                                                                        >  1:73353/1‑137 (MQ=255)                                                                                   ccGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCgctgct                                                            <  2:298019/149‑1 (MQ=255)                                                                                     gAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGctgctg                                                           >  3:159475/1‑148 (MQ=255)                                                                                      aaGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGaa                                                         <  2:295238/149‑1 (MQ=255)                                                                                      aaGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGaa                                                         >  7:275503/1‑149 (MQ=255)                                                                                           tCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGc                                                    >  2:109635/1‑149 (MQ=255)                                                                                                 aaaGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGCCGAAGGCATCGAGAAGCCAGTCGAAGACTTCGCTGCTGAAGTTGCCGCt                                                >  7:59930/1‑147 (MQ=255)                                                                                                       cAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCgc                                        <  4:171271/149‑1 (MQ=255)                                                                                                       cAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGAGAAGGCATCGAGAAGCCAGTCGAAGACTTCGCTGCTGAAGTTGCAGCTCAGGTCgc                                        >  5:246304/1‑149 (MQ=255)                                                                                                         actccTTCGTCATTAACCCGGAAATCAAGGTTGGCGATCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCgct                                       <  6:246304/146‑1 (MQ=255)                                                                                                          g‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCgctgc                                     >  2:228472/1‑149 (MQ=255)                                                                                                             cTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGc                                                    >  5:223021/1‑132 (MQ=255)                                                                                                             cTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGc                                                    <  6:223021/132‑1 (MQ=255)                                                                                                             cTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGctgctg                                   >  1:71968/1‑149 (MQ=255)                                                                                                                    aTGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAg                              >  6:9398/1‑147 (MQ=255)                                                                                                                        aCCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGt                          >  6:30220/1‑147 (MQ=255)                                                                                                                         cccGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAg                       <  5:9398/149‑1 (MQ=255)                                                                                                                            ggAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGa                      >  8:169259/1‑147 (MQ=255)                                                                                                                             gAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAg                   >  2:242677/1‑149 (MQ=255)                                                                                                                               aaTCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGcc                 <  8:143789/149‑1 (MQ=255)                                                                                                                               aaTCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGGGCAGAAATCGTTTCCTTCACCTCCTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGcc                 <  8:139696/149‑1 (MQ=255)                                                                                                                                aTCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCgc                                        >  6:306775/1‑125 (MQ=255)                                                                                                                                aTCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCgc                                        <  5:306775/125‑1 (MQ=255)                                                                                                                                 tCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAAGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGcccc               >  6:24804/1‑149 (MQ=255)                                                                                                                                  cAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCg              >  2:158948/1‑149 (MQ=255)                                                                                                                                   aaGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGt             <  1:109635/149‑1 (MQ=255)                                                                                                                                   aaGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGt             <  1:242677/149‑1 (MQ=255)                                                                                                                                     ggTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCt           <  7:8853/149‑1 (MQ=255)                                                                                                                                         ggCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTc        <  1:228472/148‑1 (MQ=255)                                                                                                                                         ggCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTc        <  5:30220/148‑1 (MQ=255)                                                                                                                                         ggCGAGCT‑GCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAg                              <  3:268181/125‑1 (MQ=255)                                                                                                                                         ggCGAGCT‑GCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAg                              <  5:165688/125‑1 (MQ=255)                                                                                                                                         ggCGAGCT‑GCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAg                              >  6:165688/1‑125 (MQ=255)                                                                                                                                         ggCGAGCT‑GCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAg                              >  4:268181/1‑125 (MQ=255)                                                                                                                                           cGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGcc     >  8:25608/1‑149 (MQ=255)                                                                                                                                                tGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCa  >  7:117138/1‑147 (MQ=255)                                                                                                                                                 ggCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGcccc   <  5:274708/145‑1 (MQ=255)                                                                                                                                                 ggCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGcccc   >  6:274708/1‑145 (MQ=255)                                                                                                                                                 ggCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCa  >  5:205165/1‑146 (MQ=255)                                                                                                                                                 |                                                                                                                                                  CAGCTGAACGCCGACAAGATCGCCGGCAAGCCGGAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAG‑CCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCA  >  NC_002947/1786767‑1787055

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