Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,862,212	G→C	9.1%	I204M (ATC→ATG)	gllB ←	4‑oxalmesaconate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,862,212	0	G	C	9.1%	195.4 / 8.3	66	I204M (ATC→ATG)	gllB	4‑oxalmesaconate hydratase
Reads supporting (aligned to +/- strand):  ref base G (28/32);  new base C (3/3);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

ACGCGTGGGAACACGCTCTGCAGGCCTTCGGCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCTTCTCGCCCGGCTTGTAGCCTTCAGCCTGCG  >  NC_002947/2862066‑2862360                                                                                                                                                    |                                                                                                                                                      aCGCGTGGGAACACGCTCTGCAGGCCTTCGGCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATCCAGTGCAt                                                                                                                                                     >  2:312676/1‑149 (MQ=255)      tGGGAACACGCTCTGCAGGCCTTCGGCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATgg‑cgg                                                                                                                                               <  7:12330/149‑1 (MQ=255)      tGGGAACACGCTCTGCAGGCCTTCGGCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATCCAGTGCATgg‑ccg                                                                                                                                               <  1:312676/149‑3 (MQ=255)          aaCACGCTCTGCAGGCCTTCGGCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATa                                                                                                                                           >  2:143682/1‑149 (MQ=255)           acacGCTCTGCAGGCCTTCGGCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATAc                                                                                                                                          >  5:335389/1‑149 (MQ=255)             acGCGCTGCCTGCCTTCGGCGTAGACGATGTTGCGGGCGCTGGGTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACtt                                                                                                                                        <  6:11594/149‑1 (MQ=255)               gCTCTGCAGGCCTTCGGCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACtttt                                                                                                                                      >  6:363281/1‑149 (MQ=255)               gCTCTGCAGGCCTTCGGCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACtttt                                                                                                                                      <  1:309185/149‑1 (MQ=255)                     cAGGCCTTCGGCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGt                                                                                                                                   <  6:199814/146‑1 (MQ=255)                     cAGGCCTTCGGCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGt                                                                                                                                   <  8:286405/146‑1 (MQ=255)                     cAGGCCTTCGGCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTcc                                                                                                                                 >  6:390512/1‑148 (MQ=255)                            tCGGCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACAcc                                                                                                                           >  3:120699/1‑147 (MQ=255)                             cGGCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGa                                                                                                                         <  4:96971/148‑1 (MQ=255)                              ggCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCTGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAg                                                                                                                        <  1:269088/148‑1 (MQ=255)                                   aGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATg                                                                                                                  >  6:149122/1‑149 (MQ=255)                                    gACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGt                                                                                                                 >  5:356319/1‑149 (MQ=255)                                    gACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGt                                                                                                                 >  4:9450/1‑149 (MQ=255)                                       gATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCa                                                                                                              <  5:390512/149‑1 (MQ=255)                                        atgatgCGTGCGCAGGACATGCCCACGTTACGCTTGGCCTGCACGCCGCCCAGCAGGGCTACACGCGTGT‑AGTACTCCCAGAGGTGCTCCTGGCCGGCCATGCACTGGATCG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAg                                                                                                             <  4:81631/149‑1 (MQ=255)                                             gcgtgcgCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAAc                                                                                                         <  2:240119/148‑1 (MQ=255)                                             gcgtgcgCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAAc                                                                                                         <  1:143682/148‑1 (MQ=255)                                                 gcgcTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATCCAGTGCATgg‑ccggc                                                                                                                                             <  2:244515/108‑5 (MQ=255)                                                 gcgcTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATCCAGTGCATgg‑ccggc                                                                                                                                             >  1:244515/1‑104 (MQ=255)                                                     tGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCggg                                                                                                >  2:326216/1‑149 (MQ=255)                                                      ggTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGgc                                                                                               >  4:196924/1‑149 (MQ=255)                                                       gTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGgc                                                                                               <  7:29151/148‑1 (MQ=255)                                                            gCCCACGTTACGCTTGGCCTGCACGCCGAGATGCAGGGCTACAAGGGTGTGAGTACTCCCAGAGATGTTCCTGGCCGGCAATGCACTGGATGG‑CGGCATACTTTTAGTCCCACACAGAGGTAATGTCCAGGAACGTGTCAGGGCGCCAt                                                                                          >  1:291275/1‑149 (MQ=255)                                                             cccACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTc                                                                                        <  6:335389/149‑1 (MQ=255)                                                             cccACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGACGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAg                                                                                                             >  6:370672/1‑128 (MQ=255)                                                             cccACGTTACGCTTGCCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAg                                                                                                             <  5:370672/128‑1 (MQ=255)                                                                       gcatGGCCTGCACGCCGCCCTGCAGGGCTACACGGGTGT‑AGTACACCCAGAGGTGTTACATGCCGGCCATGCACTGGAAGG‑CGGCATACTTTTTGTACCACACCGAGGTAATGTACAGGAACCTGTCCCGGCGCCATACGCACTGCTCg                                                                              >  4:61178/4‑149 (MQ=255)                                                                            gCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCt                                                                          <  4:89784/148‑1 (MQ=255)                                                                                 cACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGc                                                                     <  5:363281/148‑1 (MQ=255)                                                                                  aCGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCgg                                                                   <  3:265215/149‑1 (MQ=255)                                                                                  aCGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCgg                                                                   >  6:168987/1‑149 (MQ=255)                                                                                   cgccgcGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGc                                                                  <  3:9450/149‑1 (MQ=255)                                                                                    gccgcGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCt                                                                 >  4:124833/1‑149 (MQ=255)                                                                                         gcTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAg                                                            <  3:61178/149‑1 (MQ=255)                                                                                             cAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGc                                                          <  7:183798/147‑1 (MQ=255)                                                                                             cAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCg                                                         >  2:26358/1‑148 (MQ=255)                                                                                                   tACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACcg                                                  >  1:314424/1‑149 (MQ=255)                                                                                                     cacGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACcggc                                                >  4:30999/1‑149 (MQ=255)                                                                                                        gggTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGt                                             >  3:354365/1‑149 (MQ=255)                                                                                                         ggTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTg                                            <  6:197410/149‑1 (MQ=255)                                                                                                          gtgt‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTgc                                           <  5:149122/149‑1 (MQ=255)                                                                                                                 taCTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAAc                                     <  1:26358/149‑1 (MQ=255)                                                                                                                 taCTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAAc                                     >  2:304317/1‑149 (MQ=255)                                                                                                                    tCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGAt                                  >  7:43104/1‑149 (MQ=255)                                                                                                                        agagGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGAtcttc                              <  6:356319/149‑1 (MQ=255)                                                                                                                            ttggtCCTGGACGTCCATGCACTGGATGG‑CGGCAGACTTGTTGTCCCACACCTAGGTAATGTCCAGGAACGTGTCCGGTCGCCATTCGCACTGCGCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGAtcttctag                           <  2:291275/144‑3 (MQ=255)                                                                                                                              gtTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCTTCTCGccc                        >  3:204594/1‑149 (MQ=255)                                                                                                                               ttCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCTTCTCGCCCg                       >  5:290226/1‑149 (MQ=255)                                                                                                                               ttCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCTTCTCGCCCg                       >  3:111089/1‑149 (MQ=255)                                                                                                                                 ccTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCTTCTCGCCCgg                      <  4:204594/148‑1 (MQ=255)                                                                                                                                     gccggccATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCTTCTCGCCCGGCtt                   <  4:111089/147‑1 (MQ=255)                                                                                                                                      cccgccATCCAGTGCATGGCCGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCTTCTCGCCCGGCTTg                  >  2:39811/4‑148 (MQ=255)                                                                                                                                       ccgccATCCAGTGCATGGCCGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCTTCTCGCCCGGCTTGt                 <  1:39811/146‑1 (MQ=255)                                                                                                                                         gccATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCTTCTCGCCCGGCTTGTAGcc             <  3:30999/149‑1 (MQ=255)                                                                                                                                         gccATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCTTCTCGCCCGGCTTGTAGcc             <  8:43104/149‑1 (MQ=255)                                                                                                                                            aTGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCTTCTCGCCCGGCTTGAAGCCtt           >  1:136119/1‑148 (MQ=255)                                                                                                                                              gCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGAtcttc                              >  2:169213/1‑127 (MQ=255)                                                                                                                                              gCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGAtcttc                              <  1:169213/127‑1 (MQ=255)                                                                                                                                                aCTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCTTCTCGCCCGGCTTGTAGCCTTCAGc       >  5:301581/1‑148 (MQ=255)                                                                                                                                                aCTGGAGGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTTGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCTTCTCGCCCGGCTTGTAGCCTTCAGc       <  6:301581/148‑1 (MQ=255)                                                                                                                                                  tGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAg                                                            >  3:220685/1‑93 (MQ=255)                                                                                                                                                  tGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAg                                                            <  4:220685/93‑1 (MQ=255)                                                                                                                                                    gATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCTTCTCGCCCGGCTTGTAGCCTTCAGCCTgcg  <  3:124833/149‑1 (MQ=255)                                                                                                                                                    gATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCATCTCGCCCGGCTTGTAGCCTTCAGCCTgcg  >  4:138094/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      ACGCGTGGGAACACGCTCTGCAGGCCTTCGGCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCTTCTCGCCCGGCTTGTAGCCTTCAGCCTGCG  >  NC_002947/2862066‑2862360

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