Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,470,453 C→A 9.3% V351V (GTC→GTA)  ydbD → PF10971 family putative periplasmic methylglyoxal resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,470,4530CA 9.3% 793.0 / 36.5 225V351V (GTC→GTA) ydbDPF10971 family putative periplasmic methylglyoxal resistance protein
Reads supporting (aligned to +/- strand):  ref base C (95/109);  new base A (6/15);  total (101/124)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 225 total aligned reads displayed.

GTATGGCGAGATATATGCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCTGAATCCATTATTGATGGCTATATCCGCAG  >  CP009273/1470308‑1470599
                                                                                                                                                 |                                                                                                                                                  
gTATGGCGAGATATATGCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCgg                                                                                                                                                  <  1:2327674/148‑1 (MQ=255)
 tATGGCGAGATATATGCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGt                                                                                                                                                 <  2:1474044/148‑1 (MQ=255)
 tATGGCGAGATATATGCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGt                                                                                                                                                 <  2:3357412/148‑1 (MQ=255)
    ggCGAGATATATGCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCg                                                                                                                                              >  2:1339908/1‑148 (MQ=255)
    ggCGAGATATATGCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCg                                                                                                                                              >  2:2047594/1‑148 (MQ=255)
    ggCGAGATATATGCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCg                                                                                                                                              >  2:2475714/1‑148 (MQ=255)
     gCGAGATATATGCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGt                                                                                                                                             >  1:247909/1‑148 (MQ=255)
        agaTATATGCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGa                                                                                                                                          <  1:1031309/148‑1 (MQ=255)
           tatatGCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATc                                                                                                                                       <  2:3524241/148‑1 (MQ=255)
           tatatGCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATc                                                                                                                                       <  1:644982/148‑1 (MQ=255)
           tatatGCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATc                                                                                                                                       <  1:1968457/148‑1 (MQ=255)
           tatatGCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATc                                                                                                                                       <  1:2404182/148‑1 (MQ=255)
            atatGCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCt                                                                                                                                      >  2:1831450/1‑148 (MQ=255)
             tatGCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCtt                                                                                                                                     >  1:1126963/1‑148 (MQ=255)
              atGCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTg                                                                                                                                    <  2:887761/148‑1 (MQ=255)
              atGCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTg                                                                                                                                    >  1:2827680/1‑148 (MQ=255)
                gCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTgcg                                                                                                                                  <  2:2280204/148‑1 (MQ=255)
                gCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTgcg                                                                                                                                  <  1:2466693/148‑1 (MQ=255)
                gCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTAGGTCCGTTGAATCTTgcg                                                                                                                                  <  2:2428781/148‑1 (MQ=255)
                   tGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTa                                                                                                                               >  1:2751708/1‑148 (MQ=255)
                   tGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTa                                                                                                                               >  1:746569/1‑148 (MQ=255)
                   tGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTa                                                                                                                               <  1:725667/148‑1 (MQ=255)
                    gACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTAc                                                                                                                              >  2:2809376/1‑148 (MQ=255)
                       aaTGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTAGGTCCGTTGAATCTTGCGCTACGTg                                                                                                                           <  2:2217565/148‑1 (MQ=255)
                        aTGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGa                                                                                                                          >  2:708244/1‑148 (MQ=255)
                        aTGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGa                                                                                                                          >  1:131929/1‑148 (MQ=255)
                         tGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGaa                                                                                                                         >  1:862712/1‑148 (MQ=255)
                         tGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGaa                                                                                                                         >  1:24378/1‑148 (MQ=255)
                         tGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGaa                                                                                                                         >  1:1678997/1‑148 (MQ=255)
                         tGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGaa                                                                                                                         >  1:3214561/1‑148 (MQ=255)
                          gagaTTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTAGGTCCGTTGACTCTTGCGCTACGTGaaa                                                                                                                        >  2:977607/1‑148 (MQ=255)
                           agaTTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACAGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGaaaa                                                                                                                       >  2:3504429/1‑148 (MQ=255)
                             aTTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAAtc                                                                                                                     <  2:363295/148‑1 (MQ=255)
                             aTTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTAGGTCCGTTGAATCTTGCGCTACGTGAAAAtc                                                                                                                     >  2:3073019/1‑148 (MQ=255)
                                  ggAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCct                                                                                                                <  1:821625/148‑1 (MQ=255)
                                  ggAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAAAAATGACGCGACTGAGTCGTCTACCTGATGGTGACGGTCCGTTGAATCGTGCGCTCCGAGAAAATCTcccg                                                                                                                >  2:1924019/1‑147 (MQ=255)
                                   gAGTCGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCctc                                                                                                               <  2:2714609/148‑1 (MQ=255)
                                     gTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTg                                                                                                             <  2:866527/148‑1 (MQ=255)
                                     gTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTg                                                                                                             <  2:2763402/148‑1 (MQ=255)
                                      ttGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTgg                                                                                                            >  2:3301504/1‑148 (MQ=255)
                                       tGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGa                                                                                                           >  2:2547364/1‑148 (MQ=255)
                                       tGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTAGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGa                                                                                                           <  2:1554223/148‑1 (MQ=255)
                                          ccGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGatat                                                                                                        <  1:2208818/148‑1 (MQ=255)
                                           cGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCACTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATc                                                                                                       >  2:339865/1‑148 (MQ=255)
                                              tCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTc                                                                                                    <  2:1644917/148‑1 (MQ=255)
                                              tCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTc                                                                                                    >  1:2107308/1‑148 (MQ=255)
                                              tCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTAGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTc                                                                                                    <  1:2496041/148‑1 (MQ=255)
                                                aGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCAt                                                                                                  >  1:1774900/1‑148 (MQ=255)
                                                   aaTTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGAc                                                                                               <  1:1899584/148‑1 (MQ=255)
                                                    aTTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGAct                                                                                              >  1:1044727/1‑148 (MQ=255)
                                                     ttGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGActc                                                                                             >  1:1825049/1‑148 (MQ=255)
                                                      tGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGActct                                                                                            >  1:2031928/1‑148 (MQ=255)
                                                      tGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAAACTCCCTCTGGATATAGTCAAGActat                                                                                            >  1:2975948/1‑146 (MQ=255)
                                                       gggTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTg                                                                                           <  1:932706/148‑1 (MQ=255)
                                                              cAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGaaaa                                                                                    <  1:1242628/148‑1 (MQ=255)
                                                                 ccTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAagag                                                                                 <  2:615891/148‑1 (MQ=255)
                                                                 ccTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAagag                                                                                 >  1:3032564/1‑148 (MQ=255)
                                                                 ccTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTAGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAagag                                                                                 <  2:2123168/148‑1 (MQ=255)
                                                                    gttgttCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGAtt                                                                              >  2:1510080/1‑148 (MQ=255)
                                                                     ttgttCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTg                                                                             <  1:3109122/148‑1 (MQ=255)
                                                                     ttgttCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTAGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTg                                                                             <  1:702131/148‑1 (MQ=255)
                                                                        ttCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGa                                                                          <  2:2679105/148‑1 (MQ=255)
                                                                         tCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGat                                                                         <  1:1177292/148‑1 (MQ=255)
                                                                         tCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGat                                                                         >  1:1662976/1‑148 (MQ=255)
                                                                             gggTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGa                                                                     <  1:2623715/148‑1 (MQ=255)
                                                                                      gACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAgg                                                            <  2:130144/148‑1 (MQ=255)
                                                                                      gACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAgg                                                            <  2:130367/148‑1 (MQ=255)
                                                                                       aCCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGa                                                           >  2:3473413/1‑148 (MQ=255)
                                                                                        ccGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAg                                                          <  2:2151112/148‑1 (MQ=255)
                                                                                         cGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTAGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGc                                                         <  1:2790075/148‑1 (MQ=255)
                                                                                             aaCCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTta                                                     <  1:615463/148‑1 (MQ=255)
                                                                                                   tctctGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGa                                               >  2:2379907/1‑148 (MQ=255)
                                                                                                         tCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTg                                         >  1:1737645/1‑148 (MQ=255)
                                                                                                          cGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTgg                                        <  1:1381263/148‑1 (MQ=255)
                                                                                                          cGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTAGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTgg                                        >  1:1575016/1‑148 (MQ=255)
                                                                                                              ataatGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTAt                                    <  1:971548/148‑1 (MQ=255)
                                                                                                              ataatGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTAt                                    >  1:3454351/1‑148 (MQ=255)
                                                                                                                 atGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGc                                 <  2:1472486/148‑1 (MQ=255)
                                                                                                                 atGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGc                                 <  2:3337781/148‑1 (MQ=255)
                                                                                                                 atGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGc                                 <  1:2959517/148‑1 (MQ=255)
                                                                                                                  tGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGcc                                <  2:1079850/148‑1 (MQ=255)
                                                                                                                   gACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCt                               >  2:2022424/1‑148 (MQ=255)
                                                                                                                   gACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCt                               <  2:167057/148‑1 (MQ=255)
                                                                                                                    aCGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCTg                              >  2:28891/1‑148 (MQ=255)
                                                                                                                      gcgAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCTGaa                            <  1:85875/148‑1 (MQ=255)
                                                                                                                           tGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCTGAATCCat                       >  2:3027238/1‑148 (MQ=255)
                                                                                                                           tGAGTCGTCTACCTGATGGTGTAGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCTGAATCCat                       <  1:2055064/148‑1 (MQ=255)
                                                                                                                            gAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCTGAATCCatt                      >  2:315834/1‑148 (MQ=255)
                                                                                                                             aGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCTGAATCCatta                     <  2:3072137/148‑1 (MQ=255)
                                                                                                                             aGTCGTCTACCTGATGGTGTAGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCTGAATCCatta                     >  2:727368/1‑148 (MQ=255)
                                                                                                                                cgtcTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCTGAATCCATTATTg                  >  1:1499344/1‑148 (MQ=255)
                                                                                                                                     aCCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCTGAATCCATTATTGATGGc             <  1:1961893/148‑1 (MQ=255)
                                                                                                                                     aCCTGATGGTGTAGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCTGAATCCATTATTGATGGc             >  1:3375622/1‑148 (MQ=255)
                                                                                                                                      ccTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCTGAATCCATTATTGATGGCt            <  2:1268916/148‑1 (MQ=255)
                                                                                                                                        tGATGGTGTCGGTCCGTTTAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCTGAATCCATTATTGATGGCtat          <  2:2586091/148‑1 (MQ=255)
                                                                                                                                         gATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCTGAATCCATTATTGATGGCtata         >  1:126296/1‑148 (MQ=255)
                                                                                                                                         gATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCTGAATCCATTATTGATGGCtata         >  1:340107/1‑148 (MQ=255)
                                                                                                                                             gtgtCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCTGAATCCATTATTGATGGCTATATCCg     <  2:647786/148‑1 (MQ=255)
                                                                                                                                                tCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCTGAATCCATTATTGATGGCTATATCCGCAg  <  1:464441/148‑1 (MQ=255)
                                                                                                                                                tCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCTGAATCCATTATTGATGGCTATATCCGCAg  <  2:862670/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GTATGGCGAGATATATGCCTGACAATGAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGCTTTATAAAGATGCGTGGTTATTGCCTGAATCCATTATTGATGGCTATATCCGCAG  >  CP009273/1470308‑1470599

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