Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,754,840 T→A 10.9% E687V (GAA→GTA)  yfjK ← DEAD/H helicase‑like protein, CP4‑57 putative defective prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,754,8400TA10.9% 895.1 / 76.5 282E687V (GAA→GTA) yfjKDEAD/H helicase‑like protein, CP4‑57 putative defective prophage
Reads supporting (aligned to +/- strand):  ref base T (151/100);  new base A (19/12);  total (170/112)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 287 total aligned reads displayed.

GTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGG  >  CP009273/2754691‑2754989
                                                                                                                                                     |                                                                                                                                                     
   gACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCa                                                                                                                                                    >  1:2086420/1‑150 (MQ=255)
   gACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCa                                                                                                                                                    >  2:3159440/1‑150 (MQ=255)
    aCAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAg                                                                                                                                                   <  1:3736340/150‑1 (MQ=255)
     cAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGc                                                                                                                                                  >  2:4517833/1‑150 (MQ=255)
        aCGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAg                                                                                                                                               <  1:3274382/150‑1 (MQ=255)
              cTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTAcc                                                                                                                                         <  1:3059852/150‑1 (MQ=255)
               tAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCa                                                                                                                                        <  1:2783548/150‑1 (MQ=255)
                  aTCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTc                                                                                                                                     <  1:1674853/150‑1 (MQ=255)
                   tCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCg                                                                                                                                    <  1:2822108/150‑1 (MQ=255)
                    cAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCgg                                                                                                                                   >  1:3850240/1‑150 (MQ=255)
                     aaTGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGgt                                                                                                                                  >  1:3489230/1‑150 (MQ=255)
                      aTGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAACTGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGgtg                                                                                                                                 <  2:1184790/150‑1 (MQ=255)
                      aTGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGgtg                                                                                                                                 >  1:1982176/1‑150 (MQ=255)
                      aTGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGgtg                                                                                                                                 >  1:3927124/1‑150 (MQ=255)
                        gtcgtcTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTa                                                                                                                               >  1:3292034/1‑150 (MQ=255)
                             cTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTc                                                                                                                          >  1:4406417/1‑150 (MQ=255)
                              tAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCt                                                                                                                         >  1:523562/1‑150 (MQ=255)
                              tAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCt                                                                                                                         >  2:1393329/1‑150 (MQ=255)
                               aaTACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTg                                                                                                                        >  1:886485/1‑150 (MQ=255)
                               aaTACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTg                                                                                                                        >  2:1692928/1‑150 (MQ=255)
                               aaTACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTg                                                                                                                        >  2:2668503/1‑150 (MQ=255)
                               aaTACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTg                                                                                                                        >  1:4589837/1‑150 (MQ=255)
                                 tACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGaa                                                                                                                      >  2:297412/1‑150 (MQ=255)
                                  aCCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAAt                                                                                                                     <  2:1781817/150‑1 (MQ=255)
                                  aCCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAAt                                                                                                                     <  2:3018485/150‑1 (MQ=255)
                                  aCCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAAt                                                                                                                     <  2:509216/150‑1 (MQ=255)
                                   ccGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATg                                                                                                                    <  2:1610849/150‑1 (MQ=255)
                                    cGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATgg                                                                                                                   <  2:1450564/150‑1 (MQ=255)
                                         gcgcgcCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATg                                                                                                              <  1:4380986/150‑1 (MQ=255)
                                          cgcgcCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGc                                                                                                             >  2:4330950/1‑150 (MQ=255)
                                             gcCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCaa                                                                                                          >  1:3222052/1‑150 (MQ=255)
                                               cTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAAtc                                                                                                        >  2:4094079/1‑150 (MQ=255)
                                                 gATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAAtctc                                                                                                      >  1:1725093/1‑150 (MQ=255)
                                                   tAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCtctt                                                                                                    >  1:2327252/1‑150 (MQ=255)
                                                       cATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTaa                                                                                                >  1:2648736/1‑150 (MQ=255)
                                                        aTCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAg                                                                                               >  1:660161/1‑150 (MQ=255)
                                                         tCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAgg                                                                                              >  1:709492/1‑150 (MQ=255)
                                                         tCTGATCAACGTAAACAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCAACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACAAGTCGGTGTACCGTCTGAAGGGGAATGCAAATCTCATCTAAgg                                                                                              >  2:64349/1‑150 (MQ=255)
                                                          cTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGc                                                                                             >  1:1944360/1‑150 (MQ=255)
                                                           tGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCa                                                                                            >  1:4027056/1‑150 (MQ=255)
                                                              tCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTg                                                                                         <  1:632073/150‑1 (MQ=255)
                                                                aaCGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGac                                                                                       >  1:3116368/1‑150 (MQ=255)
                                                                aaCGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGac                                                                                       >  2:2412676/1‑150 (MQ=255)
                                                                aaCGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGac                                                                                       >  2:3846032/1‑150 (MQ=255)
                                                                  cGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGacac                                                                                     >  1:299445/1‑150 (MQ=255)
                                                                  cGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGacac                                                                                     <  1:853190/150‑1 (MQ=255)
                                                                      cccAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGACTGCCAATCTCTTCTAAGGCATTGACACCtgg                                                                                 >  2:3650840/1‑150 (MQ=255)
                                                                      cccAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCtgg                                                                                 <  2:2729134/150‑1 (MQ=255)
                                                                         aaTCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCACGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCtggtgg                                                                              >  1:122021/1‑150 (MQ=255)
                                                                             gATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGa                                                                          >  1:2111893/1‑150 (MQ=255)
                                                                              aTCTACTCCAACTGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGAt                                                                         <  2:1924183/150‑1 (MQ=255)
                                                                                cTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATgg                                                                       >  2:1012515/1‑150 (MQ=255)
                                                                                cTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATgg                                                                       >  2:2146782/1‑150 (MQ=255)
                                                                                  aCTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGaa                                                                     <  1:1471504/150‑1 (MQ=255)
                                                                                  aCTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGaa                                                                     >  2:76369/1‑150 (MQ=255)
                                                                                   cTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACCTCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGaaa                                                                    <  2:4023474/150‑1 (MQ=255)
                                                                                    tCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAAc                                                                   >  2:3064882/1‑150 (MQ=255)
                                                                                     ccAACTGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACt                                                                  <  1:2925187/150‑1 (MQ=255)
                                                                                        actGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTc                                                               <  2:2376515/148‑1 (MQ=255)
                                                                                        aaTGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTc                                                               <  1:2404914/150‑1 (MQ=255)
                                                                                        aaTGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTc                                                               <  1:4158778/150‑1 (MQ=255)
                                                                                            tCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAAc                                                           <  1:3422594/150‑1 (MQ=255)
                                                                                            tCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAAc                                                           <  2:2861998/150‑1 (MQ=255)
                                                                                            tCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAAc                                                           <  1:3178853/150‑1 (MQ=255)
                                                                                             ccTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACg                                                          >  2:2573117/1‑150 (MQ=255)
                                                                                               tGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGcc                                                        <  1:3683469/150‑1 (MQ=255)
                                                                                                gAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCa                                                       >  1:3696852/1‑150 (MQ=255)
                                                                                                gAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGGGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCa                                                       >  1:568815/1‑150 (MQ=255)
                                                                                                 aGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACTCCaa                                                      >  2:1598706/1‑150 (MQ=255)
                                                                                                 aGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCaa                                                      >  2:4361689/1‑150 (MQ=255)
                                                                                                 aGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCaa                                                      >  2:1164074/1‑150 (MQ=255)
                                                                                                     tCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTg                                                  >  2:3339160/1‑150 (MQ=255)
                                                                                                     tCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTg                                                  <  1:2567401/150‑1 (MQ=255)
                                                                                                       ccGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGAc                                                >  2:22095/1‑150 (MQ=255)
                                                                                                         gAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACAt                                              >  2:2666245/1‑150 (MQ=255)
                                                                                                         gAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACAt                                              >  2:616214/1‑150 (MQ=255)
                                                                                                         gAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACAt                                              <  1:3598777/150‑1 (MQ=255)
                                                                                                             aTACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACAtgtgt                                          >  1:2192893/1‑150 (MQ=255)
                                                                                                                cTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTa                                       >  1:1787730/1‑150 (MQ=255)
                                                                                                                 tGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTaa                                      >  2:1078663/1‑150 (MQ=255)
                                                                                                                  ggCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAAt                                     >  1:4625664/1‑150 (MQ=255)
                                                                                                                   gCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATc                                    >  1:1371426/1‑150 (MQ=255)
                                                                                                                       aaCTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACt                                >  2:3904716/1‑150 (MQ=255)
                                                                                                                         cTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACttt                              >  1:1302782/1‑150 (MQ=255)
                                                                                                                         cTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACttt                              <  2:3491130/150‑1 (MQ=255)
                                                                                                                          tCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTtg                             <  2:2293243/150‑1 (MQ=255)
                                                                                                                           cATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTACTCGACTTtgc                            >  2:2897766/1‑150 (MQ=255)
                                                                                                                            aTCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTtgct                           <  2:3250836/150‑1 (MQ=255)
                                                                                                                            aTCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTtgct                           <  2:2889342/150‑1 (MQ=255)
                                                                                                                                 caTCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGta                      >  2:3856284/1‑150 (MQ=255)
                                                                                                                                 caTCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGta                      <  1:4306230/150‑1 (MQ=255)
                                                                                                                                  aTCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGtat                     <  2:3008029/150‑1 (MQ=255)
                                                                                                                                  aTCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGtat                     <  2:2045800/150‑1 (MQ=255)
                                                                                                                                   tCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGtata                    <  2:831512/150‑1 (MQ=255)
                                                                                                                                    cAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAc                   >  1:791674/1‑150 (MQ=255)
                                                                                                                                            tcgGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCttt           >  2:4481069/1‑150 (MQ=255)
                                                                                                                                             cgGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTa          >  2:4193118/1‑150 (MQ=255)
                                                                                                                                                aaaaTTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAg       >  1:2701994/1‑150 (MQ=255)
                                                                                                                                                     tCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtgg  <  2:4479164/150‑1 (MQ=255)
                                                                                                                                                     tCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtgg  >  1:758396/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGG  >  CP009273/2754691‑2754989

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: