Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,862,206	G→C	10.0%	C206W (TGC→TGG)	gllB ←	4‑oxalmesaconate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,862,206	0	G	C	10.0%	180.2 / 8.6	60	C206W (TGC→TGG)	gllB	4‑oxalmesaconate hydratase
Reads supporting (aligned to +/- strand):  ref base G (25/29);  new base C (3/3);  total (28/32)
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.67e-01

CTCCGTCACGCGTGGGAACACGCTCTGCAGGCCTTCGGCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCTTCTCGCCCGGCTTGTAGCCTT  >  NC_002947/2862059‑2862351                                                                                                                                                     |                                                                                                                                                   ctcCGTCACGCGTGGGAACACGCTCTGCAGGCCTTCGGCGTAGACGATGTTTCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATg                                                                                                                                                    <  1:246745/148‑1 (MQ=255) ctcCGTCACGCGTGGGAACACGCTCTGCAGGCCTTCGGCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATg                                                                                                                                                    >  5:197410/1‑148 (MQ=255)        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                                                          >  6:168987/1‑149 (MQ=255)                                                                                         aCGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCgg                                                          <  3:265215/149‑1 (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)                                                                                                                                                     |                                                                                                                                                   CTCCGTCACGCGTGGGAACACGCTCTGCAGGCCTTCGGCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCTTCTCGCCCGGCTTGTAGCCTT  >  NC_002947/2862059‑2862351

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