Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,138,039 C→A 30.4% I431I (ATC→ATA yegE → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,138,0390CA30.4% 372.9 / 145.7 197I431I (ATC→ATAyegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (63/74);  new base A (20/40);  total (83/114)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 198 total aligned reads displayed.

CGCAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAGTGCGTGCTG  >  CP009273/2137893‑2138186
                                                                                                                                                  |                                                                                                                                                   
cGCAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCa                                                                                                                                                    <  1:2304605/148‑1 (MQ=255)
  cAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTTAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAAct                                                                                                                                                  >  1:1299147/1‑148 (MQ=255)
   aaTGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAActc                                                                                                                                                 <  1:134851/148‑1 (MQ=255)
    aTGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCActct                                                                                                                                                >  2:2949438/1‑148 (MQ=255)
      ggCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTgg                                                                                                                                              >  2:1745118/1‑148 (MQ=255)
           tgttgtCTGGGCGTTGCTTGCCGTCTCACTGGTACGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAAc                                                                                                                                         >  2:624602/1‑148 (MQ=255)
            gttgtCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACg                                                                                                                                        >  1:3131986/1‑148 (MQ=255)
              tgtCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGaa                                                                                                                                      >  2:1062262/1‑148 (MQ=255)
                tCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGCTTGAAGACATTAACGAGCTAAAACGCACCGCACAGGTCAATCAGCACCTGATGGAGCGCATCACTCTGGCTAACGAAgc                                                                                                                                    >  2:1988088/1‑148 (MQ=255)
                        gtgcttgcCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGa                                                                                                                            >  2:823928/2‑148 (MQ=255)
                         tgcttgcCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGAt                                                                                                                           >  1:2373475/1‑148 (MQ=255)
                         tgcttgcCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGAt                                                                                                                           <  2:3198527/148‑1 (MQ=255)
                         tgcttgcCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGAt                                                                                                                           >  2:862228/1‑148 (MQ=255)
                         tgcttgcCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGAt                                                                                                                           >  2:1334331/1‑148 (MQ=255)
                             tgcCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGc                                                                                                                       >  2:2817389/1‑148 (MQ=255)
                             tgcCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGc                                                                                                                       <  1:2847766/148‑1 (MQ=255)
                              gcCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCa                                                                                                                      <  1:1988289/148‑1 (MQ=255)
                                  tctcACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtg                                                                                                                  >  2:1104795/1‑148 (MQ=255)
                                  tctcACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtg                                                                                                                  >  1:2927488/1‑148 (MQ=255)
                                  tctcACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtg                                                                                                                  <  1:2019644/148‑1 (MQ=255)
                                   ctcACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgg                                                                                                                 >  1:1883546/1‑148 (MQ=255)
                                     cACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtggga                                                                                                               >  1:1446106/1‑148 (MQ=255)
                                     cACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtggga                                                                                                               <  2:3121057/148‑1 (MQ=255)
                                      aCTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggag                                                                                                              <  1:2159435/148‑1 (MQ=255)
                                        tGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggagtg                                                                                                            <  1:2314437/148‑1 (MQ=255)
                                          gTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggagtggg                                                                                                          <  1:2149581/148‑1 (MQ=255)
                                          gTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggagtggg                                                                                                          <  1:482949/148‑1 (MQ=255)
                                               ccACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTg                                                                                                     <  2:1221620/148‑1 (MQ=255)
                                                cacaCGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGa                                                                                                    <  1:1686083/148‑1 (MQ=255)
                                                  cacGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAg                                                                                                  >  1:494948/1‑148 (MQ=255)
                                                   acGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGc                                                                                                 <  1:1547598/148‑1 (MQ=255)
                                                   acGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGc                                                                                                 <  1:621166/148‑1 (MQ=255)
                                                    cGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGcc                                                                                                >  2:88502/1‑148 (MQ=255)
                                                       aTGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGaa                                                                                             >  1:788033/1‑148 (MQ=255)
                                                       aTGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGaa                                                                                             >  1:788048/1‑148 (MQ=255)
                                                        tGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAat                                                                                            <  1:2989231/148‑1 (MQ=255)
                                                         ggCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAata                                                                                           <  1:674550/148‑1 (MQ=255)
                                                         ggCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAata                                                                                           <  1:2445439/148‑1 (MQ=255)
                                                          gCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAatat                                                                                          >  1:2494903/1‑148 (MQ=255)
                                                          gCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAatat                                                                                          >  1:2038185/1‑148 (MQ=255)
                                                            aCGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATAttt                                                                                        <  2:2271565/148‑1 (MQ=255)
                                                                cgcTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAg                                                                                    <  1:845615/148‑1 (MQ=255)
                                                                    ctATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTgg                                                                                >  1:195369/1‑148 (MQ=255)
                                                                       ttttATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGAt                                                                             >  2:472374/1‑148 (MQ=255)
                                                                        tttATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATa                                                                            >  2:1597880/1‑148 (MQ=255)
                                                                        tttATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATa                                                                            <  2:81665/148‑1 (MQ=255)
                                                                        tttATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATa                                                                            <  1:3106895/148‑1 (MQ=255)
                                                                         ttATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATaa                                                                           <  1:547394/148‑1 (MQ=255)
                                                                         ttATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATaa                                                                           <  1:780162/148‑1 (MQ=255)
                                                                         ttATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATaa                                                                           <  2:1017109/148‑1 (MQ=255)
                                                                         ttATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATaa                                                                           <  1:1552088/148‑1 (MQ=255)
                                                                          tATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAg                                                                          >  2:773959/1‑148 (MQ=255)
                                                                             cGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCgg                                                                       <  2:2540228/148‑1 (MQ=255)
                                                                                tCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATg                                                                    <  2:1458845/148‑1 (MQ=255)
                                                                                tCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATg                                                                    <  2:1194012/148‑1 (MQ=255)
                                                                                 cAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGt                                                                   <  1:3236306/148‑1 (MQ=255)
                                                                                 cAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGt                                                                   >  1:187210/1‑148 (MQ=255)
                                                                                  aGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGtt                                                                  <  2:1830341/148‑1 (MQ=255)
                                                                                    aTTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCg                                                                <  1:1046120/148‑1 (MQ=255)
                                                                                    aTTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCg                                                                >  1:1215660/1‑148 (MQ=255)
                                                                                    aTTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCg                                                                <  1:3074568/148‑1 (MQ=255)
                                                                                      tGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAg                                                              >  2:2583361/1‑148 (MQ=255)
                                                                                      tGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAg                                                              <  1:824800/148‑1 (MQ=255)
                                                                                        aaGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCt                                                            >  2:1141070/1‑148 (MQ=255)
                                                                                        aaGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCt                                                            >  2:1141061/1‑148 (MQ=255)
                                                                                        aaGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCt                                                            >  2:1141038/1‑148 (MQ=255)
                                                                                         aGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTg                                                           >  1:1999082/1‑148 (MQ=255)
                                                                                          gACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGt                                                          >  1:2544718/1‑148 (MQ=255)
                                                                                             aTTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATg                                                       >  1:321225/1‑148 (MQ=255)
                                                                                                   gAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCGGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATtc                                                 <  1:1559925/148‑1 (MQ=255)
                                                                                                      cTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATtcctc                                              <  1:104984/148‑1 (MQ=255)
                                                                                                       tAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATtcctcc                                             >  1:2410674/1‑148 (MQ=255)
                                                                                                       tAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATtcctcc                                             >  1:1807552/1‑148 (MQ=255)
                                                                                                        aaaaCGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTCTGAAATTCCTCca                                            >  2:1937946/1‑148 (MQ=255)
                                                                                                         aaaCGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCcac                                           >  2:2876002/1‑148 (MQ=255)
                                                                                                          aaCGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCcaca                                          <  2:2373475/148‑1 (MQ=255)
                                                                                                           aCGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACat                                         >  2:321639/1‑148 (MQ=255)
                                                                                                           aCGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACat                                         >  2:2416151/1‑148 (MQ=255)
                                                                                                            cGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACata                                        >  1:537828/1‑148 (MQ=255)
                                                                                                                 cGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCaaa                                   <  2:2067538/148‑1 (MQ=255)
                                                                                                                      aGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGaa                              <  2:580873/148‑1 (MQ=255)
                                                                                                                      aGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGaa                              <  2:2248609/148‑1 (MQ=255)
                                                                                                                        gTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACt                            <  2:2722659/148‑1 (MQ=255)
                                                                                                                           aaTCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGc                         <  1:1745118/148‑1 (MQ=255)
                                                                                                                           aaTCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGc                         >  2:3299226/1‑148 (MQ=255)
                                                                                                                            aTCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCa                        <  1:3178664/148‑1 (MQ=255)
                                                                                                                                  aaCTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGtgtg                  <  1:2302639/148‑1 (MQ=255)
                                                                                                                                  aaCTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGtgtg                  <  1:2196847/148‑1 (MQ=255)
                                                                                                                                  aaCTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGtgtg                  <  2:3049759/148‑1 (MQ=255)
                                                                                                                                  aaCTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGtgtg                  <  2:3011430/148‑1 (MQ=255)
                                                                                                                                  aaCTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGtgtg                  <  2:2775693/148‑1 (MQ=255)
                                                                                                                                   aCTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTgg                 >  2:1367150/1‑148 (MQ=255)
                                                                                                                                    cTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGt                <  2:1315730/148‑1 (MQ=255)
                                                                                                                                       aTGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACg             <  1:3217897/148‑1 (MQ=255)
                                                                                                                                       aTGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACg             >  2:387863/1‑148 (MQ=255)
                                                                                                                                         ggAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGTGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAg           <  1:1146523/148‑1 (MQ=255)
                                                                                                                                         ggAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAg           <  1:1423874/148‑1 (MQ=255)
                                                                                                                                          gAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgt          >  2:1914907/1‑148 (MQ=255)
                                                                                                                                               cATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgtgcgtg     >  2:937121/1‑148 (MQ=255)
                                                                                                                                               cATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATAACAAACCGAACTGGCAGGTGTGGTACGAgtgcgtg     <  1:3100045/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGCAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAGTGCGTGCTG  >  CP009273/2137893‑2138186

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: