Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,862,209 C→G 9.4% Q205H (CAG→CAC)  gllB ← 4‑oxalmesaconate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,862,2090CG 9.4% 188.5 / 8.7 64Q205H (CAG→CAC) gllB4‑oxalmesaconate hydratase
Reads supporting (aligned to +/- strand):  ref base C (27/31);  new base G (3/3);  total (30/34)
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.90e-01

ACGCGTGGGAACACGCTCTGCAGGCCTTCGGCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCTTCTCGCCCGGCTTGTAGCCTTCAGC  >  NC_002947/2862066‑2862355
                                                                                                                                                |                                                                                                                                                   
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                                                                                                            >  4:9450/1‑149 (MQ=255)
                                   gACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGt                                                                                                            >  5:356319/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                                                                                                    <  1:143682/148‑1 (MQ=255)
                                            gcgtgcgCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAAc                                                                                                    <  2:240119/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                                >  2:304317/1‑149 (MQ=255)
                                                                                                                taCTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAAc                                <  1:26358/149‑1 (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)
                                                                                                                                                |                                                                                                                                                   
ACGCGTGGGAACACGCTCTGCAGGCCTTCGGCGTAGACGATGTTGCGTGCGCTGGTTATGCCCACGTTACGCTTGGCCTGCACGCCGCGCTGCAGGGCTACACGGGTGT‑AGTACTCCCAGAGGTGTTCCTGGCCGGCCATGCACTGGATGG‑CGGCATACTTTTTGTCCCACACCGAGGTAATGTCCAGGAACGTGTCCGGGCGCCATTCGCACTGCTCGGGCTGGTGCGGCTCAAAGGCGTAGACCGGCGGTGCGCCAACGATCTTCTCGCCCGGCTTGTAGCCTTCAGC  >  NC_002947/2862066‑2862355

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