Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,861,136 Δ1 bp 8.6% coding (196/1086 nt) gllD ← 4‑oxalomesaconate tautomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,861,1360G. 8.6% 178.9 / 11.1 58coding (196/1086 nt)gllD4‑oxalomesaconate tautomerase
Reads supporting (aligned to +/- strand):  ref base G (27/26);  new base . (3/2);  total (30/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 = 6.06e-01

CCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATAGGCGCCCTTG  >  NC_002947/2860995‑2861284
                                                                                                                                              |                                                                                                                                                    
ccaccaGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGatgatg                                                                                                                                               >  6:116227/1‑149 (MQ=255)
ccaccaGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGatgatg                                                                                                                                               <  4:327327/149‑1 (MQ=255)
  accaGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGc                                                                                                                                             >  6:235460/1‑149 (MQ=255)
  accaGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGc                                                                                                                                             <  3:289432/149‑1 (MQ=255)
   ccaGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCg                                                                                                                                            <  3:109560/149‑1 (MQ=255)
    caGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCAGCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGACGCGCGGATGATGGc                                                                                                                                             >  7:299505/1‑147 (MQ=255)
    caGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGa                                                                                                                                           <  5:116227/149‑1 (MQ=255)
       gCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACtt                                                                                                                                        >  5:23859/1‑149 (MQ=255)
       gCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGCCATCATCCGCGCTTCGCAGCGCGATGACGCCGAcgt                                                                                                                                        <  6:318207/149‑3 (MQ=37)
         cgcgCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTg                                                                                                                                      <  8:47378/149‑1 (MQ=255)
         cgcgCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTg                                                                                                                                      <  5:235460/149‑1 (MQ=255)
              tCAAGGGCGAACGGCCCTATGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGt                                                                                                                                 >  1:21993/1‑149 (MQ=255)
              tCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTAAGCGAACAGGTAGTCCACGTCGGCGTCAACGCGATGCCAACCACGGATAAAGGCAAAATTGCTgg                                                                                                                                  >  8:230034/1‑148 (MQ=255)
                aaGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGt                                                                                                                                 <  7:274661/147‑1 (MQ=255)
                  gggCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGc                                                                                                                             <  6:23859/149‑1 (MQ=255)
                   ggCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGc                                                                                                                             >  6:109211/1‑148 (MQ=255)
                    gCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGa                                                                                                                           <  3:75957/149‑1 (MQ=255)
                       aaCGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTc                                                                                                                        <  7:246035/149‑1 (MQ=255)
                             ccTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGcgcc                                                                                                                  >  7:103191/1‑149 (MQ=255)
                               tACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGcgccgc                                                                                                                >  3:186503/1‑149 (MQ=255)
                                aCGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCgccgcc                                                                                                               >  7:165513/1‑149 (MQ=255)
                                 cGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCgccgcc                                                                                                               <  2:18007/148‑1 (MQ=255)
                                   ccGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGAt                                                                                                            <  4:186503/149‑1 (MQ=255)
                                                   cgttgTTTTGCCCTTTGTCCACCCGCGCCTGGTCCACCACCACCTGAGCTTACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCg                                                                                             <  7:230034/144‑1 (MQ=255)
                                                            gCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGc                                                                                                                                     >  5:227454/1‑99 (MQ=255)
                                                            gCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGc                                                                                                                                     <  6:227454/99‑1 (MQ=255)
                                                             cccATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGc                                                                                   <  4:75950/148‑1 (MQ=255)
                                                             cccATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCg                                                                                  >  4:112424/1‑149 (MQ=255)
                                                              ccATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGa                                                                                 <  2:73519/149‑1 (MQ=255)
                                                                  aGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGcc                                                                              >  5:92894/1‑148 (MQ=255)
                                                                      cACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATc                                                                          <  2:21993/148‑1 (MQ=255)
                                                                             gcCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCcag                                                                  <  1:332266/149‑1 (MQ=255)
                                                                               cTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATCCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGccacca                                                                >  5:345287/1‑149 (MQ=255)
                                                                                tCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCcagcag                                                               >  7:195084/1‑149 (MQ=255)
                                                                                  gTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCcagcag                                                               <  8:165513/147‑1 (MQ=255)
                                                                                      accaccaccTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCg                                                         >  7:265223/1‑149 (MQ=255)
                                                                                        caccaccTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGt                                                       >  4:116769/1‑149 (MQ=255)
                                                                                           caccTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCAc                                                    >  7:214483/1‑149 (MQ=255)
                                                                                                   cGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCGCGCAGCggg                                            <  6:246719/149‑1 (MQ=255)
                                                                                                      aCAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGc                                           >  2:243587/1‑147 (MQ=255)
                                                                                                       cAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAg                                        <  8:195084/149‑1 (MQ=255)
                                                                                                          gTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAgggg                                     >  4:39012/1‑149 (MQ=255)
                                                                                                          gTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAgggg                                     >  7:74647/1‑149 (MQ=255)
                                                                                                               ccACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGgggcgggc                                <  3:39012/149‑1 (MQ=255)
                                                                                                                cACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGgggcgggc                                <  8:265223/148‑1 (MQ=255)
                                                                                                                cACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGgggcggg                                 >  5:191156/1‑147 (MQ=255)
                                                                                                                cACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGgggcggg                                 >  4:96608/1‑147 (MQ=255)
                                                                                                                    tCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGgtc                           <  1:243587/149‑1 (MQ=255)
                                                                                                                    tCGCCGTCATCCGCGCTTCGCAGCGCGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGgtc                           >  7:336135/1‑149 (MQ=255)
                                                                                                                      ggCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGgtcgt                         <  7:181635/149‑1 (MQ=255)
                                                                                                                        catcatccgcgCTTCGCAGCGCGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCa                       >  1:241131/8‑149 (MQ=255)
                                                                                                                        catcatccgcgCTTCGCAGCGCGATGACGCCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCCCCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCa                       >  7:153984/8‑149 (MQ=255)
                                                                                                                            atccgcgCTTCGCAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCTATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATg                     >  3:84913/4‑148 (MQ=255)
                                                                                                                                  gcTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCaaaaaa              >  7:142103/1‑149 (MQ=255)
                                                                                                                                  gcTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCaaaaaa              >  7:48859/1‑149 (MQ=255)
                                                                                                                                   cttcgcaGCGCGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATa            <  8:336135/143‑1 (MQ=255)
                                                                                                                                      cGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATAg           <  8:103191/148‑1 (MQ=255)
                                                                                                                                        aaGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATAGgcg        >  1:123399/1‑149 (MQ=255)
                                                                                                                                              ggATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATAGGCGCCCTTg  <  3:96608/149‑1 (MQ=255)
                                                                                                                                              ggATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCAGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATAGGCGCCCTTg  >  1:95956/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CCACCAGGCCGCGCTCAAGGGCGAACGGCCCTACGCCGGCCAGAATGTTGCCGCAGTTTTGCCCATAGTCCACCCGCGCCTCGTCCACCACCACCTGAGCGAACAGGTAGTCCACGTCGGCGTCATCGCGCTGCGAAGCGCGGATGATGGCGACTTTGCTGGTGAGCGAGTCGGCGCCGCCGATGCCGTCAATCTGCCGGGCATCGGGCGAGCCCATCACCGCCAGCAGCACCCGGTCACGCAGCGGGCCAGGGGCGGGCAGGTCGTCATGCAAAAAATAGGCGCCCTTG  >  NC_002947/2860995‑2861284

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: