Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,968,750	G→A	7.0%	A23A (GCC→GCT)	gph ←	phosphoglycolate phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,968,750	0	G	A	7.0%	168.8 / 2.3	57	A23A (GCC→GCT)	gph	phosphoglycolate phosphatase
Reads supporting (aligned to +/- strand):  ref base G (26/27);  new base A (2/2);  total (28/29)
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.42e-01

CAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGTTGACGTCGACGTCACTGTTGAGCTTGTAGTGCTTGGTCAGCG  >  NC_002947/1968603‑1968894                                                                                                                                                         |                                                                                                                                                 cAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑Tg                                                                                                                                                  <  1:150960/148‑1 (MQ=255)    gcgcAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGa                                                                                                                                              >  3:132409/1‑149 (MQ=255)                gAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGgcg                                                                                                                                  >  3:226764/1‑149 (MQ=255)                  aGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGcc                                                                                                                                <  3:48228/149‑1 (MQ=255)                    ccGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCgt                                                                                                                              <  4:225611/149‑1 (MQ=255)                          gCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAg                                                                                                                        >  6:195082/1‑149 (MQ=255)                                gggTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTAGCAAA‑TGGCGATGAAg                                                                                                                                         >  4:129104/1‑126 (MQ=255)                                gggTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTAGCAAA‑TGGCGATGAAg                                                                                                                                         <  3:129104/126‑1 (MQ=255)                                  ggggaTGGCGCACCCAGTGGCTCCCATCGCGCGGGCGCCCCCCCCGATCACGCCCCCCCTCAGGTTGTCGTCCCTTTTCGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAg                                                                                                                    <  1:276420/145‑1 (MQ=255)                                  gTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGc                                                                                                                >  6:110603/1‑149 (MQ=255)                                      aTGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTc                                                                                                            >  8:296292/1‑149 (MQ=255)                                       tGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTcc                                                                                                           <  7:137191/149‑1 (MQ=255)                                       tGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTc                                                                                                            >  5:330828/1‑148 (MQ=255)                                           gAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGt                                                                                                       >  4:9747/1‑149 (MQ=255)                                            aaCGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTc                                                                                                      <  1:26936/149‑1 (MQ=255)                                               gCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGaa                                                                                                   <  5:195082/149‑1 (MQ=255)                                               gCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGAACAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCg                                                                                                     >  1:246815/1‑147 (MQ=255)                                                   tGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAagag                                                                                                >  5:151574/1‑148 (MQ=255)                                                   tGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGt                                                                                               <  2:251515/149‑1 (MQ=255)                                                     gCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTAc                                                                                             >  8:205233/1‑149 (MQ=255)                                                       tcgggtGGCGCCCGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TAGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTg                                                                                           >  2:338983/6‑149 (MQ=255)                                                       tcgggtGGCGCCCGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TAGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTg                                                                                           >  2:333678/6‑149 (MQ=255)                                                         ccATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCt                                                                                         >  4:13144/1‑149 (MQ=255)                                                          cATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTc                                                                                        <  8:245852/149‑1 (MQ=255)                                                              gcgcGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCa                                                                                     >  2:181838/1‑148 (MQ=255)                                                                gcgcgcgcCACCCGAGCCCACGCCCCGGAACAGGTTCTCCTC‑‑‑‑ACCGGCCGGCAGCCCGCGTTCCGCGAGCATGCCCTCGCAAC‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATCTCGAAGAGTACTGCTCGCAAAc                                                                                  <  2:101808/149‑1 (MQ=255)                                                                 cgGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACg                                                                                 <  5:110603/149‑1 (MQ=255)                                                                         aCCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCAt                                                                              >  3:146747/1‑144 (MQ=255)                                                                         aCCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCAt                                                                              <  4:146747/144‑1 (MQ=255)                                                                          cccGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTAt                                                                        <  7:144460/149‑1 (MQ=255)                                                                                   aCGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCAt                                                                              >  2:115086/1‑134 (MQ=255)                                                                                   aCGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCAt                                                                              <  1:115086/134‑1 (MQ=255)                                                                                     gCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAg                                                             <  4:132409/149‑1 (MQ=255)                                                                                       cACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAggg                                                           <  1:181838/149‑1 (MQ=255)                                                                                       cACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAggg                                                           >  2:150908/1‑149 (MQ=255)                                                                                        aCGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGt                                                          >  5:31797/1‑149 (MQ=255)                                                                                               aGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGAt                                                    <  7:205233/148‑1 (MQ=255)                                                                                                 gTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCAt                                                 >  3:326147/1‑149 (MQ=255)                                                                                                   tGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATg                                                >  8:49239/1‑148 (MQ=255)                                                                                                     tcgtc‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGATGAGTACTGCTCGCAAACGCATGCTTATTCCTGGCGCAGGGTCTGGATCATGTAg                                             <  7:49239/149‑1 (MQ=255)                                                                                                      cgtc‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGt                                            >  5:104820/1‑149 (MQ=255)                                                                                                       gtc‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGtt                                           >  4:139578/1‑149 (MQ=255)                                                                                                              aaCGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TAGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGTTga                                         <  1:333678/148‑1 (MQ=255)                                                                                                              aaCGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TAGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGTTga                                         <  1:338983/148‑1 (MQ=255)                                                                                                                cggccggcAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGTTgacgt                                      >  4:2804/1‑149 (MQ=255)                                                                                                                 ggccggcAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGTTgacgtc                                     >  5:171909/1‑149 (MQ=255)                                                                                                                   ccggcAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGTTgacgtcga                                   >  5:302806/1‑149 (MQ=255)                                                                                                                    cggcAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGTTgacgtcgac                                  >  7:295155/1‑149 (MQ=255)                                                                                                                    cggcAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGAGGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGTTgacgtcgac                                  <  4:226282/149‑1 (MQ=255)                                                                                                                      gcAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGTTgacgtcgacgt                                <  6:104820/149‑1 (MQ=255)                                                                                                                        aGCCCGCGTTCGGCGAGCATGGCCTGGCAAAGTGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGTTgacgtcgacgt                                >  1:35803/1‑148 (MQ=255)                                                                                                                                    gCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGTTGACGTCGACGTCACTGTTGAGCTTg                  <  6:302806/149‑1 (MQ=255)                                                                                                                                      gAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGTTGACGTCGACGTCACTGTTGAGCTTGTa                <  6:171909/149‑1 (MQ=255)                                                                                                                                      gAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGTTGACGTCGACGTCACTGTTGAGCTTGTa                <  4:326147/149‑1 (MQ=255)                                                                                                                                       aGCATGGCCTGGCAAA‑TGTCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGTTGACGTCGACGTCACTGTTGAGCTTGTAg               <  6:31797/149‑1 (MQ=255)                                                                                                                                         cATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCCCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGTTGACGTCGACGTCACTGTTTAGCTTGTAGTg             <  2:246815/149‑1 (MQ=255)                                                                                                                                                   cAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGTTGACGTCGACGTCACTGTTGAGCTTGTAGTGCTTGGTCAGc   <  4:226764/149‑1 (MQ=255)                                                                                                                                                    aaa‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGTTGACGTCGACGTCACTGTTGAGCTTGTAGTGCTTGGTCAGCg  <  3:13144/149‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                 CAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTC‑‑‑‑AACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAA‑TGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTCGCAAACGCATGCTTATTCCTCGCGCAGGGTCTGGATCATGTAGTTGACGTCGACGTCACTGTTGAGCTTGTAGTGCTTGGTCAGCG  >  NC_002947/1968603‑1968894

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