Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,962,962 T→A 45.1% 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,9620TA45.1% 134.8 / inf 355intergenic (‑626/‑59)rppH/mutHRNA pyrophosphohydrolase/methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base T (114/81);  new base A (84/76);  total (198/157)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

Only 100 of 356 total aligned reads displayed.

CCCAGCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACATTGGGA  >  CP009273/2962814‑2963110
                                                                                                                                                    |                                                                                                                                                    
   aGCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTaca                                                                                                                                                  >  1:74723/1‑150 (MQ=255)
        tCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGACAAAGTTACACTGCg                                                                                                                                             >  2:3264524/1‑150 (MQ=255)
           gtgGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAat                                                                                                                                          <  1:927138/150‑1 (MQ=255)
           gtgGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAat                                                                                                                                          >  1:1881012/1‑150 (MQ=255)
              gATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATAtt                                                                                                                                       >  1:2573040/1‑150 (MQ=255)
               aTAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTc                                                                                                                                      >  2:2059653/1‑150 (MQ=255)
               aTAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTc                                                                                                                                      >  2:1040975/1‑150 (MQ=255)
               aTAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTc                                                                                                                                      >  1:2119341/1‑150 (MQ=255)
                tAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCg                                                                                                                                     >  2:4065072/1‑150 (MQ=255)
                tAACAGGGTGTTAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATAGACATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCg                                                                                                                                     <  1:304108/150‑1 (MQ=255)
                    aTGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCTGcac                                                                                                                                 >  2:4269725/1‑150 (MQ=255)
                      ggTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACAt                                                                                                                               <  2:509912/150‑1 (MQ=255)
                          tAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAAtt                                                                                                                           <  1:2793849/150‑1 (MQ=255)
                          tAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAAtt                                                                                                                           <  2:4685888/150‑1 (MQ=255)
                             gATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCt                                                                                                                        >  2:617650/1‑150 (MQ=255)
                             gATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCt                                                                                                                        >  2:4998226/1‑150 (MQ=255)
                               tCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCtgt                                                                                                                      <  2:4020722/150‑1 (MQ=255)
                                 cTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCtgttt                                                                                                                    >  2:5127240/1‑150 (MQ=255)
                                     ttATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTgtttgttt                                                                                                                <  1:3814874/150‑1 (MQ=255)
                                       aTTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGttttt                                                                                                              <  2:4338477/150‑1 (MQ=255)
                                        ttttCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGtttttt                                                                                                             <  2:290207/150‑1 (MQ=255)
                                          ttCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTaa                                                                                                           >  2:2862939/1‑150 (MQ=255)
                                           tCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAAt                                                                                                          >  1:1969675/1‑150 (MQ=255)
                                           tCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAAt                                                                                                          >  1:4546881/1‑150 (MQ=255)
                                             aGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCa                                                                                                        >  2:721132/1‑150 (MQ=255)
                                             aGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCa                                                                                                        >  1:3792299/1‑150 (MQ=255)
                                              gTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAACTCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCaa                                                                                                       <  2:1726986/150‑1 (MQ=255)
                                               tGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAg                                                                                                      >  1:3514000/1‑150 (MQ=255)
                                               tGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAg                                                                                                      >  1:1626820/1‑150 (MQ=255)
                                               tGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAg                                                                                                      >  2:3311960/1‑150 (MQ=255)
                                                gACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAgg                                                                                                     >  1:2353512/1‑150 (MQ=255)
                                                    aGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATc                                                                                                 >  1:4567767/1‑150 (MQ=255)
                                                    aGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATc                                                                                                 >  1:4697671/1‑150 (MQ=255)
                                                     gATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAATCAATATAAATCAGAGAATTGAACAACTCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGATTGTTTTTTAATCAAGGTATCa                                                                                                <  2:150113/150‑1 (MQ=255)
                                                         tGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGAc                                                                                            <  1:1744924/150‑1 (MQ=255)
                                                              aaCCCGTTGCAGTGTTGAGCAACTCGATTACCGGCAACCTAAAAAGCAAGATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTc                                                                                       <  2:5055565/150‑1 (MQ=255)
                                                               aCCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTcc                                                                                      <  1:404448/150‑1 (MQ=255)
                                                                 ccGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCa                                                                                    <  2:2405041/150‑1 (MQ=255)
                                                                 ccGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCa                                                                                    >  1:3382626/1‑150 (MQ=255)
                                                                 ccGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCa                                                                                    >  1:4942898/1‑150 (MQ=255)
                                                                        aGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGc                                                                             >  2:295378/1‑150 (MQ=255)
                                                                           gtTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCa                                                                          <  1:131865/150‑1 (MQ=255)
                                                                             tGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACt                                                                        <  1:1125790/150‑1 (MQ=255)
                                                                             tGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACt                                                                        >  2:2452930/1‑150 (MQ=255)
                                                                               cgcAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGc                                                                      >  1:2964677/1‑150 (MQ=255)
                                                                                 cAACTCGATTACCGGCAACCTAACAAGCAATATAAATCAGAGAATTGAACAACGCAGGTGGAAAAAGTTACACTGCGAATATTCGGCACATAAGGGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGctc                                                                    <  1:3877081/150‑1 (MQ=255)
                                                                                 cAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGctc                                                                    <  1:3996947/150‑1 (MQ=255)
                                                                                   aCTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGctctc                                                                  >  1:2019079/1‑150 (MQ=255)
                                                                                    cTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGctctct                                                                 >  2:10660/1‑150 (MQ=255)
                                                                                    cTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGctctct                                                                 <  2:2205013/150‑1 (MQ=255)
                                                                                       gATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTctcct                                                              <  1:3238820/150‑1 (MQ=255)
                                                                                       gATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTctcct                                                              >  2:4172496/1‑150 (MQ=255)
                                                                                        aTTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTctcctc                                                             >  2:4032915/1‑150 (MQ=255)
                                                                                          tACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTccc                                                           <  1:2265288/150‑1 (MQ=255)
                                                                                             cGGCAACCTAAAAAGCAATATAAATCCGAGAATTGAACAACGCAGGGGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGaa                                                        >  2:1295021/1‑150 (MQ=255)
                                                                                              ggCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGaaa                                                       >  2:4188254/1‑150 (MQ=255)
                                                                                              ggCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGaaa                                                       <  2:456905/150‑1 (MQ=255)
                                                                                              ggCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATAATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGaaa                                                       >  1:4331330/1‑150 (MQ=255)
                                                                                               gCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAAc                                                      <  2:878752/150‑1 (MQ=255)
                                                                                                cAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACt                                                     >  2:2565205/1‑150 (MQ=255)
                                                                                                     tAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTgaaga                                                <  2:1880989/150‑1 (MQ=255)
                                                                                                      aaaaaGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTgaagaa                                               >  2:4898820/1‑150 (MQ=255)
                                                                                                       aaaaGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAAc                                              >  1:266126/1‑150 (MQ=255)
                                                                                                          aGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgt                                           <  1:1712946/150‑1 (MQ=255)
                                                                                                           gCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgtt                                          >  1:1863871/1‑150 (MQ=255)
                                                                                                            cAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgttg                                         >  1:962092/1‑150 (MQ=255)
                                                                                                             aaTATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgttgt                                        >  1:1025010/1‑150 (MQ=255)
                                                                                                             aaTATAAATCAGAGAATTGAACAACGCATGTCGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgttgt                                        >  2:3339436/1‑150 (MQ=255)
                                                                                                              atatAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgttgtt                                       <  2:367960/150‑1 (MQ=255)
                                                                                                               tataAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTa                                      <  1:5088320/150‑1 (MQ=255)
                                                                                                                ataAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAg                                     >  2:2512481/1‑150 (MQ=255)
                                                                                                                 taAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAgc                                    >  2:1904150/1‑150 (MQ=255)
                                                                                                                   aaTCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAgcgc                                  >  1:1409590/1‑150 (MQ=255)
                                                                                                                   aaTCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAgcgc                                  >  1:3477778/1‑150 (MQ=255)
                                                                                                                   aaTCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAgcgc                                  >  1:901966/1‑150 (MQ=255)
                                                                                                                    aTCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCa                                 >  2:2718882/1‑150 (MQ=255)
                                                                                                                      cAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAg                               <  2:1546778/150‑1 (MQ=255)
                                                                                                                      cAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAg                               >  2:2780057/1‑150 (MQ=255)
                                                                                                                      cAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAg                               <  1:1242877/150‑1 (MQ=255)
                                                                                                                       agagAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGc                              >  1:149251/1‑150 (MQ=255)
                                                                                                                         agaATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGcac                            <  1:830666/150‑1 (MQ=255)
                                                                                                                          gaATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGcaca                           <  1:3261807/150‑1 (MQ=255)
                                                                                                                           aaTTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAg                          >  2:1675537/1‑150 (MQ=255)
                                                                                                                               gAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAAc                      <  1:595162/150‑1 (MQ=255)
                                                                                                                                 acaacGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACtt                    >  2:3742067/1‑150 (MQ=255)
                                                                                                                                 acaacGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACtt                    <  1:662818/150‑1 (MQ=255)
                                                                                                                                  caacGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACttt                   >  2:662936/1‑150 (MQ=255)
                                                                                                                                   aacGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATAATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTc                  >  2:3569652/1‑150 (MQ=255)
                                                                                                                                       cATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGt              >  2:3046556/1‑150 (MQ=255)
                                                                                                                                         tgtgGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTta            >  2:1410324/1‑150 (MQ=255)
                                                                                                                                          gtgGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTtat           >  2:2112002/1‑150 (MQ=255)
                                                                                                                                          gtgGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCAGTATGATTTTTAATCAAGGTATCATGACATGTACCAACCTCGCCCACTGCTCTCTCCTCCCGAAAATGAAGAACAGTTGATAGAGCAATCACAGCAACTTTCTCGTtaa           >  2:1814676/1‑149 (MQ=255)
                                                                                                                                           ttgAAAAAGATACACTGCGAATATTCGGCTCATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTtata          <  1:3572017/148‑1 (MQ=255)
                                                                                                                                           tgGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTtata          <  2:556357/150‑1 (MQ=255)
                                                                                                                                           tgGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTtata          <  1:4154675/150‑1 (MQ=255)
                                                                                                                                             gAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACa        >  2:4724213/1‑150 (MQ=255)
                                                                                                                                              aaaaaGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACAt       >  2:2321494/1‑150 (MQ=255)
                                                                                                                                              aaaaaGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACAt       >  1:3362472/1‑150 (MQ=255)
                                                                                                                                                 aaGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACATTgg    >  2:4628948/1‑150 (MQ=255)
                                                                                                                                                   gTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACATTGGga  >  1:2187023/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCCAGCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACATTGGGA  >  CP009273/2962814‑2963110

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: