Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,241,397 T→A 22.5% H479L (CAC→CTC)  treA ← periplasmic trehalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,241,3970TA22.5% 425.8 / 130.9 200H479L (CAC→CTC) treAperiplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base T (66/89);  new base A (23/22);  total (89/111)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 202 total aligned reads displayed.

TTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCA  >  CP009273/1241251‑1241544
                                                                                                                                                  |                                                                                                                                                   
ttttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTg                                                                                                                                                    >  2:520629/1‑148 (MQ=255)
 tgtCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGc                                                                                                                                                   <  2:3103957/146‑1 (MQ=255)
 tgtCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGc                                                                                                                                                   <  2:2529374/146‑1 (MQ=255)
     aGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGc                                                                                                                                               >  1:3210156/1‑148 (MQ=255)
       cGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTa                                                                                                                                             <  2:815023/148‑1 (MQ=255)
        gTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTaa                                                                                                                                            <  1:3393620/148‑1 (MQ=255)
         tCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAAt                                                                                                                                           <  2:1134537/148‑1 (MQ=255)
          cACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATg                                                                                                                                          <  2:1567494/148‑1 (MQ=255)
          cACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATg                                                                                                                                          <  2:1076265/148‑1 (MQ=255)
          cACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATg                                                                                                                                          <  2:1566541/148‑1 (MQ=255)
            cGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTc                                                                                                                                        <  2:1947913/148‑1 (MQ=255)
             gCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTcc                                                                                                                                       <  1:2216671/148‑1 (MQ=255)
             gCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTcc                                                                                                                                       <  2:3727288/148‑1 (MQ=255)
                 ttGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCg                                                                                                                                   >  1:2553737/1‑148 (MQ=255)
                     tCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCAc                                                                                                                               >  2:421914/1‑148 (MQ=255)
                         gCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACctct                                                                                                                           <  1:85883/148‑1 (MQ=255)
                         gCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACctct                                                                                                                           <  1:2016448/148‑1 (MQ=255)
                           cAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCttt                                                                                                                         <  2:1240741/148‑1 (MQ=255)
                           cAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCttt                                                                                                                         <  1:1754974/148‑1 (MQ=255)
                            aaaGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCtttt                                                                                                                        >  1:60762/1‑148 (MQ=255)
                            aaaGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCtttt                                                                                                                        <  1:721512/148‑1 (MQ=255)
                            aaaGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCtttt                                                                                                                        <  1:719703/148‑1 (MQ=255)
                              aGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTg                                                                                                                      >  1:3630802/1‑148 (MQ=255)
                              aGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTg                                                                                                                      >  2:1899728/1‑148 (MQ=255)
                               gCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGc                                                                                                                     >  2:3581125/1‑148 (MQ=255)
                               gCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGc                                                                                                                     >  1:2470771/1‑148 (MQ=255)
                                ccATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGca                                                                                                                    >  2:138224/1‑147 (MQ=255)
                                   tCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGt                                                                                                                 >  2:2571341/1‑148 (MQ=255)
                                   tCCTGTAATGGATATTCGCCACCGCCGCCCCCCGGTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGt                                                                                                                 <  2:2326315/148‑1 (MQ=255)
                                    ccTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTa                                                                                                                >  2:3147508/1‑148 (MQ=255)
                                    ccTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTa                                                                                                                >  2:3317988/1‑148 (MQ=255)
                                    ccTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTa                                                                                                                >  2:2673427/1‑148 (MQ=255)
                                       gTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGtt                                                                                                             <  1:1776990/148‑1 (MQ=255)
                                       gTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGtt                                                                                                             <  2:1764718/148‑1 (MQ=255)
                                           tGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGt                                                                                                         >  1:102403/1‑148 (MQ=255)
                                            ggATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTa                                                                                                        <  2:3766292/148‑1 (MQ=255)
                                            ggATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTa                                                                                                        >  1:720900/1‑148 (MQ=255)
                                                atTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATcc                                                                                                    <  2:3578257/148‑1 (MQ=255)
                                                atTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATcc                                                                                                    <  2:834983/148‑1 (MQ=255)
                                                atTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATcc                                                                                                    >  1:1067972/1‑148 (MQ=255)
                                                atTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATcc                                                                                                    >  1:1068123/1‑148 (MQ=255)
                                                  tCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCtt                                                                                                  >  1:827514/1‑148 (MQ=255)
                                                  tCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCtt                                                                                                  <  1:404894/148‑1 (MQ=255)
                                                   cGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTc                                                                                                 <  2:1739546/148‑1 (MQ=255)
                                                      cACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGt                                                                                              <  2:1648715/148‑1 (MQ=255)
                                                      cACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGt                                                                                              <  1:593431/148‑1 (MQ=255)
                                                           ccgccCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGa                                                                                         <  1:308058/148‑1 (MQ=255)
                                                              ccccccGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGAccc                                                                                      >  1:1218461/1‑148 (MQ=255)
                                                              ccccccGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTGTTGTAATCCTTCTGTCGCGAccc                                                                                      >  1:1427050/1‑148 (MQ=255)
                                                               cccccGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCa                                                                                     <  1:641497/148‑1 (MQ=255)
                                                               cccccGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCa                                                                                     >  2:2651704/1‑148 (MQ=255)
                                                                    gTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTa                                                                                <  1:168883/148‑1 (MQ=255)
                                                                                  gACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCa                                                                  <  1:2978128/148‑1 (MQ=255)
                                                                                    catcatATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCAtt                                                                <  2:2784320/148‑1 (MQ=255)
                                                                                    catcatATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCAtt                                                                <  1:2868766/148‑1 (MQ=255)
                                                                                      tcatATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTg                                                              >  2:1677901/1‑148 (MQ=255)
                                                                                       catATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTgg                                                             >  1:1134811/1‑148 (MQ=255)
                                                                                       catATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTgg                                                             >  1:1064429/1‑148 (MQ=255)
                                                                                         tatTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGgcg                                                           >  2:3057611/1‑148 (MQ=255)
                                                                                          atTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGgcgc                                                          <  2:2913506/148‑1 (MQ=255)
                                                                                            ttttCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCAt                                                        <  1:1397295/148‑1 (MQ=255)
                                                                                            ttttCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCAt                                                        <  2:2970110/148‑1 (MQ=255)
                                                                                            ttttCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCAt                                                        >  2:1119369/1‑148 (MQ=255)
                                                                                            tttgCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCAt                                                        <  2:1796858/148‑1 (MQ=255)
                                                                                                  accaGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCAtt                                                  >  2:3387699/1‑148 (MQ=255)
                                                                                                    caGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCAttgt                                                >  2:2161043/1‑148 (MQ=255)
                                                                                                       cTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGc                                             >  1:3157287/1‑148 (MQ=255)
                                                                                                       cTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGc                                             >  1:1252940/1‑148 (MQ=255)
                                                                                                       cTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGc                                             <  2:2809668/148‑1 (MQ=255)
                                                                                                         ttttCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGccc                                           >  2:2882913/1‑148 (MQ=255)
                                                                                                           ttCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTGTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCAc                                         <  2:723228/148‑1 (MQ=255)
                                                                                                            tctcCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACt                                        <  1:1486802/148‑1 (MQ=255)
                                                                                                             ctcCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACtt                                       <  2:3630802/148‑1 (MQ=255)
                                                                                                              tcCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACttt                                      <  1:1773445/148‑1 (MQ=255)
                                                                                                              tcCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACttt                                      <  1:1346060/148‑1 (MQ=255)
                                                                                                              tcCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACttt                                      <  1:3250219/148‑1 (MQ=255)
                                                                                                                    tCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCg                                <  2:3467587/148‑1 (MQ=255)
                                                                                                                     cATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCCTTTGGCGCATCCCATTGTTGCCCACTTTTCACCGa                               >  2:1973378/1‑148 (MQ=255)
                                                                                                                       tAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACg                             >  2:49978/1‑148 (MQ=255)
                                                                                                                       tAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACg                             >  2:891853/1‑148 (MQ=255)
                                                                                                                        aGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACgt                            >  1:1644680/1‑148 (MQ=255)
                                                                                                                         ggTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACgtg                           <  1:2742124/148‑1 (MQ=255)
                                                                                                                             tgCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGgtgt                       >  2:636231/1‑148 (MQ=255)
                                                                                                                             tgCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCCTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGgtgt                       <  1:2601731/148‑1 (MQ=255)
                                                                                                                              gCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGtt                      <  1:762583/148‑1 (MQ=255)
                                                                                                                                  aaCATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAgg                  >  1:3522944/1‑148 (MQ=255)
                                                                                                                                  aaCATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAgg                  <  1:520629/148‑1 (MQ=255)
                                                                                                                                  aaCATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAgg                  <  2:2386897/148‑1 (MQ=255)
                                                                                                                                    cATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGgcc                >  1:3562421/1‑148 (MQ=255)
                                                                                                                                      ttGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGgccgc              <  2:3691520/148‑1 (MQ=255)
                                                                                                                                          tCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCTTTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCggg          <  2:1446526/148‑1 (MQ=255)
                                                                                                                                          tCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCggg          <  1:984289/148‑1 (MQ=255)
                                                                                                                                          tCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCggg          <  1:2882913/148‑1 (MQ=255)
                                                                                                                                          tCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCGTTTTGCCCGTAGTTTTGGAATCCGTCTGTCGCGACCCACGGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCggg          <  2:448107/148‑1 (MQ=255)
                                                                                                                                              taaGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTgc      >  1:3034210/2‑148 (MQ=255)
                                                                                                                                              gAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTgc      >  1:1236815/1‑148 (MQ=255)
                                                                                                                                              gAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTgc      >  1:1236796/1‑148 (MQ=255)
                                                                                                                                              gAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTgc      >  1:1236366/1‑148 (MQ=255)
                                                                                                                                               aaGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCACATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTgca     >  2:1469643/1‑148 (MQ=255)
                                                                                                                                                  tGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGGGGTGTTCAGGCCGCCGGGTTgcagca  >  2:2303405/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCA  >  CP009273/1241251‑1241544

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