Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,962,962 T→A 100% intergenic (‑626/‑59) rppH ← / → mutH RNA pyrophosphohydrolase/methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,962,9620TA96.9% 826.2 / 9.0 253intergenic (‑626/‑59)rppH/mutHRNA pyrophosphohydrolase/methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base T (7/1);  new base A (131/114);  total (138/115)
Fisher's exact test for biased strand distribution p-value = 7.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 253 total aligned reads displayed.

CAGCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACATTGG  >  CP009273/2962816‑2963108
                                                                                                                                                  |                                                                                                                                                  
 aGCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATaca                                                                                                                                                >  2:2792586/1‑150 (MQ=255)
  gCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAATACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATacac                                                                                                                                               <  2:4813114/150‑1 (MQ=255)
    tttCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTg                                                                                                                                             <  2:4428906/150‑1 (MQ=255)
     gtCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGc                                                                                                                                            >  1:3658165/2‑150 (MQ=255)
      tCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCg                                                                                                                                           >  1:2274715/1‑150 (MQ=255)
       cTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGa                                                                                                                                          <  1:1125410/150‑1 (MQ=255)
        tgtgGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGaa                                                                                                                                         >  2:1721400/1‑150 (MQ=255)
         gtgGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAat                                                                                                                                        <  1:1544315/150‑1 (MQ=255)
         gtgGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAat                                                                                                                                        <  1:3989002/150‑1 (MQ=255)
         gtgGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAAGATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAat                                                                                                                                        <  2:1221999/150‑1 (MQ=255)
          tgGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAata                                                                                                                                       >  2:3535694/1‑150 (MQ=255)
           ggATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAatat                                                                                                                                      >  2:4028000/1‑150 (MQ=255)
            gATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATAtt                                                                                                                                     >  1:1253146/1‑150 (MQ=255)
            gATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATAtt                                                                                                                                     >  2:2770802/1‑150 (MQ=255)
               aaCATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCgg                                                                                                                                  <  1:990688/150‑1 (MQ=255)
                   tGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGcaca                                                                                                                              >  2:2579086/1‑150 (MQ=255)
                       gtAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAAt                                                                                                                          <  2:2420781/150‑1 (MQ=255)
                          aGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGc                                                                                                                       >  1:4144677/1‑150 (MQ=255)
                          aGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGc                                                                                                                       >  2:3009067/1‑150 (MQ=255)
                          aGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGc                                                                                                                       >  1:3242253/1‑150 (MQ=255)
                           gATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCt                                                                                                                      >  1:1094944/1‑150 (MQ=255)
                              ccTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCtgtt                                                                                                                   <  1:1594691/150‑1 (MQ=255)
                                tGTTTATTTTCAGTGACCAGATTTTGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCtgtttg                                                                                                                 <  2:2660658/150‑1 (MQ=255)
                                tGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCtgtttg                                                                                                                 <  1:2280535/150‑1 (MQ=255)
                                tGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCtgtttg                                                                                                                 <  2:2349014/150‑1 (MQ=255)
                                 gTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCtgtttgt                                                                                                                <  2:3763572/150‑1 (MQ=255)
                                    tATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGtttt                                                                                                             >  1:712466/1‑150 (MQ=255)
                                     aTTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGttttt                                                                                                            >  1:4996032/1‑150 (MQ=255)
                                      ttttCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGtttttt                                                                                                           <  1:2711814/150‑1 (MQ=255)
                                         tCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAAt                                                                                                        <  2:5052072/150‑1 (MQ=255)
                                              gACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAgg                                                                                                   <  1:955624/150‑1 (MQ=255)
                                                   gATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCa                                                                                              <  2:2281913/150‑1 (MQ=255)
                                                   gATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCa                                                                                              >  1:3884041/1‑150 (MQ=255)
                                                   gATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCa                                                                                              >  2:3251727/1‑150 (MQ=255)
                                                    aTTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCAt                                                                                             <  1:2179703/150‑1 (MQ=255)
                                                        ggAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACa                                                                                         <  1:919707/150‑1 (MQ=255)
                                                        ggAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACa                                                                                         <  1:4847099/150‑1 (MQ=255)
                                                          aaaaCCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATg                                                                                       >  1:4493400/1‑150 (MQ=255)
                                                            aaCCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTc                                                                                     >  1:3528293/1‑150 (MQ=255)
                                                             aCCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATACTTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTcc                                                                                    >  2:884010/1‑150 (MQ=255)
                                                             aCCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTcc                                                                                    <  2:643680/150‑1 (MQ=255)
                                                              cccGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTccc                                                                                   <  1:2976027/150‑1 (MQ=255)
                                                              cccGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTccc                                                                                   >  2:4856703/1‑150 (MQ=255)
                                                                  ttGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAAcc                                                                               >  1:194540/1‑150 (MQ=255)
                                                                   tGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGTCATGTCCCAAcca                                                                              >  2:4411877/1‑149 (MQ=255)
                                                                   tGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCt                                                                              <  1:3485396/150‑1 (MQ=255)
                                                                    gCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTc                                                                             <  1:3551790/150‑1 (MQ=255)
                                                                     cAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCg                                                                            >  2:3357858/1‑150 (MQ=255)
                                                                        tgtTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGccc                                                                         <  2:1673401/150‑1 (MQ=255)
                                                                         gtTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGccaa                                                                        <  2:5140386/150‑3 (MQ=255)
                                                                         gtTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCa                                                                        <  2:3918905/150‑1 (MQ=255)
                                                                             cgcAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGc                                                                    <  2:1801706/150‑1 (MQ=255)
                                                                              gcAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGct                                                                   >  2:157744/1‑150 (MQ=255)
                                                                               cAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGctc                                                                  >  2:739037/1‑150 (MQ=255)
                                                                                   tCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCtctctc                                                              >  1:555237/1‑150 (MQ=255)
                                                                                    cGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTctcc                                                             >  1:2657336/1‑150 (MQ=255)
                                                                                    cGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTctcc                                                             >  1:5132315/1‑150 (MQ=255)
                                                                                     gATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACTCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTctcct                                                            <  2:2213443/150‑1 (MQ=255)
                                                                                     gATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTctcct                                                            <  1:1769450/150‑1 (MQ=255)
                                                                                      aTTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTctcctc                                                           >  2:3399027/1‑150 (MQ=255)
                                                                                       ttACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCtcctcc                                                          <  1:528468/150‑1 (MQ=255)
                                                                                       ttACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCtcctcc                                                          >  1:680394/1‑150 (MQ=255)
                                                                                        tACCGGCAACCTAAAGAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATACTTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTccc                                                         >  1:2332299/1‑150 (MQ=255)
                                                                                        tACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTccc                                                         <  2:4425383/150‑1 (MQ=255)
                                                                                        tACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTccc                                                         >  2:2933987/1‑150 (MQ=255)
                                                                                        tACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTccc                                                         <  1:2824140/150‑1 (MQ=255)
                                                                                        tACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTccc                                                         >  1:4340385/1‑150 (MQ=255)
                                                                                         aCCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCg                                                        >  2:3566769/1‑150 (MQ=255)
                                                                                          ccGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGa                                                       >  1:1397673/1‑150 (MQ=255)
                                                                                           cGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAGCGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGaa                                                      >  1:3815072/1‑150 (MQ=255)
                                                                                           cGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGaa                                                      >  1:854779/1‑150 (MQ=255)
                                                                                           cGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGaa                                                      <  2:2006000/150‑1 (MQ=255)
                                                                                             gCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAAc                                                    >  1:428062/1‑150 (MQ=255)
                                                                                             gCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAAc                                                    <  1:1698819/150‑1 (MQ=255)
                                                                                                 ccTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTgaa                                                >  1:4426546/1‑150 (MQ=255)
                                                                                                    aaaaaGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTgaagaa                                             >  2:2178904/1‑150 (MQ=255)
                                                                                                         gCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgtt                                        >  2:815636/1‑150 (MQ=255)
                                                                                                         gCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgtt                                        <  1:2931136/150‑1 (MQ=255)
                                                                                                         gCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgtt                                        <  1:287268/150‑1 (MQ=255)
                                                                                                         gCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgtt                                        >  2:1188413/1‑150 (MQ=255)
                                                                                                          cAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgttg                                       <  2:343775/150‑1 (MQ=255)
                                                                                                          cAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgttg                                       <  1:3419309/150‑1 (MQ=255)
                                                                                                            atatAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgttgtt                                     >  2:4199830/1‑150 (MQ=255)
                                                                                                             tataAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTa                                    >  1:5141247/1‑150 (MQ=255)
                                                                                                             tataAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTa                                    >  2:3738692/1‑150 (MQ=255)
                                                                                                                       agaATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGcac                          >  1:592844/1‑150 (MQ=255)
                                                                                                                         aaTTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAg                        <  2:3469865/150‑1 (MQ=255)
                                                                                                                         aaTTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAg                        >  2:61420/1‑150 (MQ=255)
                                                                                                                         aaTTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAg                        >  2:1378814/1‑150 (MQ=255)
                                                                                                                            tGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCaa                     >  1:1463913/1‑150 (MQ=255)
                                                                                                                            tGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCaa                     >  1:1797553/1‑150 (MQ=255)
                                                                                                                             gAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAAc                    <  2:1094944/150‑1 (MQ=255)
                                                                                                                                caacGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACttt                 <  2:1434686/150‑1 (MQ=255)
                                                                                                                                 aacGCATGTGGAAAAAGATACACTGCGAATATTCGGCCCATACTTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTc                >  2:3152030/1‑150 (MQ=255)
                                                                                                                                 aacGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTc                <  1:1052123/150‑1 (MQ=255)
                                                                                                                                    gCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTgg             >  1:2424490/1‑150 (MQ=255)
                                                                                                                                    gCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTgg             >  1:4696669/1‑150 (MQ=255)
                                                                                                                                     cATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGt            >  1:4766372/1‑150 (MQ=255)
                                                                                                                                      aTGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGtt           >  1:3928710/1‑150 (MQ=255)
                                                                                                                                               aaGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACATTgg  <  2:465479/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CAGCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACATTGG  >  CP009273/2962816‑2963108

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: