Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,138,039 C→A 100% 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,0390CA99.1% 693.2 / ‑3.2 216I431I (ATC→ATA) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (107/106);  total (108/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 216 total aligned reads displayed.

GCAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAGTGCGTG  >  CP009273/2137894‑2138183
                                                                                                                                                 |                                                                                                                                                
gCAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAAc                                                                                                                                                >  2:3296302/1‑148 (MQ=255)
 cAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAAct                                                                                                                                               >  2:2231786/1‑148 (MQ=255)
  aaTGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAActc                                                                                                                                              >  2:1468216/1‑148 (MQ=255)
   aTGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAActct                                                                                                                                             >  1:1118041/1‑148 (MQ=255)
     ggCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTgg                                                                                                                                           <  1:404871/148‑1 (MQ=255)
     ggCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTgg                                                                                                                                           >  1:2325228/1‑148 (MQ=255)
      gCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGc                                                                                                                                          <  1:3467039/148‑1 (MQ=255)
         aTGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTaa                                                                                                                                       <  1:2117303/148‑1 (MQ=255)
            ttgtCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGa                                                                                                                                    >  2:376918/1‑148 (MQ=255)
                       ttgcttgcCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGa                                                                                                                         >  2:2936400/1‑148 (MQ=255)
                               cGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATc                                                                                                                 <  2:3468062/148‑1 (MQ=255)
                               cGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATc                                                                                                                 <  2:78276/148‑1 (MQ=255)
                                gTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCt                                                                                                                <  2:1402947/148‑1 (MQ=255)
                                  ctcACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgg                                                                                                              >  2:469841/1‑148 (MQ=255)
                                   tcACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtggg                                                                                                             >  1:1064902/1‑148 (MQ=255)
                                     aCTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggag                                                                                                           >  1:2625472/1‑148 (MQ=255)
                                     aCTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggag                                                                                                           >  1:3046340/1‑148 (MQ=255)
                                     aCTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggag                                                                                                           >  2:3473503/1‑148 (MQ=255)
                                       tGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggagtg                                                                                                         <  1:3373168/148‑1 (MQ=255)
                                       tGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggagtg                                                                                                         >  2:552631/1‑148 (MQ=255)
                                       tGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggagtg                                                                                                         >  2:1666801/1‑148 (MQ=255)
                                            cgcCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTgggagtgggagt                                                                                                    <  1:168022/148‑1 (MQ=255)
                                              ccACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTg                                                                                                  <  2:262034/148‑1 (MQ=255)
                                              ccACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTg                                                                                                  <  1:2644908/148‑1 (MQ=255)
                                              ccACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTg                                                                                                  <  2:2501566/148‑1 (MQ=255)
                                               cacaCGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGa                                                                                                 <  1:202061/148‑1 (MQ=255)
                                                acacGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGaa                                                                                                <  1:653601/148‑1 (MQ=255)
                                                acacGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGaa                                                                                                <  1:821750/148‑1 (MQ=255)
                                                acacGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGaa                                                                                                >  1:1316359/1‑148 (MQ=255)
                                                acacGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGaa                                                                                                <  1:2415066/148‑1 (MQ=255)
                                                  acGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGc                                                                                              <  2:2941824/148‑1 (MQ=255)
                                                   cGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGcc                                                                                             <  1:3100548/148‑1 (MQ=255)
                                                    ggATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAAAGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCg                                                                                            >  1:1487883/1‑148 (MQ=255)
                                                       tGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAat                                                                                         <  2:2363357/148‑1 (MQ=255)
                                                        ggCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAata                                                                                        <  1:3356303/148‑1 (MQ=255)
                                                         gCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAatat                                                                                       <  2:1755444/148‑1 (MQ=255)
                                                            cgccgcTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATAtttt                                                                                    >  2:818405/1‑148 (MQ=255)
                                                             gccgcTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATAttttt                                                                                   <  1:1552502/148‑1 (MQ=255)
                                                             gccgcTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATAttttt                                                                                   <  2:876745/148‑1 (MQ=255)
                                                               cgcTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAg                                                                                 <  2:3149618/148‑1 (MQ=255)
                                                                 ctctATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCt                                                                               <  1:522449/148‑1 (MQ=255)
                                                                  tctATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTg                                                                              >  1:2909300/1‑148 (MQ=255)
                                                                     aTTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGa                                                                           >  1:704486/1‑148 (MQ=255)
                                                                       tttATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATa                                                                         >  2:1312449/1‑148 (MQ=255)
                                                                          aTCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGc                                                                      >  2:1510974/1‑148 (MQ=255)
                                                                            cGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCgg                                                                    >  2:1668470/1‑148 (MQ=255)
                                                                               tCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATg                                                                 >  2:1698372/1‑148 (MQ=255)
                                                                                cAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGCTAAGCGGATGt                                                                <  1:745067/148‑1 (MQ=255)
                                                                                cAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGt                                                                <  2:1348403/148‑1 (MQ=255)
                                                                                   aTTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCg                                                             <  2:3141246/148‑1 (MQ=255)
                                                                                    ttGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGa                                                            >  1:1439480/1‑148 (MQ=255)
                                                                                    ttGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGAGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGTATAAGCGGATGTTCGa                                                            >  2:3414117/1‑148 (MQ=255)
                                                                                     tGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAg                                                           >  1:781278/1‑148 (MQ=255)
                                                                                      gAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGc                                                          >  1:2463879/1‑148 (MQ=255)
                                                                                        aGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTg                                                        >  2:679649/1‑148 (MQ=255)
                                                                                        aGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTg                                                        <  2:2528021/148‑1 (MQ=255)
                                                                                        aGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTg                                                        >  2:687463/1‑148 (MQ=255)
                                                                                        aGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTg                                                        >  2:485967/1‑148 (MQ=255)
                                                                                          aCATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTa                                                      <  2:222037/148‑1 (MQ=255)
                                                                                          aCATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTa                                                      >  1:5779/1‑148 (MQ=255)
                                                                                            aTTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATg                                                    >  2:2714412/1‑148 (MQ=255)
                                                                                              tAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGaa                                                  >  1:215762/1‑148 (MQ=255)
                                                                                              tAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGaa                                                  >  1:991916/1‑148 (MQ=255)
                                                                                              tAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGaa                                                  >  1:2745879/1‑148 (MQ=255)
                                                                                               aaCGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGaaa                                                 <  2:110085/148‑1 (MQ=255)
                                                                                                aCGAGCTAAAACGCACCGACCAGGTCCCTCAGCAACTGAGGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAAt                                                <  1:2672746/148‑1 (MQ=255)
                                                                                                   aGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATtcc                                             >  2:3300773/1‑148 (MQ=255)
                                                                                                     cTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATtcctc                                           <  1:978762/148‑1 (MQ=255)
                                                                                                     cTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATtcctc                                           <  2:308881/148‑1 (MQ=255)
                                                                                                     cTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATtcctc                                           <  2:308796/148‑1 (MQ=255)
                                                                                                        aaaCGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCcac                                        >  2:2283999/1‑148 (MQ=255)
                                                                                                         aaCGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCcaca                                       >  2:265029/1‑148 (MQ=255)
                                                                                                         aaCGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCcaca                                       >  2:1575390/1‑148 (MQ=255)
                                                                                                          aCGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACat                                      >  1:1574142/1‑148 (MQ=255)
                                                                                                            gCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACatat                                    >  2:3463991/1‑148 (MQ=255)
                                                                                                             cACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATc                                   >  1:1002672/1‑148 (MQ=255)
                                                                                                               ccGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCaa                                 >  1:936593/1‑148 (MQ=255)
                                                                                                                    cAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGa                            <  2:515568/148‑1 (MQ=255)
                                                                                                                    cAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGa                            >  1:468804/1‑148 (MQ=255)
                                                                                                                     aGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGaa                           >  1:2871744/1‑148 (MQ=255)
                                                                                                                      ggTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAAc                          >  2:2164301/1‑148 (MQ=255)
                                                                                                                      ggTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAAACGAAc                          >  2:2536752/1‑148 (MQ=255)
                                                                                                                       gTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACt                         <  1:429129/148‑1 (MQ=255)
                                                                                                                         cAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTgg                       >  2:979862/1‑148 (MQ=255)
                                                                                                                            tCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAg                    >  2:2531809/1‑148 (MQ=255)
                                                                                                                              aGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGgt                  <  2:1454412/148‑1 (MQ=255)
                                                                                                                               gCATCTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGgtg                 <  1:3326470/148‑1 (MQ=255)
                                                                                                                                cAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGgtgt                >  1:917339/1‑148 (MQ=255)
                                                                                                                                cAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGgtgt                >  1:2649467/1‑148 (MQ=255)
                                                                                                                                 aaCTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGtgtg               >  2:3054120/1‑148 (MQ=255)
                                                                                                                                   cTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGt             >  1:1860793/1‑148 (MQ=255)
                                                                                                                                   cTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGt             >  2:1842695/1‑148 (MQ=255)
                                                                                                                                      aTGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACg          >  1:1127990/1‑148 (MQ=255)
                                                                                                                                       tGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGa         <  2:418281/148‑1 (MQ=255)
                                                                                                                                       tGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGa         <  1:245571/148‑1 (MQ=255)
                                                                                                                                         gAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgt       <  1:1650048/148‑1 (MQ=255)
                                                                                                                                            cgcATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGTGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgtgcg    <  1:1916623/148‑1 (MQ=255)
                                                                                                                                            cgcATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgtgcg    <  1:3299977/148‑1 (MQ=255)
                                                                                                                                              cATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgtgcgtg  <  1:3473503/148‑1 (MQ=255)
                                                                                                                                              cATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgtgcgtg  <  1:3035906/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAGTGCGTG  >  CP009273/2137894‑2138183

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