Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,993,997 C→G 7.5% T227T (ACC→ACG)  bkdAB → branched‑chain alpha‑keto acid dehydrogenase complex subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,993,9970CG 7.5% 198.1 / 3.2 64T227T (ACC→ACG) bkdABbranched‑chain alpha‑keto acid dehydrogenase complex subunit beta
Reads supporting (aligned to +/- strand):  ref base C (31/28);  new base G (3/2);  total (34/30)
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.85e-01

TGTACAACGGCCCGTTCGATGGCCACCACGACCGCCCTGTAACCCCGTGGTCGAAGCACCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTATCGTCGAGTCGGTGAAAAAGACCGGCCGTTGCGTGGTGGTGCA  >  NC_002947/4993849‑4994143
                                                                                                                                                    |                                                                                                                                                  
tGTACAACGGCCCGTTCGATGGCCACCACGACCGCCCTGTAACCCCGTGGTCGAAGCACCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGAcc                                                                                                                                                    >  3:308053/1‑149 (MQ=255)
     aaCGGCCCGTTCGATGGCCACCACGACCGCCCTGTAACCCCGTGGTCGAAGCACCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTCCCCGCCCTGGCAATGACGTGACCGTGCTGACCTACgg                                                                                                                                               >  2:227988/1‑149 (MQ=255)
            cGTTCGATGGCCACCACGACCGCCCTGTAACCCCGTGGTCGAAGCACCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACTTGACCGTGCTGACCTACGGCACCACg                                                                                                                                        <  8:22615/149‑1 (MQ=255)
                 gATGGCCACCACGACCGCCCTGTAACCCCGTGGTCGAAGCACCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTa                                                                                                                                   <  5:65158/149‑1 (MQ=255)
                   tGGCCACCACGACCTCCCTGTAACCCCGTGGTAGACGCACCCGCACAGCGCCGAGCCCGACGGAGATTACACCGTACCGCTGGAACAGGCCGCCATAACCCGCCCTGGCAATGACGTGACCGGGCGGACCTACGGCCCCACGGGGTCCg                                                                                                                                 >  8:145038/1‑149 (MQ=255)
                         ccacGACCGCCCTGTAACCCCGTGGTCGAAGCACCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGcc                                                                                                                            >  2:76962/1‑148 (MQ=255)
                           acGACCGCCCTGTAACCCCGTGGTCGAAGCACCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAg                                                                                                                         <  2:116349/149‑1 (MQ=255)
                                cGCCCTGTAACCCCGTGGTCGAAGCACCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGgc                                                                                                                    <  8:78327/149‑1 (MQ=255)
                                         aCCCCGTGGTCGAAGCACCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCAGCCgatg                                                                                                            >  6:125451/1‑146 (MQ=255)
                                            ccGTGGTCGAAGCACCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAg                                                                                                        >  1:281372/1‑149 (MQ=255)
                                            ccGTGGTCGAAGCACCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGCCCCAGGTGGCCGCCGAAGAAAg                                                                                                        >  8:231463/1‑149 (MQ=255)
                                            ccGTGGTCGAAGCACCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCATGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAg                                                                                                        >  6:187501/1‑149 (MQ=255)
                                             cGTGGTCGAAGCACCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAgc                                                                                                       >  2:206701/1‑149 (MQ=255)
                                               tGGTCGAAGCACCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCATGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAgcgg                                                                                                     <  5:187501/149‑1 (MQ=255)
                                               tGGGCGAACCACCCGCCCAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAgcgg                                                                                                     <  4:246195/149‑1 (MQ=255)
                                                   cGAAGCACCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACCGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTc                                                                                                 <  1:188648/149‑1 (MQ=255)
                                                    gAAGCACCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCGGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTc                                                                                                 <  5:120718/148‑1 (MQ=255)
                                                      aGCACCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGAt                                                                                              <  3:168818/149‑1 (MQ=255)
                                                         aCCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGcc                                                                                           >  2:92347/1‑149 (MQ=255)
                                                         aCCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGcc                                                                                           >  4:172842/1‑149 (MQ=255)
                                                          cccGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCAATACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCg                                                                                          <  2:135651/149‑1 (MQ=255)
                                                            cGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGaa                                                                                        <  1:76962/149‑1 (MQ=255)
                                                             gCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAg                                                                                       <  2:281372/149‑1 (MQ=255)
                                                                 aGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGAt                                                                                   >  1:216328/1‑149 (MQ=255)
                                                                     ccGTGCCCGACGGGTATTACAGCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGAc                                                                               <  3:18425/149‑1 (MQ=255)
                                                                         gCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTgc                                                                           >  8:37172/1‑149 (MQ=255)
                                                                           ccGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCt                                                                             <  5:52848/145‑1 (MQ=255)
                                                                           ccGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCt                                                                             >  6:52848/1‑145 (MQ=255)
                                                                            cGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCa                                                                        <  7:231463/149‑1 (MQ=255)
                                                                                     ttaCACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGc                                                               >  4:206924/1‑149 (MQ=255)
                                                                                      taCACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGcc                                                              >  7:102825/1‑149 (MQ=255)
                                                                                      taCACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGcc                                                              >  6:313455/1‑149 (MQ=255)
                                                                                      taCACCGTACCGATGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGAACGTGCTGACCTACGGCACCACGGAGTAAGTGGCCCAGGTGGCCGCCCAAGAAAGCGGCGTAGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGcc                                                              >  3:253507/1‑149 (MQ=255)
                                                                                        cacCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGc                                                            >  6:295227/1‑149 (MQ=255)
                                                                                          cgggaCCGCTGGACCAGGCCGCCATTACCCGCCCTGGCAATGACGTGCCCGTGCTGCAGTAGGGCTCCACGGTGGACGTGGCCCAGGTGGCCTCCGAAGAAAGCGGCGTCGATGCCGATGTGATCGAGCTGCGCAGCCTGTGGCCGctg                                                          <  7:145038/145‑1 (MQ=255)
                                                                                           cGTACCGCTGGACAAGGCCGCCATTCCCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGctgg                                                         <  1:227988/149‑1 (MQ=255)
                                                                                                 gCTGGACAAGGCCGCCATTACACGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGctggacctg                                                    >  6:316847/1‑148 (MQ=255)
                                                                                                  cTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGctggacctgga                                                  <  3:172842/149‑1 (MQ=255)
                                                                                                   tGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGctggac                                                       <  2:250267/143‑1 (MQ=255)
                                                                                                   tGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGctggac                                                       <  4:101581/143‑1 (MQ=255)
                                                                                                   tGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGctggac                                                       >  3:101581/1‑143 (MQ=255)
                                                                                                   tGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGctggac                                                       >  1:250267/1‑143 (MQ=255)
                                                                                                      aCAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTATGGCCGCTGGACCTGGacac                                               >  3:180994/1‑148 (MQ=255)
                                                                                                        aaGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCCGCCTGTGGCCGCTGGACCTGGACACttt                                            >  8:53813/1‑147 (MQ=255)
                                                                                                                cATTACCCGCCCTGGCAATGACGTGACCGTGCGGACCACCGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGCCGTCGATGCCGAAGTGATCGACCTGCGCAGCCGGTGGCCGCTGGACCTGGACACTATCGTCGAGt                                    <  6:328978/149‑1 (MQ=255)
                                                                                                                   tACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTATCGTCGAGTCgg                                 >  7:236144/1‑149 (MQ=255)
                                                                                                                    aCCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTATCGTCGAGTCGGt                                >  4:21366/1‑149 (MQ=255)
                                                                                                                     cccGCCCTGGCTATGACGTTACGGTGCTGACGTACGGCACCATGGTGTTCGTGGCCGAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCGGAAGTGATCGACCTGCGCAGACTGTGGCCGCTGGACCTGGACACTATCGTAGAGTCGGt                                <  4:253507/148‑1 (MQ=255)
                                                                                                                        gCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTATCGTCGAGTCGGTGaa                             >  6:75914/1‑148 (MQ=255)
                                                                                                                              gCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTATCGTCGAGTCGGTGAAAAAGAcc                      <  2:216328/149‑1 (MQ=255)
                                                                                                                               cAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTATCGTCGAGTCGGTGAAAAAGACCg                     <  1:206701/149‑1 (MQ=255)
                                                                                                                                   gACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTATCGTCGAGTCGGTGAAAAAGACCGGCCg                 >  7:29302/1‑149 (MQ=255)
                                                                                                                                     cGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTATCGTCGAGTCGGTGAAAAAGACCGGCCGtt               <  1:12406/149‑1 (MQ=255)
                                                                                                                                      gTGGCCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTATCGTCGAGTCGGTGAAAAAGa                        >  6:160315/1‑139 (MQ=255)
                                                                                                                                      gTGGCCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTATCGTCGAGTCGGTGAAAAAGa                        <  5:160315/139‑1 (MQ=255)
                                                                                                                                       acaCCGTGGTGCCGTAGGTCAGCACGGTGTACGTAGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTATCGTCGAGTCGGTGAAAAAGACCGGCCGTTGc             >  7:142387/3‑149 (MQ=255)
                                                                                                                                         aCCGTGGTGCCGTAGGTCAGCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTATCGTCGAGTCGGTGAAAAAGACCGGCCGTTGCgt           >  1:63610/1‑149 (MQ=255)
                                                                                                                                         aCCGTGGTGCCGTAGGTCAGCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTATCGTCGAGTCGGGGAAAAAGACCGGCCGTTGCgt           >  8:273045/1‑149 (MQ=255)
                                                                                                                                             tGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTATCGTCGAGTCGGTGAAAAAGACCGGCCGTTGCgtggtg       <  5:316847/149‑1 (MQ=255)
                                                                                                                                               cTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTATCGTCGAGTCGGTGAAAAAGACCGGCCGTTGCgtggtggt     <  3:21366/149‑1 (MQ=255)
                                                                                                                                                tGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTAtcgtc                                        >  5:2616/1‑113 (MQ=255)
                                                                                                                                                tGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTAtcgtc                                        <  6:2616/113‑1 (MQ=255)
                                                                                                                                                 gACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTATCGTCGAGTCGGTGAAAAAGACCGGCCGTTGCgt           >  4:256879/1‑141 (MQ=255)
                                                                                                                                                 gACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTATCGTCGAGTCGGTGAAAAAGACCGGCCGTTGCgt           <  3:256879/141‑1 (MQ=255)
                                                                                                                                                  aCCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTATCGTCGAGTCGGTGAAAAAGACCGGCCGTTGCGTGGTGGTGCa  >  3:88468/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGTACAACGGCCCGTTCGATGGCCACCACGACCGCCCTGTAACCCCGTGGTCGAAGCACCCGCACAGCGCCGTGCCCGACGGTTATTACACCGTACCGCTGGACAAGGCCGCCATTACCCGCCCTGGCAATGACGTGACCGTGCTGACCTACGGCACCACGGTGTACGTGGCCCAGGTGGCCGCCGAAGAAAGCGGCGTCGATGCCGAAGTGATCGACCTGCGCAGCCTGTGGCCGCTGGACCTGGACACTATCGTCGAGTCGGTGAAAAAGACCGGCCGTTGCGTGGTGGTGCA  >  NC_002947/4993849‑4994143

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