Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,736,589 Δ1 bp 6.2% coding (913/1674 nt) hutU ← urocanate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,736,5890C. 6.2% 310.4 / 14.4 97coding (913/1674 nt)hutUurocanate hydratase
Reads supporting (aligned to +/- strand):  ref base C (45/45);  new base . (3/3);  total (49/48)
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.19e-01

ACAGTGGGCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACCATGTCCGGGCGGAC  >  NC_002947/5736442‑5736736
                                                                                                                                                    |                                                                                                                                                   
aCAGTGGGCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACgg                                                                                                                                                   <  1:1991829/150‑1 (MQ=255)
    tGGGCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCAt                                                                                                                                               >  2:1099405/1‑150 (MQ=255)
       gCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGa                                                                                                                                            <  1:1902910/150‑1 (MQ=255)
         ggATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACt                                                                                                                                          <  2:888679/150‑1 (MQ=255)
         ggATGGAGGCCGGGACCAAGCCCGGGAAGTCGAAGGCAGTGGCCACGCCCTCCTCCTGGGCCATCTGCCGGATGTTTGTGCCGTAAGCGAAGCTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCGGGCACGTGCAGGCCAtgctgg                                                                                                                                          <  2:327720/150‑6 (MQ=255)
           aTGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGc                                                                                                                                        >  2:2376620/1‑150 (MQ=255)
            tGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCt                                                                                                                                       <  1:1984144/150‑1 (MQ=255)
                  ccGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCGTGCACGTGCAGGCCAtgctg                                                                                                                                           >  2:1258496/1‑136 (MQ=255)
                  ccGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCGTGCACGTGCAGGCCAtgctg                                                                                                                                           <  1:1258496/140‑5 (MQ=255)
                    gggACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGcc                                                                                                                               >  1:2094975/1‑150 (MQ=255)
                         gAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGAc                                                                                                                          >  1:2312911/1‑150 (MQ=255)
                              cTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTg                                                                                                                     >  1:1040580/1‑150 (MQ=255)
                                   aaGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGcc                                                                                                                <  2:1624614/150‑1 (MQ=255)
                                      tCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGAc                                                                                                                          >  1:302011/1‑137 (MQ=255)
                                      tCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGAc                                                                                                                          <  2:302011/137‑1 (MQ=255)
                                         aaGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGcc                                                                                                                               <  1:672754/129‑1 (MQ=255)
                                         aaGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGcc                                                                                                                               >  2:672754/1‑129 (MQ=255)
                                          aGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCa                                                                                                                   >  2:848143/1‑140 (MQ=255)
                                          aGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCa                                                                                                                   <  1:848143/140‑1 (MQ=255)
                                           ggCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGt                                                                                                        <  1:1628896/150‑1 (MQ=255)
                                              aTTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAg                                                                                                                  >  2:2237628/1‑137 (MQ=255)
                                              aTTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAg                                                                                                                  <  1:2237628/137‑1 (MQ=255)
                                              aTTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGGCCGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGAAAGGCAAACCACTGCTTGGCGGCGTAGGGCACTCCAGTCGGTTCGGTCTg                                                                                                     >  2:22057/1‑150 (MQ=255)
                                               ttGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTgc                                                                                                    >  2:2247520/1‑150 (MQ=255)
                                               ttGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTgc                                                                                                    <  1:2376620/150‑1 (MQ=255)
                                                    cACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACg                                                                                               >  1:1560773/1‑150 (MQ=255)
                                                          ctcctcctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGAc                                                                                                                          <  2:778499/117‑1 (MQ=255)
                                                          ctcctcctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGAc                                                                                                                          >  1:778499/1‑117 (MQ=255)
                                                           tcctcctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAACATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCgg                                                                                        <  1:761820/150‑1 (MQ=255)
                                                             ctcctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTa                                                                                      >  1:667859/1‑150 (MQ=255)
                                                              tcctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTAc                                                                                     <  1:2247520/150‑1 (MQ=255)
                                                               cctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACt                                                                                    >  1:1015642/1‑150 (MQ=255)
                                                                ctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCGTGCACGTGCAGGCCAtg                                                                                                                                              <  1:775324/91‑2 (MQ=39)
                                                                ctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCGTGCACGTGCAGGCCAtg                                                                                                                                              >  2:775324/1‑90 (MQ=39)
                                                                ctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGGGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTg                                                                                   >  2:1613103/1‑150 (MQ=255)
                                                                ctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCgg                                                                                        >  2:342793/1‑145 (MQ=255)
                                                                ctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCgg                                                                                        <  1:342793/145‑1 (MQ=255)
                                                                ctTGGCCATCTGGCGGATGTTGTTGCCGTAATAGAAGGTCGGGATGCACTGCTTATGGAAGTCCAACATGACCTGCACGTGCACGGCCATCGACTACATGGCGGCCTTGACCACTGCAGCCGGTACGGTCTGCGCACGATAGCGGTACTg                                                                                   >  2:992224/1‑150 (MQ=255)
                                                                    gCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATTGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTcc                                                                               <  1:629348/150‑1 (MQ=255)
                                                                       aTCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAg                                                                            <  2:2312911/150‑1 (MQ=255)
                                                                       aTCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAg                                                                            <  2:667859/150‑1 (MQ=255)
                                                                       aTCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATCTCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAg                                                                            <  1:22057/150‑1 (MQ=255)
                                                                          tGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATCCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTc                                                                         <  2:1291756/150‑1 (MQ=255)
                                                                            gCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCa                                                                             <  2:932582/144‑1 (MQ=255)
                                                                            gCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCa                                                                             >  1:932582/1‑144 (MQ=255)
                                                                              ggATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGc                                                                     <  1:1285403/150‑1 (MQ=255)
                                                                                      ttgCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGc                                                             <  1:1230196/150‑1 (MQ=255)
                                                                                       tgCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCa                                                            <  2:1015642/150‑1 (MQ=255)
                                                                                          cGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGt                                                         >  1:1585154/1‑150 (MQ=255)
                                                                                              aTCGAAGGTCGGGATGTCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGcc                                                     <  2:613386/150‑1 (MQ=255)
                                                                                               tCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTggcg                                                                                                                                  >  1:1579024/1‑72 (MQ=255)
                                                                                               tCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTggcg                                                                                                                                  <  2:1579024/72‑1 (MQ=255)
                                                                                               tCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCg                                                    >  1:485058/1‑150 (MQ=255)
                                                                                                cGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAg                                                                            >  2:2171116/1‑125 (MQ=255)
                                                                                                cGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAg                                                                            <  1:2171116/125‑1 (MQ=255)
                                                                                                 gAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGtt                                                  <  2:1585154/150‑1 (MQ=255)
                                                                                                 gAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGtt                                                  >  1:771062/1‑150 (MQ=255)
                                                                                                        gggATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTggg                                           <  2:2199295/150‑1 (MQ=255)
                                                                                                        gggATACCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTggg                                           >  2:2058077/1‑150 (MQ=255)
                                                                                                         ggATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCATGTGCACGGCCATc                                                                                                                                              <  1:2394011/50‑1 (MQ=255)
                                                                                                         ggATGCCCTGATTCTGGAAGTCCAGCATGGCCTGCATGTGCACGGCCATc                                                                                                                                              >  2:2394011/1‑50 (MQ=38)
                                                                                                            tGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTcc                                                                               <  2:2158730/110‑1 (MQ=255)
                                                                                                            tGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTcc                                                                               >  1:2158730/1‑110 (MQ=255)
                                                                                                            tGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGt                                       <  2:2094975/150‑1 (MQ=255)
                                                                                                              cccTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGtt                                                                                 <  2:1711403/106‑1 (MQ=255)
                                                                                                              cccTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGtt                                                                                 >  1:1711403/1‑106 (MQ=255)
                                                                                                               ccTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCt                                                                <  1:665431/122‑1 (MQ=255)
                                                                                                               ccTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCt                                                                >  2:665431/1‑122 (MQ=255)
                                                                                                               ccTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGt                                          >  2:80326/1‑144 (MQ=255)
                                                                                                               ccTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTggg                                    >  2:1777248/1‑150 (MQ=255)
                                                                                                               ccTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTggg                                    >  2:623934/1‑150 (MQ=255)
                                                                                                               ccTGCTTCTGGAAGCCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGt                                          <  1:80326/144‑1 (MQ=255)
                                                                                                                cTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGgc                                   >  2:1913121/1‑150 (MQ=255)
                                                                                                                cTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGgc                                   >  2:376755/1‑150 (MQ=255)
                                                                                                                 tGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGgcg                                  >  2:1685913/1‑150 (MQ=255)
                                                                                                                   cTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGct                                >  2:1139286/1‑150 (MQ=255)
                                                                                                                       tGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGcc                                                     <  2:192713/125‑1 (MQ=255)
                                                                                                                       tGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGcc                                                     >  1:192713/1‑125 (MQ=255)
                                                                                                                       tGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCa                                                >  1:1085717/1‑130 (MQ=255)
                                                                                                                       tGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCa                                                <  2:1085717/130‑1 (MQ=255)
                                                                                                                         gAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGAc                                                                                                                          >  2:326176/1‑54 (MQ=255)
                                                                                                                         gAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGAc                                                                                                                          <  1:326176/54‑1 (MQ=255)
                                                                                                                            gTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCtggtctggtc                       <  1:1139286/150‑1 (MQ=255)
                                                                                                                             tCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTg                                      <  1:1681669/134‑1 (MQ=255)
                                                                                                                             tCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTg                                      >  2:1681669/1‑134 (MQ=255)
                                                                                                                               cAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATc                                                                                            >  1:1000589/1‑78 (MQ=255)
                                                                                                                               cAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATc                                                                                            <  2:1000589/78‑1 (MQ=255)
                                                                                                                               cAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTggg                                           >  2:1406490/1‑127 (MQ=255)
                                                                                                                               cAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTggg                                           <  1:1406490/127‑1 (MQ=255)
                                                                                                                               cAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTCCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCtggtctggtcacc                    >  1:1897244/1‑147 (MQ=255)
                                                                                                                                tcgaTGGCCGTGCACGTGCAGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGTTACTGTTCCCAGGTCAAGCCGGATGGACTGAAGCCCGTCATTTGGTCGTGGGCGCTGGTCTGGTCGGTg                   >  1:2340172/4‑150 (MQ=255)
                                                                                                                                  cATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTgg                                            >  2:1577855/1‑123 (MQ=255)
                                                                                                                                  cATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTgg                                            <  1:1577855/123‑1 (MQ=255)
                                                                                                                                  cATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTg                                      >  1:1980360/1‑129 (MQ=255)
                                                                                                                                  cATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTg                                      <  2:1980360/129‑1 (MQ=255)
                                                                                                                                     ggCCTGCACGTGCACGCTCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACCAt              <  2:751529/150‑1 (MQ=255)
                                                                                                                                        cTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACCATGTc           >  1:829199/1‑150 (MQ=255)
                                                                                                                                            cacgtgcacgGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGCCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACCATGTCCgg        >  1:1221253/1‑150 (MQ=255)
                                                                                                                                                  cacgGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACCATGTCCGGGCGGAc  <  1:623934/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACAGTGGGCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACCATGTCCGGGCGGAC  >  NC_002947/5736442‑5736736

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: