Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,501,358 G→T 5.4% R6S (CGC→AGC)  fryC ← putative enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,501,3580GT 5.4% 1211.1 / 25.4 313R6S (CGC→AGC) fryCputative enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base G (120/176);  new base T (7/10);  total (127/186)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 317 total aligned reads displayed.

ACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCG  >  CP009273/2501210‑2501506
                                                                                                                                                    |                                                                                                                                                    
acacCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCg                                                                                                                                                      >  2:2999750/1‑149 (MQ=255)
 cacCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGt                                                                                                                                                     <  2:3513410/149‑1 (MQ=255)
   ccGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCtttt                                                                                                                                                   <  1:2656598/149‑1 (MQ=255)
   ccGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGttt                                                                                                                                                   <  2:3815958/149‑1 (MQ=255)
     gCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGttttt                                                                                                                                                 >  1:2478391/1‑149 (MQ=255)
      cccATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCTtttttt                                                                                                                                                >  2:4779712/1‑149 (MQ=255)
       ccATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTa                                                                                                                                               >  2:1963215/1‑149 (MQ=255)
        cATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTaa                                                                                                                                              <  1:79919/149‑1 (MQ=255)
        cATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTaa                                                                                                                                              <  1:2230899/149‑1 (MQ=255)
          tAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATg                                                                                                                                            <  1:1237242/149‑1 (MQ=255)
           aaTTAAGGTCGGCACCATGCGAGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATgg                                                                                                                                           >  2:2123057/1‑149 (MQ=255)
            aTTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGc                                                                                                                                          >  1:253384/1‑149 (MQ=255)
             ttAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGcc                                                                                                                                         <  1:880298/149‑1 (MQ=255)
             ttAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGcc                                                                                                                                         <  1:3546/149‑1 (MQ=255)
              tAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCa                                                                                                                                        >  1:3978467/1‑149 (MQ=255)
               aaGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCAt                                                                                                                                       >  1:2492334/1‑149 (MQ=255)
                  gTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAAc                                                                                                                                    <  1:4212778/149‑1 (MQ=255)
                     ggCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAAt                                                                                                                                 <  2:229351/149‑1 (MQ=255)
                     ggCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAAt                                                                                                                                 <  2:3716574/149‑1 (MQ=255)
                        aCCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTAt                                                                                                                              >  1:4572540/1‑149 (MQ=255)
                        aCCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTAt                                                                                                                              <  2:340199/149‑1 (MQ=255)
                           aTGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATccc                                                                                                                           <  2:3963523/149‑1 (MQ=255)
                           aTGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATccc                                                                                                                           >  2:441432/1‑149 (MQ=255)
                             gCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTg                                                                                                                         <  1:4779881/149‑1 (MQ=255)
                               gTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTc                                                                                                                       <  1:4729674/149‑1 (MQ=255)
                               gTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTc                                                                                                                       >  2:3224109/1‑149 (MQ=255)
                                tGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCg                                                                                                                      >  1:2383314/1‑149 (MQ=255)
                                tGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCg                                                                                                                      >  2:3664823/1‑149 (MQ=255)
                                tGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGGTGCACTGCGATGGTTAAAGGCCATAACAATAATCCCTGTCg                                                                                                                      >  2:4188683/1‑149 (MQ=255)
                                   aaTCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAat                                                                                                                   <  1:1178424/149‑1 (MQ=255)
                                     tCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAatat                                                                                                                 >  2:3095177/1‑149 (MQ=255)
                                       ccTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTa                                                                                                               <  2:1177242/149‑1 (MQ=255)
                                       ccTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTa                                                                                                               <  2:2520670/149‑1 (MQ=255)
                                       ccTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTa                                                                                                               <  2:456471/149‑1 (MQ=255)
                                         tGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACt                                                                                                             <  1:108047/149‑1 (MQ=255)
                                            cATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACtgtt                                                                                                          <  1:1323641/149‑1 (MQ=255)
                                            cATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACtgtt                                                                                                          <  2:2057118/149‑1 (MQ=255)
                                              tCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACtgttgt                                                                                                        <  1:716950/149‑1 (MQ=255)
                                              tCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACtgttgt                                                                                                        >  1:2194158/1‑149 (MQ=255)
                                              tCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACtgttgt                                                                                                        <  2:1313838/149‑1 (MQ=255)
                                               cacaTGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTgttgtt                                                                                                       >  1:4135920/1‑149 (MQ=255)
                                                    gctgcGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCagag                                                                                                  <  1:2001749/149‑1 (MQ=255)
                                                       gcGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCa                                                                                               <  2:1687778/149‑1 (MQ=255)
                                                         gggAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAat                                                                                             <  1:4268383/149‑1 (MQ=255)
                                                          ggAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAatc                                                                                            >  1:358630/1‑149 (MQ=255)
                                                             gTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAatcatc                                                                                         <  1:3756384/149‑1 (MQ=255)
                                                             gTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAatcatc                                                                                         >  1:3438029/1‑149 (MQ=255)
                                                              ttCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCAtct                                                                                        >  1:4486548/1‑149 (MQ=255)
                                                                    ccAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGa                                                                                  >  1:3220378/1‑149 (MQ=255)
                                                                     cAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGAt                                                                                 <  2:4359212/149‑1 (MQ=255)
                                                                           cTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCGGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCttt                                                                           <  2:1527594/149‑1 (MQ=255)
                                                                             ggTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTTATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTtga                                                                         <  1:1820850/149‑1 (MQ=255)
                                                                             ggTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTtga                                                                         <  1:626552/149‑1 (MQ=255)
                                                                                   agagGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGc                                                                   >  2:1039509/1‑149 (MQ=255)
                                                                                   agagGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGc                                                                   <  1:3410790/149‑1 (MQ=255)
                                                                                   agagGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGc                                                                   <  1:3015605/149‑1 (MQ=255)
                                                                                                ttCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGttttt                                                      <  2:3790598/149‑1 (MQ=255)
                                                                                                 tCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATTATGCCCGCAGCGTTTTTa                                                     >  1:384910/1‑149 (MQ=255)
                                                                                                  cTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTaa                                                    <  2:198915/149‑1 (MQ=255)
                                                                                                  cTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTaa                                                    >  1:4758907/1‑149 (MQ=255)
                                                                                                   ttAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAAt                                                   <  2:3148714/149‑1 (MQ=255)
                                                                                                   ttAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAAt                                                   >  1:3666090/1‑149 (MQ=255)
                                                                                                    tAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAAtt                                                  <  1:2876277/149‑1 (MQ=255)
                                                                                                       ctgccgctgccgCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCg                                               <  1:4034808/149‑1 (MQ=255)
                                                                                                           cgctgccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCt                                           <  1:2746761/149‑1 (MQ=255)
                                                                                                           cgctgccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCt                                           <  1:1507182/149‑1 (MQ=255)
                                                                                                           cgctgccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCAGAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCt                                           <  1:2085969/149‑1 (MQ=255)
                                                                                                             ctgccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGc                                         <  1:934717/149‑1 (MQ=255)
                                                                                                             ctgccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGc                                         <  1:3324394/149‑1 (MQ=255)
                                                                                                              tgccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCa                                        >  2:1513908/1‑149 (MQ=255)
                                                                                                              tgccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCa                                        <  1:3963314/149‑1 (MQ=255)
                                                                                                              tgccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCa                                        >  2:2380637/1‑149 (MQ=255)
                                                                                                                   cACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTg                                   <  1:2636905/149‑1 (MQ=255)
                                                                                                                   cACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTg                                   >  1:2817386/1‑149 (MQ=255)
                                                                                                                   cACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTg                                   >  2:1075740/1‑149 (MQ=255)
                                                                                                                      cGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGAtt                                <  2:4052995/149‑1 (MQ=255)
                                                                                                                      cGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGAtt                                >  2:2231006/1‑149 (MQ=255)
                                                                                                                      cGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGAtt                                <  1:2058719/149‑1 (MQ=255)
                                                                                                                       ggATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGAttt                               >  2:1697002/1‑149 (MQ=255)
                                                                                                                        gATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTc                              >  2:2725576/1‑149 (MQ=255)
                                                                                                                        gATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTc                              <  2:1598870/149‑1 (MQ=255)
                                                                                                                        gATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTc                              >  2:664221/1‑149 (MQ=255)
                                                                                                                         aTGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCa                             <  1:2999750/149‑1 (MQ=255)
                                                                                                                          tGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCAt                            <  1:3210972/149‑1 (MQ=255)
                                                                                                                            cgcCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATaa                          <  2:2606732/149‑1 (MQ=255)
                                                                                                                             gcCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATaaa                         >  1:2619504/1‑149 (MQ=255)
                                                                                                                               cAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATTATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACg                       <  1:423711/149‑1 (MQ=255)
                                                                                                                                aGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGt                      <  1:524918/149‑1 (MQ=255)
                                                                                                                                aGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGt                      <  1:1450117/149‑1 (MQ=255)
                                                                                                                                 ggCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTc                     >  2:4458319/1‑149 (MQ=255)
                                                                                                                                 ggCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTc                     >  1:1249412/1‑149 (MQ=255)
                                                                                                                                        cGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTTCTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGAt              <  1:3183929/149‑1 (MQ=255)
                                                                                                                                         gggtGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGAtt             <  2:2596490/146‑1 (MQ=255)
                                                                                                                                          gTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTc            >  2:3048062/1‑149 (MQ=255)
                                                                                                                                           ttGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCg           >  1:3513039/1‑149 (MQ=255)
                                                                                                                                            tGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATCTGACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGa          <  2:765689/149‑1 (MQ=255)
                                                                                                                                              cACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGaaa        <  1:1023492/149‑1 (MQ=255)
                                                                                                                                               aCTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAAc       >  1:4331747/1‑149 (MQ=255)
                                                                                                                                                cTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAAcg      >  2:4498305/1‑149 (MQ=255)
                                                                                                                                                  gCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAAcgtt    <  1:1691580/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCG  >  CP009273/2501210‑2501506

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: