Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,266,719	G→C	7.4%	S162T (AGC→ACC) ‡	dedD →	cell division protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,266,719	0	G	C	7.4%	207.0 / 5.0	67	S162T (AGC→ACC) ‡	dedD	cell division protein
Reads supporting (aligned to +/- strand):  ref base G (36/26);  new base C (3/2);  total (39/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

GTGGAGGCCAAAGCTGCTGCAACACCAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTCAAGGGCTTCGTGGTGCGCTTCCAGC  >  NC_002947/2266577‑2266865                                                                                                                                               |                                                                                                                                                   gTGGAGGCCAAAGCTGCTGCAACACCAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAggg                                                                                                                                              <  2:237714/149‑1 (MQ=255)  tGGAGGCCAAAGCTGCTGCAACACCAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGt                                                                                                                                             <  1:84475/149‑1 (MQ=255)    gAGGCCAAAGCTGCTGCAACACCAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTa                                                                                                                                           >  5:30045/1‑149 (MQ=255)     aGGCCAAAGCTGCTGCAACACCAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGGCTACGCAGGGtttg                                                                                                                                          <  8:78951/149‑3 (MQ=255)          aaaGCTGCTGCAACACCAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGGCTACGCAGGGtttgctgca                                                                                                                                     <  1:283543/149‑8 (MQ=255)          aaaGCTGCTGCAACACCAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGc                                                                                                                                     >  1:117571/1‑149 (MQ=255)                    cAACACCAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCg                                                                                                                            <  1:160230/148‑1 (MQ=255)                    cAACACCAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCg                                                                                                                            >  4:14878/1‑148 (MQ=255)                     aaCACCAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTa                                                                                                                          >  7:102304/1‑149 (MQ=255)                      acacCAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTAc                                                                                                                         >  2:372451/1‑149 (MQ=255)                       cacCAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTAcg                                                                                                                        >  8:14126/1‑149 (MQ=255)                       cacCAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATACAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTCCAGCAAACCCTGCGAAGCCAGGGTTACAACGCCTATATCCGTAcg                                                                                                                        >  6:95719/1‑149 (MQ=255)                          cAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTAcggcg                                                                                                                     >  2:87120/1‑149 (MQ=255)                          cAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTAcggcg                                                                                                                     >  6:208990/1‑149 (MQ=255)                          cAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTAcggcg                                                                                                                     >  8:59638/1‑149 (MQ=255)                               ccTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTgg                                                                                                                <  5:208990/149‑1 (MQ=255)                                cTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCggg                                                                                                                   <  7:253400/145‑1 (MQ=255)                                cTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCggg                                                                                                                   >  8:253400/1‑145 (MQ=255)                                   ccGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCAt                                                                                                             >  6:62712/1‑148 (MQ=255)                                           tCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGt                                                                                                    <  4:74731/149‑1 (MQ=255)                                           tCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGt                                                                                                    >  3:253691/1‑149 (MQ=255)                                              aaGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTa                                                                                                 >  1:349732/1‑149 (MQ=255)                                               aGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTAc                                                                                                <  8:102304/149‑1 (MQ=255)                                                 aTCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGt                                                                                              <  7:14126/149‑1 (MQ=255)                                                  tCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAATCCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTg                                                                                             <  5:95719/149‑1 (MQ=255)                                                    gACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTgg                                                                                            <  2:349732/148‑1 (MQ=255)                                                         cAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGt                                                                                      <  3:14878/149‑1 (MQ=255)                                                            tGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGa                                                                                   <  2:117571/149‑1 (MQ=255)                                                               ccTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATc                                                                                 >  2:66027/1‑148 (MQ=255)                                                               ccTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCg                                                                                <  5:62712/149‑1 (MQ=255)                                                                     cGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAgcgcg                                                                          >  3:370051/1‑149 (MQ=255)                                                                          gTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCgccgag                                                                     >  4:250465/1‑149 (MQ=255)                                                                           ttGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCgccgagg                                                                    >  1:33299/1‑149 (MQ=255)                                                                                 gATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAgc                                                              >  3:104711/1‑149 (MQ=255)                                                                                 gATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAgc                                                              >  2:19322/1‑149 (MQ=255)                                                                                     cAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCgcgt                                                          <  8:222729/149‑1 (MQ=255)                                                                                      aGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGa                                                                                   <  2:242159/123‑1 (MQ=255)                                                                                      aGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGa                                                                                   >  1:242159/1‑123 (MQ=255)                                                                                       gCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAgc                                                              <  4:236344/143‑1 (MQ=255)                                                                                       gCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAgc                                                              >  3:236344/1‑143 (MQ=255)                                                                                       gCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCgcgtgc                                                        <  6:30045/149‑1 (MQ=255)                                                                                           gCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGa                                                    >  2:48841/1‑149 (MQ=255)                                                                                             cAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGtt                                                                                     <  2:361164/114‑1 (MQ=255)                                                                                             cAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGtt                                                                                     >  1:361164/1‑114 (MQ=255)                                                                                              aGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGc                                                 <  1:372451/149‑1 (MQ=255)                                                                                                    tCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATc                                             <  2:187202/147‑1 (MQ=255)                                                                                                    tCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGAGACGAGGCCGAGCGCGTGCGTGACGCCATCaa                                           >  8:179206/1‑149 (MQ=255)                                                                                                      cAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAAcc                                         <  4:104711/149‑1 (MQ=255)                                                                                                      cAACCGCGCACGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGACCGAGCGCGTGCGTGACGCCATCAAcc                                         >  5:88590/1‑149 (MQ=255)                                                                                                       aacaccgcAGGAGCCGAAAAGCTGCACCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTATACGAGGGGACGTTGATCGAGCGCCACGACGCCGAGCGCGTGCGTGACGCCAACAACCg                                        >  2:269264/6‑149 (MQ=255)                                                                                                       aaCCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCg                                        >  8:329555/1‑149 (MQ=255)                                                                                                                 gAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGCCGCCATCAACCGTCAGAACAGc                              <  7:59638/149‑1 (MQ=255)                                                                                                                     cGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTc                           >  8:343371/1‑148 (MQ=255)                                                                                                                     cGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTc                           >  4:155426/1‑148 (MQ=255)                                                                                                                      gAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTc                           >  1:173523/1‑147 (MQ=255)                                                                                                                       aaaaGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTCAAg                        >  7:326696/1‑149 (MQ=255)                                                                                                                               gttgtaaCCCTGGCTACGCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTCAAGGGCTTCgt                >  5:152713/6‑149 (MQ=255)                                                                                                                                  gtaaCCCTGGCTACGCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGGCTCAAGGGCTTCgtggt             >  2:356003/3‑149 (MQ=255)                                                                                                                                  caAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTCAAGGGCTTCgtggt             <  1:66027/149‑1 (MQ=255)                                                                                                                                  caAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTCAAGGGCTTCgtggt             >  6:47440/1‑149 (MQ=255)                                                                                                                                   taaCCCTGGCTACGCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTCAAGGGCTTCgtggtg            >  4:64564/2‑149 (MQ=255)                                                                                                                                   aaaCCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTCAAGGGCTTCgtggtg            <  4:253691/149‑1 (MQ=255)                                                                                                                                       ccTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTCAAGGGCTTCGTGGTGCGCt        >  7:313456/1‑149 (MQ=255)                                                                                                                                       ccTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTCAAGGGCTTCGTGGTGCGCt        >  6:86254/1‑149 (MQ=255)                                                                                                                                         tGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTCAAGGGCTTCGTGGTGCGCTTc      >  2:278215/1‑149 (MQ=255)                                                                                                                                             tAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTCAAGGGCTTCGTGGTGCGCTTCCAGc  <  2:33299/149‑1 (MQ=255)                                                                                                                                             tAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGACCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTCAAGGGCTTCGTGGTGCGCTTCCAGc  <  6:88590/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                   GTGGAGGCCAAAGCTGCTGCAACACCAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTCAAGGGCTTCGTGGTGCGCTTCCAGC  >  NC_002947/2266577‑2266865

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