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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,241,3970TA13.7% 597.6 / 69.2 208H479L (CAC→CTC) treAperiplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base T (101/78);  new base A (15/14);  total (116/92)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 213 total aligned reads displayed.

ATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCA  >  CP009273/1241250‑1241544
                                                                                                                                                   |                                                                                                                                                   
aTTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGt                                                                                                                                                     >  2:672884/1‑148 (MQ=255)
aTTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGt                                                                                                                                                     >  2:672752/1‑148 (MQ=255)
 ttttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTg                                                                                                                                                    >  1:1040315/1‑148 (MQ=255)
 ttttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTg                                                                                                                                                    <  1:2324285/148‑1 (MQ=255)
  tttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGc                                                                                                                                                   <  2:10250/148‑1 (MQ=255)
  tttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGc                                                                                                                                                   <  2:1637787/148‑1 (MQ=255)
  tttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGc                                                                                                                                                   >  1:1613201/1‑148 (MQ=255)
           cACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATg                                                                                                                                          >  2:2340102/1‑148 (MQ=255)
           cACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATg                                                                                                                                          >  2:614121/1‑148 (MQ=255)
           cACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATg                                                                                                                                          <  2:3236356/148‑1 (MQ=255)
            aCGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGt                                                                                                                                         >  2:1506404/1‑148 (MQ=255)
             cGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTc                                                                                                                                        <  1:693332/148‑1 (MQ=255)
               ccATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCa                                                                                                                                      >  1:55756/1‑148 (MQ=255)
               ccATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCa                                                                                                                                      <  2:2269956/148‑1 (MQ=255)
                cATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCAt                                                                                                                                     <  2:36635/148‑1 (MQ=255)
                cATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCAt                                                                                                                                     <  1:3091270/148‑1 (MQ=255)
                cATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCAt                                                                                                                                     >  1:3009169/1‑148 (MQ=255)
                   tGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGc                                                                                                                                  >  1:2396811/1‑148 (MQ=255)
                      tCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCAc                                                                                                                               <  2:35347/148‑1 (MQ=255)
                       ccAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCAcc                                                                                                                              >  1:1280689/1‑148 (MQ=255)
                       ccAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCAcc                                                                                                                              <  2:1794431/148‑1 (MQ=255)
                         aGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACctc                                                                                                                            >  1:1939014/1‑148 (MQ=255)
                          gCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACctct                                                                                                                           <  2:2413272/148‑1 (MQ=255)
                          gCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACctct                                                                                                                           >  2:191748/1‑148 (MQ=255)
                          gCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACctct                                                                                                                           <  1:3480997/148‑1 (MQ=255)
                           ccAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCtt                                                                                                                          <  2:1390794/148‑1 (MQ=255)
                            cAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCttt                                                                                                                         <  1:1699400/148‑1 (MQ=255)
                             aaaGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCtttt                                                                                                                        >  2:1898438/1‑148 (MQ=255)
                              aaGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCttttt                                                                                                                       >  1:2849512/1‑148 (MQ=255)
                               aGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTg                                                                                                                      <  1:1160377/148‑1 (MQ=255)
                                gCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGc                                                                                                                     >  2:2468727/1‑148 (MQ=255)
                                 ccATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGcc                                                                                                                    >  2:2633187/1‑148 (MQ=255)
                                 ccATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGcc                                                                                                                    <  2:2994941/148‑1 (MQ=255)
                                   aTCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCg                                                                                                                  <  1:2827037/148‑1 (MQ=255)
                                   aTCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCg                                                                                                                  >  1:74138/1‑148 (MQ=255)
                                    tCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGt                                                                                                                 <  2:3090660/148‑1 (MQ=255)
                                     ccTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTa                                                                                                                <  2:2441737/148‑1 (MQ=255)
                                     ccTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTa                                                                                                                >  2:2082013/1‑148 (MQ=255)
                                             ggATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTa                                                                                                        <  2:1245703/148‑1 (MQ=255)
                                              gATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTaa                                                                                                       <  2:2038428/148‑1 (MQ=255)
                                               atatTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAAt                                                                                                      >  1:3155030/1‑148 (MQ=255)
                                               atatTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAAt                                                                                                      <  1:774192/148‑1 (MQ=255)
                                               atatTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAAt                                                                                                      >  2:460783/1‑148 (MQ=255)
                                                tatTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATc                                                                                                     >  1:2437644/1‑148 (MQ=255)
                                                   tCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCtt                                                                                                  <  2:2066482/148‑1 (MQ=255)
                                                    cGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTc                                                                                                 <  2:2704471/148‑1 (MQ=255)
                                                       cACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGt                                                                                              <  2:1319982/148‑1 (MQ=255)
                                                          cgccgcCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTcgc                                                                                           >  1:910580/1‑148 (MQ=255)
                                                                 ccccGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCATCTGTCGCGACCCAc                                                                                    >  1:855291/1‑148 (MQ=255)
                                                                   ccGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTg                                                                                  >  1:513894/1‑148 (MQ=255)
                                                                   ccGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTg                                                                                  >  1:2906222/1‑148 (MQ=255)
                                                                     gTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTa                                                                                >  2:3424359/1‑148 (MQ=255)
                                                                     gTTCCGGCGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTa                                                                                >  2:2509490/1‑148 (MQ=255)
                                                                      ttCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCGGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTaa                                                                               <  1:2479318/148‑1 (MQ=255)
                                                                      ttCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTaa                                                                               >  2:1561970/1‑148 (MQ=255)
                                                                            tggtgCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGc                                                                         >  2:2861400/1‑148 (MQ=255)
                                                                                    aCATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCAt                                                                 >  1:3133708/1‑148 (MQ=255)
                                                                                     catcatATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCAtt                                                                >  1:1366278/1‑148 (MQ=255)
                                                                                      atcatATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCAttt                                                               >  1:1085147/1‑148 (MQ=255)
                                                                                       tcatATTTTTCCACCAGCTTTTTCTCGCGGTCATAGGTGTGCTGAACATTTGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGGCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTg                                                              <  1:2737951/148‑1 (MQ=255)
                                                                                       tcatATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTg                                                              >  1:2073524/1‑148 (MQ=255)
                                                                                        catATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTgg                                                             >  2:2826992/1‑148 (MQ=255)
                                                                                         atatTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGgc                                                            >  1:65329/1‑148 (MQ=255)
                                                                                         atatTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGgc                                                            >  1:65239/1‑148 (MQ=255)
                                                                                          tatTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGgcg                                                           <  1:1397270/148‑1 (MQ=255)
                                                                                          tatTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGgcg                                                           >  1:2578400/1‑148 (MQ=255)
                                                                                            tttttCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCa                                                         >  1:1118366/1‑148 (MQ=255)
                                                                                                 ccaccaGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCa                                                    <  1:746867/148‑1 (MQ=255)
                                                                                                  caccaGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCAt                                                   >  1:3255874/1‑148 (MQ=255)
                                                                                                   accaGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCAtt                                                  >  1:2878808/1‑148 (MQ=255)
                                                                                                    ccaGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCAttg                                                 <  1:2907420/148‑1 (MQ=255)
                                                                                                     caGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCAttgt                                                <  2:1452853/148‑1 (MQ=255)
                                                                                                     caGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCAttgt                                                <  2:2736897/148‑1 (MQ=255)
                                                                                                      aGCTTTTTCTCCCGGTCATATTTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCTCCAACTCTTTTTGCACGTAGTTTTGTAATCAATCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCAttgtt                                               >  2:1574420/1‑148 (MQ=255)
                                                                                                       gCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCAttgttg                                              >  2:2180187/1‑148 (MQ=255)
                                                                                                        cTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGc                                             <  1:3398773/148‑1 (MQ=255)
                                                                                                             tctcCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACt                                        >  1:2246444/1‑148 (MQ=255)
                                                                                                                cccGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACtttt                                     >  2:2965436/1‑148 (MQ=255)
                                                                                                                  cGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCa                                   >  1:3456603/1‑148 (MQ=255)
                                                                                                                  cGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCa                                   <  1:445678/148‑1 (MQ=255)
                                                                                                                  cGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCa                                   >  1:1697682/1‑148 (MQ=255)
                                                                                                                  cGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCa                                   >  2:180966/1‑148 (MQ=255)
                                                                                                                  cGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCa                                   <  2:3223047/148‑1 (MQ=255)
                                                                                                                      cATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGa                               >  2:527843/1‑148 (MQ=255)
                                                                                                                        tAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACg                             <  1:3202452/148‑1 (MQ=255)
                                                                                                                        tAGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACg                             >  2:938155/1‑148 (MQ=255)
                                                                                                                         aGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACgt                            >  1:942835/1‑148 (MQ=255)
                                                                                                                         aGGTGTGCTGAACATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACgt                            <  2:1745352/148‑1 (MQ=255)
                                                                                                                           gtgtGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACgtgg                          <  2:379585/148‑1 (MQ=255)
                                                                                                                           gtgtGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTATAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACgtgg                          >  2:3466463/1‑148 (MQ=255)
                                                                                                                                   aaCATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAgg                  <  1:2468727/148‑1 (MQ=255)
                                                                                                                                   aaCATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCCCTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAgg                  <  1:2613881/148‑1 (MQ=255)
                                                                                                                                   aaCATTGGTCAGGAAGAGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAgg                  <  1:2614027/148‑1 (MQ=255)
                                                                                                                                        tGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGgccgcc             >  1:713566/1‑148 (MQ=255)
                                                                                                                                           tCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCggg          <  1:2868746/148‑1 (MQ=255)
                                                                                                                                           tCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCATCTGTCGCGACCCACTGTAACGGGGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCggg          <  2:855291/148‑1 (MQ=255)
                                                                                                                                            cAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGt         >  2:1046288/1‑148 (MQ=255)
                                                                                                                                              ggAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTg       <  2:3286725/148‑1 (MQ=255)
                                                                                                                                                  gTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTgcagc   >  1:1897675/1‑148 (MQ=255)
                                                                                                                                                   tGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTgcagca  >  1:1996976/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCA  >  CP009273/1241250‑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: