Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,266,720	C→G	7.5%	S162R (AGC→AGG) ‡	dedD →	cell division protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,266,720	0	C	G	7.5%	205.7 / 6.6	67	S162R (AGC→AGG) ‡	dedD	cell division protein
Reads supporting (aligned to +/- strand):  ref base C (36/26);  new base G (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 = 2.80e-01

GTGGAGGCCAAAGCTGCTGCAACACCAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTCAAGGGCTTCGTGGTGCGCTTCCAGCCTG  >  NC_002947/2266577‑2266868                                                                                                                                                |                                                                                                                                                     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                                                                                                                        >  8:59638/1‑149 (MQ=255)                          cAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTAcggcg                                                                                                                        >  6:208990/1‑149 (MQ=255)                               ccTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTgg                                                                                                                   <  5:208990/149‑1 (MQ=255)                                cTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCggg                                                                                                                      >  8:253400/1‑145 (MQ=255)                                cTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCggg                                                                                                                      <  7:253400/145‑1 (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                              >  4:155426/1‑148 (MQ=255)                                                                                                                     cGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTc                              >  8:343371/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                >  6:47440/1‑149 (MQ=255)                                                                                                                                  caAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTCAAGGGCTTCgtggt                <  1:66027/149‑1 (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)                                                                                                                                                ccAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGTGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTCAAGGGCTTCGTGGTGCGCTTCCAGCCTg  <  7:179206/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     GTGGAGGCCAAAGCTGCTGCAACACCAGCACCTGCCGCGCCATCGAAGATCGACGCCAATGGCCTGCCCGTCAGTTGGTCGATCCAGCTGGCCAGTCTGTCCAACCGCGCAGGAGCCGAAAAGCTGCAGCAAACCCTGCGTAGCCAGGGTTACAACGCCTATATCCGTACGGCGGGTGGCATGAACCGCGTTTACGTGGGGCCGTTGATCGAGCGCGCCGAGGCCGAGCGCGTGCGTGACGCCATCAACCGTCAGAACAGCCTCAAGGGCTTCGTGGTGCGCTTCCAGCCTG  >  NC_002947/2266577‑2266868

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