Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,138,039 C→A 23.9% 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,0390CA23.9% 506.0 / 116.7 213I431I (ATC→ATAyegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (68/94);  new base A (16/35);  total (84/129)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 217 total aligned reads displayed.

CCGCAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAGTGCGTGCTG  >  CP009273/2137892‑2138186
                                                                                                                                                   |                                                                                                                                                   
ccGCAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACAGAACAGGTCAATCAGCAACTGATGGAGCGCATc                                                                                                                                                     >  1:70798/1‑148 (MQ=255)
   cAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCAct                                                                                                                                                  <  2:2317415/148‑1 (MQ=255)
     aTGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCActct                                                                                                                                                >  1:756412/1‑148 (MQ=255)
      tGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTg                                                                                                                                               <  2:625409/148‑1 (MQ=255)
      tGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTg                                                                                                                                               <  2:1590848/148‑1 (MQ=255)
      tGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTg                                                                                                                                               <  2:2411367/148‑1 (MQ=255)
           aTGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTaa                                                                                                                                          <  2:1344129/148‑1 (MQ=255)
                gtCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAg                                                                                                                                     <  1:2133181/148‑1 (MQ=255)
                   tGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAgcgg                                                                                                                                  >  2:729274/1‑148 (MQ=255)
                    gggCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAgcggg                                                                                                                                 >  2:1415317/1‑148 (MQ=255)
                       cGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTNACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAgcgggcgg                                                                                                                              <  1:1804899/148‑1 (MQ=255)
                       cGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAgcgggcgg                                                                                                                              <  2:2528230/148‑1 (MQ=255)
                        gTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGcgggcggg                                                                                                                             <  1:2018628/148‑1 (MQ=255)
                         ttgcttgcCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGa                                                                                                                            >  2:2144280/1‑148 (MQ=255)
                         ttgcttgcCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGa                                                                                                                            >  2:439712/1‑148 (MQ=255)
                         ttgcttgcCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGa                                                                                                                            <  1:1501419/148‑1 (MQ=255)
                          tgcttgcCGTCTCACTGGTGCGCCACACGGATGGTACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGAt                                                                                                                           >  1:1340087/1‑148 (MQ=255)
                          tgcttgcCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGAt                                                                                                                           >  1:654366/1‑148 (MQ=255)
                                ccGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCAt                                                                                                                     <  2:808761/148‑1 (MQ=255)
                                  gTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCt                                                                                                                   <  1:3212645/148‑1 (MQ=255)
                                   tctcACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtg                                                                                                                  >  2:3515075/1‑148 (MQ=255)
                                    ctcACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgg                                                                                                                 <  2:272951/148‑1 (MQ=255)
                                      cACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtggga                                                                                                               <  1:2689504/148‑1 (MQ=255)
                                       aCTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggag                                                                                                              <  1:1082408/148‑1 (MQ=255)
                                       aCTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggag                                                                                                              >  1:1873218/1‑148 (MQ=255)
                                         tGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggagtg                                                                                                            <  2:1712061/148‑1 (MQ=255)
                                           gTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggagtggg                                                                                                          <  1:479159/148‑1 (MQ=255)
                                           gTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggagtggg                                                                                                          <  1:2268480/148‑1 (MQ=255)
                                             gcgcCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggagtgggag                                                                                                        <  1:287248/148‑1 (MQ=255)
                                              cgcCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTgggagtgggagt                                                                                                       <  1:1869166/148‑1 (MQ=255)
                                                 cacaCGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGa                                                                                                    <  1:2474877/148‑1 (MQ=255)
                                                 cacaCGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGa                                                                                                    <  1:805875/148‑1 (MQ=255)
                                                        aTGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGaa                                                                                             >  2:697834/1‑148 (MQ=255)
                                                        aTGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGaa                                                                                             <  1:2532330/148‑1 (MQ=255)
                                                         tGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAat                                                                                            >  1:2074733/1‑148 (MQ=255)
                                                         tGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAat                                                                                            >  1:281644/1‑148 (MQ=255)
                                                          ggCACGCCGCTCTAGTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAata                                                                                           <  2:3165816/148‑1 (MQ=255)
                                                            cACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATAtt                                                                                         >  2:1169785/1‑148 (MQ=255)
                                                             aCGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATAttt                                                                                        <  2:788968/148‑1 (MQ=255)
                                                              cgccgcTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATAtttt                                                                                       <  1:1158575/148‑1 (MQ=255)
                                                               gccgcTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATAttttt                                                                                      <  2:2658996/148‑1 (MQ=255)
                                                               gccgcTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATAttttt                                                                                      <  2:879020/148‑1 (MQ=255)
                                                                 cgcTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAg                                                                                    <  2:3051824/148‑1 (MQ=255)
                                                                  gcTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGc                                                                                   <  1:1225584/148‑1 (MQ=255)
                                                                    tctATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTg                                                                                 <  1:3209839/148‑1 (MQ=255)
                                                                         tttATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATa                                                                            <  2:3441172/148‑1 (MQ=255)
                                                                         tttATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATa                                                                            >  1:2352628/1‑148 (MQ=255)
                                                                         tttATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATa                                                                            <  2:155951/148‑1 (MQ=255)
                                                                           tATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAg                                                                          <  2:657296/148‑1 (MQ=255)
                                                                            aTCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGc                                                                         <  1:1628783/148‑1 (MQ=255)
                                                                             tCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCg                                                                        >  2:1119267/1‑148 (MQ=255)
                                                                              cGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCgg                                                                       <  2:2604862/148‑1 (MQ=255)
                                                                              cGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCgg                                                                       <  2:375857/148‑1 (MQ=255)
                                                                                 tCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATg                                                                    >  1:1410775/1‑148 (MQ=255)
                                                                                   aGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGtt                                                                  <  1:961363/148‑1 (MQ=255)
                                                                                   aGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGtt                                                                  <  1:2682227/148‑1 (MQ=255)
                                                                                     aTTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCg                                                                >  2:775346/1‑148 (MQ=255)
                                                                                      ttGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGa                                                               >  1:1562065/1‑148 (MQ=255)
                                                                                         aaGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCt                                                            <  1:220076/148‑1 (MQ=255)
                                                                                         aaGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCt                                                            <  1:2545075/148‑1 (MQ=255)
                                                                                         aaGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCt                                                            <  2:1694490/148‑1 (MQ=255)
                                                                                          aGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTg                                                           >  2:2776789/1‑148 (MQ=255)
                                                                                           gACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGt                                                          <  2:2332469/148‑1 (MQ=255)
                                                                                                tAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGaa                                                     <  2:1058983/148‑1 (MQ=255)
                                                                                                tAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGaa                                                     <  2:561158/148‑1 (MQ=255)
                                                                                                  aCGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAAt                                                   <  1:620472/148‑1 (MQ=255)
                                                                                                  aCGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAAt                                                   >  2:3484194/1‑148 (MQ=255)
                                                                                                  aCGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAAt                                                   >  2:3481938/1‑148 (MQ=255)
                                                                                                        tAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATtcctcc                                             <  1:3427191/148‑1 (MQ=255)
                                                                                                        tAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATtcctcc                                             >  2:2236336/1‑148 (MQ=255)
                                                                                                         aaaaCGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCca                                            <  1:669995/148‑1 (MQ=255)
                                                                                                         aaaaCGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCca                                            >  1:952852/1‑148 (MQ=255)
                                                                                                          aaaCGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCcac                                           <  2:2042466/148‑1 (MQ=255)
                                                                                                          aaaCGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCcac                                           >  1:3023470/1‑148 (MQ=255)
                                                                                                           aaCGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCcaca                                          <  1:1906611/148‑1 (MQ=255)
                                                                                                            aCGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACat                                         >  2:2134358/1‑148 (MQ=255)
                                                                                                              gCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACatat                                       >  2:1183877/1‑148 (MQ=255)
                                                                                                                aCCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCa                                     >  1:2348533/1‑148 (MQ=255)
                                                                                                                  cGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCaaa                                   <  2:1351219/148‑1 (MQ=255)
                                                                                                                   gAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAAc                                  <  1:3138849/148‑1 (MQ=255)
                                                                                                                     aCAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCg                                >  1:2475958/1‑148 (MQ=255)
                                                                                                                      cAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGa                               >  1:2150056/1‑148 (MQ=255)
                                                                                                                       aGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGaa                              <  1:2316161/148‑1 (MQ=255)
                                                                                                                        ggccAATCAGCAACTGAGGGAGCGCATCACTCTGGCTAACGCAGCGGGCGGGATTGTCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAAc                             <  2:1854004/145‑1 (MQ=255)
                                                                                                                         gTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGACGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACt                            >  2:883932/1‑148 (MQ=255)
                                                                                                                         gTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACt                            <  2:3271051/148‑1 (MQ=255)
                                                                                                                          tCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTg                           >  1:1973719/1‑148 (MQ=255)
                                                                                                                                      tGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTa               >  1:3041500/1‑148 (MQ=255)
                                                                                                                                      tGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTa               <  1:1682476/148‑1 (MQ=255)
                                                                                                                                       gATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTAc              >  2:475164/1‑148 (MQ=255)
                                                                                                                                        aTGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACg             >  2:1025575/1‑148 (MQ=255)
                                                                                                                                         tGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGa            >  2:2873868/1‑148 (MQ=255)
                                                                                                                                         tGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGa            >  2:3356080/1‑148 (MQ=255)
                                                                                                                                         tGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGa            >  1:3253807/1‑148 (MQ=255)
                                                                                                                                          ggAGCGCATCACTCTGGCTAACGAAGCGTGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGACGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAg           <  1:1039867/148‑1 (MQ=255)
                                                                                                                                          ggAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATGTGGGAGTGGGAGTTGAAGCCGAATTTTTTGAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAATCCGAACTGGCAGGTGTGGTACGAg           <  2:913577/148‑1 (MQ=255)
                                                                                                                                          ggAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAg           <  2:2624538/148‑1 (MQ=255)
                                                                                                                                          ggAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAg           >  1:3543710/1‑148 (MQ=255)
                                                                                                                                              cgcATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGNGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgtgcg       <  1:1801639/148‑1 (MQ=255)
                                                                                                                                              cgcATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgtgcg       <  1:1935891/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCGCAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAGTGCGTGCTG  >  CP009273/2137892‑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: