Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,138,039 C→A 26.0% 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,0390CA26.0% 520.5 / 147.1 238I431I (ATC→ATAyegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (69/107);  new base A (26/36);  total (95/143)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 239 total aligned reads displayed.

ACCGCAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAGTGCGTGCT  >  CP009273/2137891‑2138185
                                                                                                                                                    |                                                                                                                                                  
aCCGCAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATc                                                                                                                                                    >  1:3351592/1‑149 (MQ=255)
   gCAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCAct                                                                                                                                                 >  2:1479595/1‑149 (MQ=255)
   gCAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCAct                                                                                                                                                 >  2:1479508/1‑149 (MQ=255)
    cAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCActc                                                                                                                                                <  1:830633/149‑1 (MQ=255)
    cAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCActc                                                                                                                                                <  1:3545459/149‑1 (MQ=255)
    aaaTGGCGATGTTGTCTGGGCTTTGCTTGCCGTCTCACTGGTTCGCCACACTGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCActc                                                                                                                                                <  2:2479835/148‑1 (MQ=255)
        ggCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGc                                                                                                                                            <  2:3876525/149‑1 (MQ=255)
        ggCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGc                                                                                                                                            <  1:860397/149‑1 (MQ=255)
           gATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTaa                                                                                                                                         <  1:1349600/149‑1 (MQ=255)
             tgttgtCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACg                                                                                                                                       >  1:2822243/1‑149 (MQ=255)
                tgtCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTTAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAg                                                                                                                                    <  2:270774/149‑1 (MQ=255)
                 gtCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAgc                                                                                                                                   <  1:928367/149‑1 (MQ=255)
                     gggCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAgcgggc                                                                                                                               >  2:1795322/1‑149 (MQ=255)
                       gCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAgcgggcgg                                                                                                                             >  2:2841230/1‑149 (MQ=255)
                             cttgcCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTgg                                                                                                                       <  2:3435671/149‑1 (MQ=255)
                                 ccGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATc                                                                                                                   <  1:1910673/149‑1 (MQ=255)
                                 ccGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATc                                                                                                                   <  1:1362837/149‑1 (MQ=255)
                                   gTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtg                                                                                                                 >  2:608356/1‑149 (MQ=255)
                                      tcACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtggga                                                                                                              <  1:1595479/149‑1 (MQ=255)
                                      tcACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtggga                                                                                                              >  1:3068411/1‑149 (MQ=255)
                                       cACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggag                                                                                                             <  2:1490154/149‑1 (MQ=255)
                                       cACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggag                                                                                                             <  2:2594736/149‑1 (MQ=255)
                                       cACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggag                                                                                                             <  2:2412329/149‑1 (MQ=255)
                                        aCTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggagt                                                                                                            <  2:1026062/149‑1 (MQ=255)
                                         cTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggagtg                                                                                                           >  1:3861776/1‑149 (MQ=255)
                                            gTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggagtggga                                                                                                        <  2:3389578/149‑1 (MQ=255)
                                            gTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggagtggga                                                                                                        >  2:3585792/1‑149 (MQ=255)
                                             tGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCtgggagtgggag                                                                                                       <  2:2759586/149‑1 (MQ=255)
                                              gcgcCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTgggagtgggagt                                                                                                      <  2:1947540/149‑1 (MQ=255)
                                               cgcCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGtt                                                                                                     <  2:4105388/149‑1 (MQ=255)
                                               cgcCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACACGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGtt                                                                                                     <  1:3288775/149‑1 (MQ=255)
                                                  cacaCGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGaa                                                                                                  <  2:2334817/149‑1 (MQ=255)
                                                   acacGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAg                                                                                                 <  2:1331180/149‑1 (MQ=255)
                                                     acGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGcc                                                                                               <  1:1765810/149‑1 (MQ=255)
                                                      cGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCg                                                                                              <  2:849712/149‑1 (MQ=255)
                                                         aTGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAat                                                                                           <  1:642166/149‑1 (MQ=255)
                                                          tGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAata                                                                                          <  1:2808585/149‑1 (MQ=255)
                                                          tGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAata                                                                                          <  2:1882603/149‑1 (MQ=255)
                                                           ggCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAatat                                                                                         <  1:2354111/149‑1 (MQ=255)
                                                            gCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATAtt                                                                                        >  1:1861245/1‑149 (MQ=255)
                                                              aCGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATAtttt                                                                                      >  2:4008741/1‑149 (MQ=255)
                                                              aCGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATAtttt                                                                                      >  1:3964529/1‑149 (MQ=255)
                                                                gccgcTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGATCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTa                                                                                    <  2:1555049/149‑1 (MQ=255)
                                                                gacgcTCTATTTTATCGCTAAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGTTTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTa                                                                                    <  1:3727213/147‑1 (MQ=255)
                                                                  cgcTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGc                                                                                  >  1:2698521/1‑149 (MQ=255)
                                                                   gcTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCt                                                                                 <  2:65580/149‑1 (MQ=255)
                                                                     tctATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTgg                                                                               >  1:1267554/1‑149 (MQ=255)
                                                                        aTTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGAt                                                                            <  2:2613748/149‑1 (MQ=255)
                                                                          tttATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATaa                                                                          >  2:2901268/1‑149 (MQ=255)
                                                                          tttATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATaa                                                                          <  1:2256969/149‑1 (MQ=255)
                                                                            tATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGc                                                                        >  1:456028/1‑149 (MQ=255)
                                                                            tATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGc                                                                        <  1:3788787/149‑1 (MQ=255)
                                                                            tATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGc                                                                        >  2:485539/1‑149 (MQ=255)
                                                                             aTCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCg                                                                       <  1:1458394/149‑1 (MQ=255)
                                                                                   cAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGtt                                                                 <  2:3956405/149‑1 (MQ=255)
                                                                                     gATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCg                                                               <  1:108473/149‑1 (MQ=255)
                                                                                        tGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGc                                                            >  2:2308608/1‑149 (MQ=255)
                                                                                         gAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCt                                                           <  2:1350588/149‑1 (MQ=255)
                                                                                           aGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGt                                                         <  2:1602010/149‑1 (MQ=255)
                                                                                            gACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTa                                                        >  1:1944155/1‑149 (MQ=255)
                                                                                             aCATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTAt                                                       >  2:1262096/1‑149 (MQ=255)
                                                                                             aCATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTAt                                                       >  2:2135753/1‑149 (MQ=255)
                                                                                             aCATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTAt                                                       >  2:4081121/1‑149 (MQ=255)
                                                                                                ttAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGaa                                                    <  1:3372555/149‑1 (MQ=255)
                                                                                                ttAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCGGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGaa                                                    >  2:2506183/1‑149 (MQ=255)
                                                                                                  aaCGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAAt                                                  <  2:2278419/149‑1 (MQ=255)
                                                                                                   aCGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAAtt                                                 <  1:2756108/149‑1 (MQ=255)
                                                                                                    cGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATtc                                                >  1:1547473/1‑149 (MQ=255)
                                                                                                             aCGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACata                                       <  1:2210550/149‑1 (MQ=255)
                                                                                                              cGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACatat                                      <  1:2313222/149‑1 (MQ=255)
                                                                                                               gCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATc                                     <  1:4063224/149‑1 (MQ=255)
                                                                                                               gCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATc                                     <  1:3670758/149‑1 (MQ=255)
                                                                                                                 aCCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCaa                                   >  1:3049520/1‑149 (MQ=255)
                                                                                                                  ccGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCaaa                                  <  2:3618735/149‑1 (MQ=255)
                                                                                                                   cGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAAc                                 <  1:378159/149‑1 (MQ=255)
                                                                                                                    gAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAAcc                                >  1:831791/1‑149 (MQ=255)
                                                                                                                    gAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAAcc                                >  2:1166960/1‑149 (MQ=255)
                                                                                                                       cAGGTCAATCAGCAACTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGaa                             >  1:59994/1‑149 (MQ=255)
                                                                                                                        aGGTCAATCAGCAACTGATTGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAAc                            <  1:3287726/149‑1 (MQ=255)
                                                                                                                        aGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAAc                            >  1:1893427/1‑149 (MQ=255)
                                                                                                                        aGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAAc                            <  2:3384548/149‑1 (MQ=255)
                                                                                                                        aGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAAc                            <  2:3546428/149‑1 (MQ=255)
                                                                                                                          gTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTg                          <  2:3068411/149‑1 (MQ=255)
                                                                                                                                 aGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGgtg                   <  2:3712755/149‑1 (MQ=255)
                                                                                                                                    aaCTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTgg                >  1:3284981/1‑149 (MQ=255)
                                                                                                                                     aCTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGt               >  2:4061328/1‑149 (MQ=255)
                                                                                                                                      cTGATGGAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTa              <  2:2822243/149‑1 (MQ=255)
                                                                                                                                       tGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTAc             >  2:213727/1‑149 (MQ=255)
                                                                                                                                          tGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTCCGAg          >  2:3169420/1‑149 (MQ=255)
                                                                                                                                          tGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAg          <  1:2841232/149‑1 (MQ=255)
                                                                                                                                          tGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAg          >  2:3517871/1‑149 (MQ=255)
                                                                                                                                           ggAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgt         >  1:2920047/1‑149 (MQ=255)
                                                                                                                                            gAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgtg        >  1:3449318/1‑149 (MQ=255)
                                                                                                                                            gAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgtg        >  1:90539/1‑149 (MQ=255)
                                                                                                                                            gAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgtg        <  1:2035134/149‑1 (MQ=255)
                                                                                                                                            gAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgtg        <  1:485539/149‑1 (MQ=255)
                                                                                                                                            gAGCGCATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgtg        >  1:1027292/1‑149 (MQ=255)
                                                                                                                                             aGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgtgc       >  2:532772/1‑149 (MQ=255)
                                                                                                                                              gcgcATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgtgcg      <  2:1718687/149‑1 (MQ=255)
                                                                                                                                               cgcATAACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAgtgcgt     <  2:3881441/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ACCGCAATGGCGATGTTGTCTGGGCGTTGCTTGCCGTCTCACTGGTGCGCCACACGGATGGCACGCCGCTCTATTTTATCGCTCAGATTGAAGACATTAACGAGCTAAAACGCACCGAACAGGTCAATCAGCAACTGATGGAGCGCATCACTCTGGCTAACGAAGCGGGCGGGATTGGCATCTGGGAGTGGGAGTTGAAGCCGAATATTTTTAGCTGGGATAAGCGGATGTTCGAGCTGTATGAAATTCCTCCACATATCAAACCGAACTGGCAGGTGTGGTACGAGTGCGTGCT  >  CP009273/2137891‑2138185

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: