Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,612,781 C→A 31.3% A399E (GCG→GAG)  yjjJ → putative protein kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,612,7810CA31.3% 571.2 / 279.2 313A399E (GCG→GAG) yjjJputative protein kinase
Reads supporting (aligned to +/- strand):  ref base C (116/99);  new base A (60/38);  total (176/137)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 322 total aligned reads displayed.

GTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACAGCC  >  CP009273/4612632‑4612930
                                                                                                                                                     |                                                                                                                                                     
gTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGc                                                                                                                                                       >  2:1307910/1‑150 (MQ=255)
gTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGc                                                                                                                                                       >  2:4073820/1‑150 (MQ=255)
gTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGc                                                                                                                                                       >  2:272426/1‑150 (MQ=255)
  tttatttaTCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGa                                                                                                                                                     >  2:2177194/1‑150 (MQ=255)
      tttatCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATccc                                                                                                                                                 <  2:4105084/150‑1 (MQ=255)
      tttatCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATccc                                                                                                                                                 <  1:3077674/150‑1 (MQ=255)
      tttatCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATccc                                                                                                                                                 >  2:713643/1‑150 (MQ=255)
       ttatCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCg                                                                                                                                                <  2:3627006/150‑1 (MQ=255)
       ttatCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCg                                                                                                                                                <  2:4440550/150‑1 (MQ=255)
       ttatCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCg                                                                                                                                                >  2:4338963/1‑150 (MQ=255)
       ttatCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCg                                                                                                                                                <  2:1984339/150‑1 (MQ=255)
       ttatCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCg                                                                                                                                                <  2:2038871/150‑1 (MQ=255)
       ttatCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCg                                                                                                                                                >  2:1051433/1‑150 (MQ=255)
        tatCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGc                                                                                                                                               >  2:3424928/1‑150 (MQ=255)
         atCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCt                                                                                                                                              <  2:1619774/150‑1 (MQ=255)
            tGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGgc                                                                                                                                           >  1:1709343/1‑150 (MQ=255)
             gAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGgcg                                                                                                                                          >  1:4244794/1‑150 (MQ=255)
                ccGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCgcag                                                                                                                                       <  2:2916924/150‑1 (MQ=255)
                 cGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCgcagc                                                                                                                                      >  2:2535365/1‑150 (MQ=255)
                 cGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCgcagc                                                                                                                                      >  2:140414/1‑150 (MQ=255)
                  gCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCgcagca                                                                                                                                     >  1:2892433/1‑150 (MQ=255)
                   ccATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGcagcag                                                                                                                                    <  1:2610531/150‑1 (MQ=255)
                    cagttGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGt                                                                                                                                   <  1:2481350/147‑1 (MQ=255)
                     aTTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGtt                                                                                                                                  >  1:3915078/1‑150 (MQ=255)
                     aTTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGtt                                                                                                                                  >  2:602843/1‑150 (MQ=255)
                       ttGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCt                                                                                                                                >  1:4092547/1‑150 (MQ=255)
                       ttGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCt                                                                                                                                <  1:773615/150‑1 (MQ=255)
                       ttGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCt                                                                                                                                >  2:3763798/1‑150 (MQ=255)
                       ttGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCt                                                                                                                                >  2:2270700/1‑150 (MQ=255)
                        tGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTg                                                                                                                               <  1:4338963/150‑1 (MQ=255)
                        tGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTg                                                                                                                               >  2:2996426/1‑150 (MQ=255)
                        tGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTg                                                                                                                               >  2:532882/1‑150 (MQ=255)
                            cTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCaa                                                                                                                           <  1:1053846/150‑1 (MQ=255)
                             tGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCaaa                                                                                                                          >  2:2817406/1‑150 (MQ=255)
                              gACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAAc                                                                                                                         >  1:3486937/1‑150 (MQ=255)
                               aCGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACg                                                                                                                        >  1:1625788/1‑150 (MQ=255)
                                    cGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGc                                                                                                                   >  2:4767977/1‑150 (MQ=255)
                                     gTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGcc                                                                                                                  >  1:1091458/1‑150 (MQ=255)
                                          cGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAgaga                                                                                                             >  1:3610681/1‑150 (MQ=255)
                                           gACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGAt                                                                                                            <  1:3472595/150‑1 (MQ=255)
                                            aCATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATc                                                                                                           <  2:4533433/150‑1 (MQ=255)
                                              aTGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCcg                                                                                                         >  1:3838646/1‑150 (MQ=255)
                                              aTGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGGTCTGGCAAACGGTCGCCAGAGATCcg                                                                                                         >  1:2561399/1‑150 (MQ=255)
                                              aTGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGGAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCcg                                                                                                         <  2:4236337/150‑1 (MQ=255)
                                                gctgcCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCcgcg                                                                                                       >  2:2798221/1‑150 (MQ=255)
                                                     cGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCa                                                                                                  <  1:4786181/150‑1 (MQ=255)
                                                       aTGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGc                                                                                                <  2:4318726/150‑1 (MQ=255)
                                                        tGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCg                                                                                               >  1:2108121/1‑150 (MQ=255)
                                                         ggTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGa                                                                                              >  2:3314432/1‑150 (MQ=255)
                                                               tGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGtt                                                                                        <  1:3562123/150‑1 (MQ=255)
                                                                      aaCAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCcac                                                                                 >  2:1151043/1‑150 (MQ=255)
                                                                        cAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACAt                                                                               >  1:4266822/1‑150 (MQ=255)
                                                                                 aaTGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGa                                                                      <  1:2852580/150‑1 (MQ=255)
                                                                                     ctgcGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATg                                                                  <  2:3817620/150‑1 (MQ=255)
                                                                                        cGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCg                                                               >  2:4293969/1‑150 (MQ=255)
                                                                                         gTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCgg                                                              <  1:1325337/150‑1 (MQ=255)
                                                                                          tGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGa                                                             <  2:1220920/150‑1 (MQ=255)
                                                                                              gctgcCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGaaaaa                                                         >  1:650649/1‑150 (MQ=255)
                                                                                               ctgcCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGaaaaaa                                                        >  2:4590092/1‑150 (MQ=255)
                                                                                               ctgcCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGaaaaaa                                                        <  2:2309504/150‑1 (MQ=255)
                                                                                                tgcCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAAt                                                       >  2:3054826/1‑150 (MQ=255)
                                                                                                tgcCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAAt                                                       >  2:265217/1‑150 (MQ=255)
                                                                                                       gAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCaa                                                <  1:908213/150‑1 (MQ=255)
                                                                                                          gTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATc                                             <  1:1917109/150‑1 (MQ=255)
                                                                                                           tGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCg                                            <  1:4767977/150‑1 (MQ=255)
                                                                                                           tGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCCCATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCg                                            >  2:2669024/1‑150 (MQ=255)
                                                                                                           tGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCg                                            >  1:1435532/1‑150 (MQ=255)
                                                                                                            gAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAGCGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGa                                           <  2:250302/150‑1 (MQ=255)
                                                                                                                tttGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAagag                                       >  2:1793560/1‑150 (MQ=255)
                                                                                                                    aTCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAg                                   >  1:710875/1‑150 (MQ=255)
                                                                                                                       ttAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTg                                <  1:2571428/150‑1 (MQ=255)
                                                                                                                        tAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGc                               >  1:2241623/1‑150 (MQ=255)
                                                                                                                          aCGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGccc                             <  2:3237154/150‑1 (MQ=255)
                                                                                                                           cGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCg                            <  2:4693469/150‑1 (MQ=255)
                                                                                                                            gTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGc                           >  1:2718765/1‑150 (MQ=255)
                                                                                                                                gTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATgg                       <  1:3088682/150‑1 (MQ=255)
                                                                                                                                 tAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATggg                      >  1:2170473/1‑150 (MQ=255)
                                                                                                                                  aaaaGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGc                     >  2:834563/1‑150 (MQ=255)
                                                                                                                                  aaaaGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGc                     >  1:79033/1‑150 (MQ=255)
                                                                                                                                  aaaaGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGc                     <  1:4577677/150‑1 (MQ=255)
                                                                                                                                  aaaaGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGc                     >  1:921141/1‑150 (MQ=255)
                                                                                                                                   aaaGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCg                    >  1:2777817/1‑150 (MQ=255)
                                                                                                                                      gcgcTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGt                 >  2:1918687/1‑150 (MQ=255)
                                                                                                                                        gcTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTaa               <  1:1771691/150‑1 (MQ=255)
                                                                                                                                        gcTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTaa               <  1:3739000/150‑1 (MQ=255)
                                                                                                                                        gcTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTaa               <  1:545611/150‑1 (MQ=255)
                                                                                                                                           tGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTa            <  2:1516329/150‑1 (MQ=255)
                                                                                                                                           tGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTa            >  2:4018673/1‑150 (MQ=255)
                                                                                                                                           tGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTa            >  2:2581744/1‑150 (MQ=255)
                                                                                                                                            ggTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTAc           <  1:2451686/150‑1 (MQ=255)
                                                                                                                                             gTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTAcc          <  1:3943650/150‑1 (MQ=255)
                                                                                                                                             cttAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCAAGGGCGGGTAACTAcc          >  2:2251720/2‑150 (MQ=255)
                                                                                                                                              ttAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCt         >  1:2820215/1‑150 (MQ=255)
                                                                                                                                                aaCGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACCTTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTa       >  1:1531392/1‑150 (MQ=255)
                                                                                                                                                aaCGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTa       >  1:3361070/1‑150 (MQ=255)
                                                                                                                                                 aCGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTAc      >  1:3241188/1‑150 (MQ=255)
                                                                                                                                                  cGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCCCATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACa     >  1:3987355/1‑150 (MQ=255)
                                                                                                                                                  cGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACCTTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACa     >  1:1255898/1‑150 (MQ=255)
                                                                                                                                                  cGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACa     <  1:2177194/150‑1 (MQ=255)
                                                                                                                                                     cGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACAGcc  >  1:3186098/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACAGCC  >  CP009273/4612632‑4612930

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: