Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,737,998 G→A 5.4% L476L (CTG→TTG)  hutU ← urocanate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,737,9980GA 5.4% 462.8 / 10.1 147L476L (CTG→TTG) hutUurocanate hydratase
Reads supporting (aligned to +/- strand):  ref base G (70/69);  new base A (5/3);  total (75/72)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 162 total aligned reads displayed.

GTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGT  >  NC_002947/5737849‑5738147
                                                                                                                                                     |                                                                                                                                                     
gTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTcagcag                                                                                                                                                       <  2:4053406/150‑1 (MQ=255)
gTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTcagcag                                                                                                                                                       >  2:3429289/1‑150 (MQ=255)
 tCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAgg                                                                                                                                                      <  2:3294865/150‑1 (MQ=255)
   gTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGc                                                                                                                                                    >  1:3774241/1‑150 (MQ=255)
    ttGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCg                                                                                                                                                   >  2:3099047/1‑150 (MQ=255)
       gTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTg                                                                                                                                                >  2:364151/1‑150 (MQ=255)
        tCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGa                                                                                                                                               >  2:1155687/1‑150 (MQ=255)
          aGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGagc                                                                                                                                             >  1:4949153/1‑150 (MQ=255)
             aCCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGagcagc                                                                                                                                          <  2:5280665/150‑1 (MQ=255)
             aCCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGagcagc                                                                                                                                          <  2:4382334/150‑1 (MQ=255)
                     gATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTc                                                                                                                                  >  1:886991/1‑150 (MQ=255)
                     gATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTc                                                                                                                                  <  1:4042639/150‑1 (MQ=255)
                        gcgcTGGTCGGCATCTTCCGAACCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGa                                                                                                                               <  2:571286/150‑1 (MQ=255)
                        gcgcTGGTCGGCATCTTCCGAACCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGa                                                                                                                               <  1:3038833/150‑1 (MQ=255)
                         cgcTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTGGAGCAGCGGCCAGTCGGac                                                                                                                              >  2:4556017/1‑149 (MQ=255)
                          gcTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGaaa                                                                                                                             >  2:843330/1‑150 (MQ=255)
                          gcTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGaaa                                                                                                                             >  1:2839646/1‑150 (MQ=255)
                           cTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAAc                                                                                                                            <  1:1201244/150‑1 (MQ=255)
                            tCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACa                                                                                                                           >  1:773737/1‑150 (MQ=255)
                              gtcggcATCTTCCGAACCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGc                                                                                                                         >  1:3789756/3‑150 (MQ=255)
                              ggcggcCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGc                                                                                                                         <  1:3616891/150‑1 (MQ=255)
                                cggcCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGt                                                                                                                       <  1:4238729/150‑1 (MQ=255)
                                cggcCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGt                                                                                                                       <  1:4238653/150‑1 (MQ=255)
                                 ggcCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTc                                                                                                                      <  1:4214664/150‑1 (MQ=255)
                                  gcCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTcc                                                                                                                     <  2:4819430/150‑1 (MQ=255)
                                     tCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGaa                                                                                                                  >  2:3338062/1‑150 (MQ=255)
                                     tCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGaa                                                                                                                  >  2:4240856/1‑150 (MQ=255)
                                       ttccgaaCCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGcc                                                                                                                >  1:5306458/7‑150 (MQ=255)
                                        tCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCa                                                                                                               >  2:2118177/1‑150 (MQ=255)
                                        tCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCa                                                                                                               <  2:1777024/150‑1 (MQ=255)
                                          ggTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATc                                                                                                             <  2:5375208/150‑1 (MQ=255)
                                           gaaCCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGTAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGCCGTct                                                                                                            >  2:4666671/3‑149 (MQ=255)
                                             aCCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACg                                                                                                          >  2:2426772/1‑150 (MQ=255)
                                               cGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCa                                                                                                        <  1:2084740/150‑1 (MQ=255)
                                                  cGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCg                                                                                                                                                   >  1:3915655/1‑104 (MQ=255)
                                                  cGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCg                                                                                                                                                   <  2:3915655/104‑1 (MQ=255)
                                                       cagATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGCCGTcttt                                                                                                          <  2:2569331/138‑4 (MQ=255)
                                                              ccATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGCCGTCTTTCATGGATTCGGTTTcac                                                                                         >  1:3207414/1‑148 (MQ=255)
                                                               cATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcgcg                                                                                        <  2:3576568/150‑1 (MQ=255)
                                                               cATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcgcg                                                                                        >  1:3228503/1‑150 (MQ=255)
                                                               cATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcgcg                                                                                        >  1:3228499/1‑150 (MQ=255)
                                                               cATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcgcg                                                                                        <  2:1065557/150‑1 (MQ=255)
                                                               cATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcgcg                                                                                        <  2:3218485/150‑1 (MQ=255)
                                                                  gCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGtt                                                                                     <  1:1278841/150‑1 (MQ=255)
                                                                  gCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGtt                                                                                     <  1:2193900/150‑1 (MQ=255)
                                                                    ccGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCg                                                                                   <  1:2525114/150‑1 (MQ=255)
                                                                    ccGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCg                                                                                   >  1:4386620/1‑150 (MQ=255)
                                                                     cGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCgg                                                                                  >  2:2277360/1‑150 (MQ=255)
                                                                      gAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCggg                                                                                 <  1:1161176/150‑1 (MQ=255)
                                                                            tgcGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGAt                                                                           >  1:4808556/1‑150 (MQ=255)
                                                                            tgcGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGAt                                                                           >  1:2082210/1‑150 (MQ=255)
                                                                              cGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATAc                                                                         <  1:2537654/150‑1 (MQ=255)
                                                                              cGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATAc                                                                         <  1:2537486/150‑1 (MQ=255)
                                                                                agaAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCg                                                                       >  1:2572102/1‑150 (MQ=255)
                                                                                agaAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCg                                                                       >  1:4722035/1‑150 (MQ=255)
                                                                                    gCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGcc                                                                   <  1:4642104/150‑1 (MQ=255)
                                                                                    gCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGcc                                                                   >  2:4566559/1‑150 (MQ=255)
                                                                                      ccATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCgg                                                                 >  2:1939091/1‑150 (MQ=255)
                                                                                      ccATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCgg                                                                 >  1:2047648/1‑150 (MQ=255)
                                                                                       cATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGa                                                                >  1:2704778/1‑150 (MQ=255)
                                                                                        aTGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAg                                                               >  2:3063647/1‑150 (MQ=255)
                                                                                            ccACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGa                                                           <  2:3989838/150‑1 (MQ=255)
                                                                                            acACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGTGCAGCGGCCAGTCGGACACAGCGTCCGAACCATCCCGCATGGCCTCGGGTTCGGGGTTCGGGCGGGATTGCGTGCCGGTGTCgg                                                           >  2:3962883/2‑149 (MQ=255)
                                                                                             cACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAg                                                          <  2:4722866/150‑1 (MQ=255)
                                                                                              aCGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAgg                                                         >  2:3708379/1‑150 (MQ=255)
                                                                                               cGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGTGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAggt                                                        <  1:412726/150‑1 (MQ=255)
                                                                                                gCCACCGCCATGGTGCAGCGATTCCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGTGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAggtg                                                       <  2:4284498/150‑1 (MQ=255)
                                                                                                gCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAggtg                                                       >  2:4878649/1‑150 (MQ=255)
                                                                                                gCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAggtg                                                       <  1:1754748/150‑1 (MQ=255)
                                                                                                gCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAggtg                                                       >  2:4878646/1‑150 (MQ=255)
                                                                                                  cACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAggtggt                                                     >  1:336055/1‑150 (MQ=255)
                                                                                                    ccGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCa                                                   >  2:2448987/1‑150 (MQ=255)
                                                                                                      gCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGGAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACg                                                 >  1:1040255/1‑150 (MQ=255)
                                                                                                       ccATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACgg                                                <  2:3406473/150‑1 (MQ=255)
                                                                                                              gCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATca                                         <  2:5344781/150‑1 (MQ=255)
                                                                                                               cAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATcac                                        <  2:5118700/150‑1 (MQ=255)
                                                                                                                aGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcg                                                                                          <  2:1670095/99‑1 (MQ=255)
                                                                                                                aGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcg                                                                                          >  1:1670095/1‑99 (MQ=255)
                                                                                                                aGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATcacc                                       <  2:2557749/150‑1 (MQ=255)
                                                                                                                        ccAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGc                               >  1:1460337/1‑150 (MQ=255)
                                                                                                                        ccAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGc                               <  1:5208431/150‑1 (MQ=255)
                                                                                                                        ccAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGc                               >  2:802018/1‑150 (MQ=255)
                                                                                                                         taGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGcc                              >  1:3893343/2‑150 (MQ=255)
                                                                                                                           ggTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGa                            <  2:3258607/150‑1 (MQ=255)
                                                                                                                             tGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACa                          <  2:1341183/150‑1 (MQ=255)
                                                                                                                               gcgcCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGTAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCa                                                   <  2:762044/123‑1 (MQ=255)
                                                                                                                               gcgcCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGc                        >  2:5202126/1‑150 (MQ=255)
                                                                                                                                 gccgccTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTc                      >  1:1006310/1‑150 (MQ=255)
                                                                                                                                 gccgccTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTc                      >  2:3020100/1‑150 (MQ=255)
                                                                                                                                 gccgccTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTc                      >  1:5180529/1‑150 (MQ=255)
                                                                                                                                 gccgccTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTc                      >  1:3244292/1‑150 (MQ=255)
                                                                                                                                       tGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCt                >  1:4140663/1‑150 (MQ=255)
                                                                                                                                            tgtTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGa           >  1:988268/1‑150 (MQ=255)
                                                                                                                                              ttCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGAcc         >  1:707215/1‑150 (MQ=255)
                                                                                                                                                 agcagGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCAtt      >  2:3530451/1‑150 (MQ=255)
                                                                                                                                                  gcagGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGa                            >  1:4275013/1‑127 (MQ=255)
                                                                                                                                                  gcagGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCAttt     >  2:4299682/1‑150 (MQ=255)
                                                                                                                                                  gcagGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCAttt     >  2:2390071/1‑150 (MQ=255)
                                                                                                                                                     gggCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGt  >  2:3223975/1‑150 (MQ=255)
                                                                                                                                                     gggCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGt  <  1:3031916/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGT  >  NC_002947/5737849‑5738147

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