Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,501,358 G→T 15.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,3580GT15.4% 793.1 / 86.5 247R6S (CGC→AGC) fryCputative enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base G (81/128);  new base T (16/22);  total (97/150)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 253 total aligned reads displayed.

CACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCG  >  CP009273/2501209‑2501506
                                                                                                                                                     |                                                                                                                                                    
cacaCCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCg                                                                                                                                                      >  2:1671001/1‑150 (MQ=255)
cacaCCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCg                                                                                                                                                      <  2:57069/150‑1 (MQ=255)
   acCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGttt                                                                                                                                                   <  2:949628/150‑1 (MQ=255)
     cGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGttttt                                                                                                                                                 >  2:615409/1‑150 (MQ=255)
     cGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGttttt                                                                                                                                                 >  1:2036606/1‑150 (MQ=255)
      gCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGtttttt                                                                                                                                                <  1:3685949/150‑1 (MQ=255)
       cccATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGATGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTa                                                                                                                                               <  2:1676911/150‑1 (MQ=255)
        ccATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTaa                                                                                                                                              >  1:4326251/1‑150 (MQ=255)
        ccATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTaa                                                                                                                                              <  1:38323/150‑1 (MQ=255)
         cATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAAt                                                                                                                                             >  2:1365848/1‑150 (MQ=255)
            aaTTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGc                                                                                                                                          <  1:3012622/150‑1 (MQ=255)
            aaTTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGc                                                                                                                                          >  2:455406/1‑150 (MQ=255)
            aaTTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGc                                                                                                                                          >  2:2557739/1‑150 (MQ=255)
              ttAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGGTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCa                                                                                                                                        >  1:773178/1‑150 (MQ=255)
                aaGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATa                                                                                                                                      <  2:4058501/150‑1 (MQ=255)
                 aGGTCGGCACCATGCGTGAAATCCCTGACATCACCTGCTGCGGGAGGTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATaa                                                                                                                                     <  2:3042304/150‑1 (MQ=255)
                     cGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAAt                                                                                                                                 <  1:4494917/150‑1 (MQ=255)
                     cGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAAt                                                                                                                                 >  1:1302573/1‑150 (MQ=255)
                       gCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTa                                                                                                                               >  2:406379/1‑150 (MQ=255)
                        cACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTAt                                                                                                                              >  1:3313834/1‑150 (MQ=255)
                        cACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTAt                                                                                                                              <  1:1058177/150‑1 (MQ=255)
                        cACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTAt                                                                                                                              <  1:1493207/150‑1 (MQ=255)
                        cACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTAt                                                                                                                              >  2:2936498/1‑150 (MQ=255)
                         aCCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATc                                                                                                                             <  2:17538/150‑1 (MQ=255)
                          ccATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATcc                                                                                                                            <  2:3169342/150‑1 (MQ=255)
                              gCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCTGt                                                                                                                        <  2:2861070/150‑1 (MQ=255)
                              gCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCTGt                                                                                                                        <  2:2861080/150‑1 (MQ=255)
                              gCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGt                                                                                                                        <  2:2266362/150‑1 (MQ=255)
                               cGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTc                                                                                                                       <  2:2600486/150‑1 (MQ=255)
                                 tGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGa                                                                                                                     <  1:4162206/150‑1 (MQ=255)
                                  gAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCTGTCGaa                                                                                                                    >  1:2964573/1‑150 (MQ=255)
                                       cccTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTa                                                                                                               <  1:4233354/150‑1 (MQ=255)
                                       cccTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTa                                                                                                               <  2:3505045/150‑1 (MQ=255)
                                          tGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACtg                                                                                                            >  2:1613619/1‑150 (MQ=255)
                                           gACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACtgt                                                                                                           <  2:72708/150‑1 (MQ=255)
                                             cATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACtgttg                                                                                                         <  1:1075236/150‑1 (MQ=255)
                                             cATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACtgttg                                                                                                         <  2:2068212/150‑1 (MQ=255)
                                             cATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACtgttg                                                                                                         <  2:216232/150‑1 (MQ=255)
                                             cATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACtgttg                                                                                                         <  2:2352568/150‑1 (MQ=255)
                                                  caTGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCag                                                                                                    >  2:2500976/1‑150 (MQ=255)
                                                    tgctgcGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCagag                                                                                                  <  1:56095/150‑1 (MQ=255)
                                                     gctgcGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAgg                                                                                                 <  1:3779284/150‑1 (MQ=255)
                                                       tgcGGGAGTTCACCCCATAAGCTGGTTTTAGAGGCCTGAGGAATATTAACAGCCGCTGCCGCACCAGATGCGCCAGGCACAACTGTTGCACTGCGTATTTTAATGGCCAAAACAATTATCCCTGTCGAATATAACTGTTGTTCAGAggaa                                                                                               >  2:1377260/1‑148 (MQ=255)
                                                         cGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGGTTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAat                                                                                             >  1:3055380/1‑150 (MQ=255)
                                                         cGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAat                                                                                             <  1:1782013/150‑1 (MQ=255)
                                                             aGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAatcatc                                                                                         <  2:2457994/150‑1 (MQ=255)
                                                              gTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCAtct                                                                                        <  2:1918993/150‑1 (MQ=255)
                                                                 cACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCtt                                                                                     <  1:4262969/150‑1 (MQ=255)
                                                                     ccAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGAt                                                                                 >  2:3703090/1‑150 (MQ=255)
                                                                         aaGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCt                                                                             >  1:1937045/1‑150 (MQ=255)
                                                                          aGCTGGTTATATAGGCCTGCGGATTCTAAACTGCCGCAGCAGCACCGGAAGCGCCAGGCACAACGGTTGCACTGCGTTTTATAATGGCCAAAACAAATATCCCTCTCGAATATTACTGTAGTACAGATGCACACATCTCTACGATTTCtt                                                                            >  2:3856668/1‑150 (MQ=255)
                                                                           gCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCttt                                                                           <  1:2784643/150‑1 (MQ=255)
                                                                           gCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCttt                                                                           <  1:475445/150‑1 (MQ=255)
                                                                            cTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTtg                                                                          <  2:1155213/150‑1 (MQ=255)
                                                                             tGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATGTCTTtga                                                                         >  1:2458813/1‑150 (MQ=255)
                                                                              ggTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTtgat                                                                        >  2:1977683/1‑150 (MQ=255)
                                                                               gTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTtgatg                                                                       >  2:827751/1‑150 (MQ=255)
                                                                               gTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTtgatg                                                                       >  1:3185699/1‑150 (MQ=255)
                                                                                ttttAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTtgatga                                                                      <  2:65664/150‑1 (MQ=255)
                                                                                 tttAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTgatgat                                                                     <  1:3239536/150‑1 (MQ=255)
                                                                                   tAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGc                                                                   <  1:2509863/150‑1 (MQ=255)
                                                                                   tAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGc                                                                   >  2:4469007/1‑150 (MQ=255)
                                                                                     gagGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGccc                                                                 >  1:3953975/1‑150 (MQ=255)
                                                                                      agGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCg                                                                <  1:270389/150‑1 (MQ=255)
                                                                                             cGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGtt                                                         >  2:3120070/1‑150 (MQ=255)
                                                                                             cGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGtt                                                         >  2:1428148/1‑150 (MQ=255)
                                                                                              ggATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGttt                                                        >  2:2205033/1‑150 (MQ=255)
                                                                                              ggATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGttt                                                        <  1:1481873/150‑1 (MQ=255)
                                                                                               gATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATCACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGtttt                                                       <  2:1735531/150‑1 (MQ=255)
                                                                                                 ttCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTa                                                     <  2:3585067/150‑1 (MQ=255)
                                                                                                  tCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTaa                                                    <  2:2627316/150‑1 (MQ=255)
                                                                                                    ttAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAAtt                                                  >  1:3260420/1‑150 (MQ=255)
                                                                                                     tAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAAGATAACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTg                                                 >  2:3447207/1‑150 (MQ=255)
                                                                                                             gctgccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCGTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGc                                         >  2:3274122/1‑150 (MQ=255)
                                                                                                                 ccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAg                                     <  1:961570/150‑1 (MQ=255)
                                                                                                                 ccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAg                                     <  1:3233248/150‑1 (MQ=255)
                                                                                                                  cgcACCGGATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGt                                    >  1:3323426/1‑150 (MQ=255)
                                                                                                                   gcACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTg                                   <  2:3199909/150‑1 (MQ=255)
                                                                                                                     aCCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGAt                                 <  1:1322043/150‑1 (MQ=255)
                                                                                                                      ccGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGAtt                                >  1:1196657/1‑150 (MQ=255)
                                                                                                                        ggATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTc                              >  2:3786371/1‑150 (MQ=255)
                                                                                                                        ggATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTc                              >  1:1335047/1‑150 (MQ=255)
                                                                                                                         gATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCa                             <  2:3563722/150‑1 (MQ=255)
                                                                                                                         gATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCa                             <  2:4176307/150‑1 (MQ=255)
                                                                                                                         gATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCa                             <  2:3065599/150‑1 (MQ=255)
                                                                                                                          aTGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCAt                            >  2:2922887/1‑150 (MQ=255)
                                                                                                                            gcgcCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATaa                          <  1:4512966/150‑1 (MQ=255)
                                                                                                                            gcgcCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATaa                          <  2:1562591/150‑1 (MQ=255)
                                                                                                                                cAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGt                      >  2:2407789/1‑150 (MQ=255)
                                                                                                                                 aGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTc                     <  1:3434227/150‑1 (MQ=255)
                                                                                                                                 aGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTc                     <  1:1479483/150‑1 (MQ=255)
                                                                                                                                         cGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGAtt             <  1:3370742/150‑1 (MQ=255)
                                                                                                                                          ggTTGCACTGCGTTTTTTAATGGCCATATCAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTc            <  1:4023546/150‑1 (MQ=255)
                                                                                                                                           gTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCg           <  2:4326251/150‑1 (MQ=255)
                                                                                                                                           gTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCg           >  1:1826285/1‑150 (MQ=255)
                                                                                                                                            ttGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGa          <  1:3122135/150‑1 (MQ=255)
                                                                                                                                             tGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGaa         <  1:2207996/150‑1 (MQ=255)
                                                                                                                                              gCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGaaa        <  1:1457094/150‑1 (MQ=255)
                                                                                                                                                aCTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAAcg      <  2:639764/150‑1 (MQ=255)
                                                                                                                                                aCTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAAcg      >  2:1574540/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCG  >  CP009273/2501209‑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: