Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 277,970 G→A 6.4% intergenic (‑276/‑19) yagA ← / → yagE CP4‑6 prophage; putative DNA‑binding transcriptional regulator/2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273277,9700GA 6.4% 1040.0 / 35.3 298intergenic (‑276/‑19)yagA/yagECP4‑6 prophage; putative DNA‑binding transcriptional regulator/2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base G (147/131);  new base A (4/15);  total (152/146)
Fisher's exact test for biased strand distribution p-value = 8.44e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 300 total aligned reads displayed.

CTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCT  >  CP009273/277821‑278119
                                                                                                                                                     |                                                                                                                                                     
  ttGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGAt                                                                                                                                                     <  1:378947/150‑1 (MQ=255)
   tGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGAtt                                                                                                                                                    >  1:4319740/1‑150 (MQ=255)
     gggAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCGTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATAATTCa                                                                                                                                                  <  2:3038540/150‑1 (MQ=255)
     gggAGAGGGCCGGGGTGAGGGGAAATATGTGATCCAGATTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCa                                                                                                                                                  <  1:4195514/150‑1 (MQ=255)
      ggAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAg                                                                                                                                                 <  2:2824355/150‑1 (MQ=255)
           gggCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAgg                                                                                                                                            >  1:4743732/1‑150 (MQ=255)
            ggCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGga                                                                                                                                           <  1:3471977/150‑1 (MQ=255)
             gCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGgag                                                                                                                                          <  2:3004782/150‑1 (MQ=255)
             gCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGgag                                                                                                                                          <  2:3564866/150‑1 (MQ=255)
             gCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCGTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGgag                                                                                                                                          <  1:3859334/150‑1 (MQ=255)
                  ggTGAGGGGCAATATGTGATCCAGCGTAAATTTCCCGCACTCCCACTTCCCTTCCGATTTCCCTCTCCATGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTCGAACAAAATGATTCAGCAAGGAGAtctc                                                                                                                                     <  2:4398711/150‑1 (MQ=255)
                     gAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATg                                                                                                                                  <  2:1664062/150‑1 (MQ=255)
                      aGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGTTTTACCACTCCTTGTTCTGCGTCTTAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGc                                                                                                                                 <  1:2370791/150‑1 (MQ=255)
                      aGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATAATTCAGCAAGGAGATCTCATGc                                                                                                                                 <  2:776539/150‑1 (MQ=255)
                        gggCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCg                                                                                                                               >  1:138787/1‑150 (MQ=255)
                        gggCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCg                                                                                                                               >  1:141988/1‑150 (MQ=255)
                        gggCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCg                                                                                                                               <  2:4604075/150‑1 (MQ=255)
                          gCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCa                                                                                                                             >  1:2031215/1‑150 (MQ=255)
                              tatGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTcc                                                                                                                         <  2:3644855/150‑1 (MQ=255)
                              tatGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTcc                                                                                                                         <  2:2458620/150‑1 (MQ=255)
                                  tgATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCgt                                                                                                                     >  2:1240542/1‑150 (MQ=255)
                                  tgATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCgt                                                                                                                     >  1:560855/1‑150 (MQ=255)
                                  tgATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCgt                                                                                                                     >  2:3461735/1‑150 (MQ=255)
                                    aTCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCgttg                                                                                                                   <  2:5285802/150‑1 (MQ=255)
                                       cAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTc                                                                                                                <  2:879939/150‑1 (MQ=255)
                                        aGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAGATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCa                                                                                                               >  2:1794349/1‑150 (MQ=255)
                                         gCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAc                                                                                                              <  1:2703894/150‑1 (MQ=255)
                                          cTTAAATTTCCCGCACTCCCTCTTCCCTTCNGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACg                                                                                                             <  2:4272011/150‑1 (MQ=255)
                                          cTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACg                                                                                                             >  2:3784283/1‑150 (MQ=255)
                                              aaTTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGaa                                                                                                         <  2:4258949/150‑1 (MQ=255)
                                                tttCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGGTGTTCACGGGCATc                                                                                                       >  1:636057/1‑150 (MQ=255)
                                                 ttCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCa                                                                                                      <  1:2309417/150‑1 (MQ=255)
                                                   cccGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCAtt                                                                                                    <  2:4927887/150‑1 (MQ=255)
                                                           cccTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGt                                                                                            >  1:3745920/1‑150 (MQ=255)
                                                             ctctTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGtct                                                                                          >  2:1389939/1‑150 (MQ=255)
                                                             ctctTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGtct                                                                                          >  1:4711098/1‑150 (MQ=255)
                                                             ctctTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGtct                                                                                          >  1:4541770/1‑150 (MQ=255)
                                                              tctTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGtctc                                                                                         <  2:1718315/150‑1 (MQ=255)
                                                               cttcccttccGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTcc                                                                                        >  1:102076/1‑150 (MQ=255)
                                                                ttcccttccGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGACTCATTCCCCCTGTCTcca                                                                                       >  1:370171/1‑150 (MQ=255)
                                                                ttcccttccGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTcca                                                                                       <  2:1453948/150‑1 (MQ=255)
                                                                   ccttccGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTccacca                                                                                    <  2:5005134/150‑1 (MQ=255)
                                                                        cGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCAttttt                                                                               >  2:4868437/1‑150 (MQ=255)
                                                                         gATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGACTCATTCCCCCTGTCTCCACCATGTTTa                                                                              >  1:772336/1‑150 (MQ=255)
                                                                               ccTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTAccgccg                                                                        >  1:1643107/1‑150 (MQ=255)
                                                                                ctctCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGa                                                                       >  2:4105338/1‑150 (MQ=255)
                                                                                ctctCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGa                                                                       <  1:1888011/150‑1 (MQ=255)
                                                                                  ctcCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACg                                                                     <  1:1705674/150‑1 (MQ=255)
                                                                                        gTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGc                                                               >  2:3352408/1‑150 (MQ=255)
                                                                                         ttCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATATAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCt                                                              <  2:3943185/150‑1 (MQ=255)
                                                                                           cTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCg                                                            <  2:4944469/150‑1 (MQ=255)
                                                                                              cGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATa                                                         <  1:3363225/150‑1 (MQ=255)
                                                                                              cGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATa                                                         <  1:1135720/150‑1 (MQ=255)
                                                                                              cGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATa                                                         >  1:3050333/1‑150 (MQ=255)
                                                                                               gTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATaa                                                        >  1:2184700/1‑150 (MQ=255)
                                                                                               gTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATaa                                                        <  1:1082451/150‑1 (MQ=255)
                                                                                               gTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATaa                                                        <  2:1523251/150‑1 (MQ=255)
                                                                                               gTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATaa                                                        >  1:3930583/1‑150 (MQ=255)
                                                                                                  aTAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGcc                                                     <  2:3134250/150‑1 (MQ=255)
                                                                                                  aTAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGcc                                                     >  1:2532369/1‑150 (MQ=255)
                                                                                                    aGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCgg                                                   >  1:922171/1‑150 (MQ=255)
                                                                                                     gTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGACTCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCggg                                                  >  1:1262702/1‑150 (MQ=255)
                                                                                                        tGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCAc                                               >  2:621873/1‑150 (MQ=255)
                                                                                                          aTCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCAccg                                             >  1:1987308/1‑150 (MQ=255)
                                                                                                          aTCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCAccg                                             <  2:2085386/150‑1 (MQ=255)
                                                                                                           tCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCAccgc                                            >  2:6759/1‑150 (MQ=255)
                                                                                                           tCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCAAGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGACTCCACCATATTAACCGCCGACGCCCAGCTCGATAAGCCGGGCAccgc                                            >  1:3095288/1‑150 (MQ=255)
                                                                                                            cGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCCccgcc                                           >  2:4452054/1‑150 (MQ=255)
                                                                                                                 aataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCt                                      >  2:2335652/1‑150 (MQ=255)
                                                                                                                  ataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTg                                     >  2:2188201/1‑150 (MQ=255)
                                                                                                                   taaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGa                                    >  1:2964497/1‑150 (MQ=255)
                                                                                                                     aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATc                                  <  2:102115/150‑1 (MQ=255)
                                                                                                                     aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATc                                  <  2:3271034/150‑1 (MQ=255)
                                                                                                                      acgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTAACGCCGACGGCAAGCTCGATAATCCGGGCACCGCCGCGCTGAtag                                 >  1:2878396/1‑148 (MQ=255)
                                                                                                                       cgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATcga                                >  1:4048263/1‑150 (MQ=255)
                                                                                                                           cgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGAt                            >  2:5152957/1‑150 (MQ=255)
                                                                                                                             cTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCt                          <  1:2268818/150‑1 (MQ=255)
                                                                                                                             cTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCt                          >  1:467268/1‑150 (MQ=255)
                                                                                                                             cTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCt                          <  1:2269057/150‑1 (MQ=255)
                                                                                                                              tGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTg                         <  2:4319740/150‑1 (MQ=255)
                                                                                                                                 tCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATc                      >  1:3508909/1‑150 (MQ=255)
                                                                                                                                  cTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCa                     >  1:1493316/1‑150 (MQ=255)
                                                                                                                                    ataATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCaaa                   >  2:1944188/1‑150 (MQ=255)
                                                                                                                                     taATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAg                  >  1:2107797/1‑150 (MQ=255)
                                                                                                                                      aaTGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGc                 >  1:673953/1‑150 (MQ=255)
                                                                                                                                      aaTGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGc                 >  2:3199265/1‑150 (MQ=255)
                                                                                                                                      aaTGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGc                 >  2:3849554/1‑150 (MQ=255)
                                                                                                                                       aTGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCa                >  1:4251898/1‑150 (MQ=255)
                                                                                                                                         gTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAgg              <  2:284919/150‑1 (MQ=255)
                                                                                                                                         gTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTAACAGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAgt              >  2:1900318/1‑149 (MQ=255)
                                                                                                                                             aaCAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGtt          >  1:3991586/1‑150 (MQ=255)
                                                                                                                                              aCAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTg         >  1:1742499/1‑150 (MQ=255)
                                                                                                                                              aCAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTg         <  2:2665060/150‑1 (MQ=255)
                                                                                                                                              aCAAAATAATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTg         <  2:5316915/150‑1 (MQ=255)
                                                                                                                                                 aaaTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCCCCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACg      >  2:2593547/1‑150 (MQ=255)
                                                                                                                                                 aaaTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACg      <  2:3595920/150‑1 (MQ=255)
                                                                                                                                                  aaTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACgg     <  1:5260006/150‑1 (MQ=255)
                                                                                                                                                   aTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGc    >  2:736485/1‑150 (MQ=255)
                                                                                                                                                    tGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGcc   <  2:560855/150‑1 (MQ=255)
                                                                                                                                                     gATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCt  <  2:4785099/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTTGGGGAGAGGGCCGGGGTGAGGGGCAATATGTGATCCAGCTTAAATTTCCCGCACTCCCTCTTCCCTTCCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCT  >  CP009273/277821‑278119

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: