Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,067,329 A→G 10.2% intergenic (‑279/+6) epd ← / ← yggC D‑erythrose 4‑phosphate dehydrogenase/putative PanK family P‑loop kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,067,3290AG10.2% 890.6 / 72.9 284intergenic (‑279/+6)epd/yggCD‑erythrose 4‑phosphate dehydrogenase/putative PanK family P‑loop kinase
Reads supporting (aligned to +/- strand):  ref base A (136/119);  new base G (19/10);  total (155/129)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 289 total aligned reads displayed.

CTTCAGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGTTAATCCTCCCGC  >  CP009273/3067181‑3067478
                                                                                                                                                    |                                                                                                                                                     
    aGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTg                                                                                                                                                  <  2:2168258/150‑1 (MQ=255)
    aGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTg                                                                                                                                                  >  1:4217740/1‑150 (MQ=255)
     gCTAGGATAAGCTAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGt                                                                                                                                                 <  1:3751786/150‑1 (MQ=255)
       tAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCa                                                                                                                                               <  2:812583/150‑1 (MQ=255)
         ggATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGc                                                                                                                                             >  2:2756732/1‑150 (MQ=255)
           aTAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCtt                                                                                                                                           <  1:115155/150‑1 (MQ=255)
           aTAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCtt                                                                                                                                           <  2:2785734/150‑1 (MQ=255)
             aaGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTtgt                                                                                                                                         <  2:2550656/150‑1 (MQ=255)
             aaGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTtgt                                                                                                                                         <  2:287367/150‑1 (MQ=255)
                cGAAACGTGGAATAAAAGGAATGTTTGTCCATCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATAACATTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAAAACTACAATAACTGGAATAAATTGTCAGCTTTTGaa                                                                                                                                      >  1:485802/1‑150 (MQ=255)
                cGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATGAATTGTCAGCTTGTGaa                                                                                                                                      >  1:2680086/1‑150 (MQ=255)
                cGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGaa                                                                                                                                      >  2:4444401/1‑150 (MQ=255)
                cGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGaa                                                                                                                                      >  1:2461572/1‑150 (MQ=255)
                  aaaCGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGt                                                                                                                                    >  1:147160/1‑150 (MQ=255)
                  aaaCGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGt                                                                                                                                    >  1:924476/1‑150 (MQ=255)
                   aaCGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTg                                                                                                                                   <  2:3702268/150‑1 (MQ=255)
                    aCGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGa                                                                                                                                  <  1:2654786/150‑1 (MQ=255)
                       tGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATAc                                                                                                                               >  2:1045352/1‑150 (MQ=255)
                        ggAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATAcc                                                                                                                              <  1:2237374/150‑1 (MQ=255)
                        ggAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATAcc                                                                                                                              <  1:236228/150‑1 (MQ=255)
                          aaTAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCgg                                                                                                                            >  2:4319291/1‑150 (MQ=255)
                           aTAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGc                                                                                                                           >  1:2272566/1‑150 (MQ=255)
                             aaaaGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCct                                                                                                                         >  1:636689/1‑150 (MQ=255)
                              aaaGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCctt                                                                                                                        <  1:3413197/150‑1 (MQ=255)
                               attgAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGTATCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCcttc                                                                                                                       <  1:1788585/147‑1 (MQ=255)
                               aaGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCcttc                                                                                                                       >  1:770531/1‑150 (MQ=255)
                               aaGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCcttc                                                                                                                       <  1:517999/150‑1 (MQ=255)
                                 ggAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCcttctt                                                                                                                     <  1:966493/150‑1 (MQ=255)
                                  gAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCttcttc                                                                                                                    <  2:3831104/150‑1 (MQ=255)
                                  gAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCttcttc                                                                                                                    >  2:3882733/1‑150 (MQ=255)
                                     tGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGt                                                                                                                 <  2:3311133/150‑1 (MQ=255)
                                     tGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGt                                                                                                                 <  1:3016304/150‑1 (MQ=255)
                                      gTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTc                                                                                                                >  1:4491823/1‑150 (MQ=255)
                                       tttGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCa                                                                                                               >  1:2962008/1‑150 (MQ=255)
                                          gTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCAtct                                                                                                            >  2:3982502/1‑150 (MQ=255)
                                          gTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCAtct                                                                                                            <  2:456619/150‑1 (MQ=255)
                                            ccAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATctct                                                                                                          <  2:3099556/150‑1 (MQ=255)
                                             cAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTcta                                                                                                         <  2:1475162/150‑1 (MQ=255)
                                              aGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTctac                                                                                                        >  2:3375242/1‑150 (MQ=255)
                                               gCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTctact                                                                                                       <  1:2925807/150‑1 (MQ=255)
                                                ccGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTctacta                                                                                                      >  1:1747560/1‑150 (MQ=255)
                                                 cGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTAt                                                                                                     <  2:1658592/150‑1 (MQ=255)
                                                  gaagaaGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATc                                                                                                    >  2:1065015/1‑150 (MQ=255)
                                                  gaagaaGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATc                                                                                                    >  2:2821462/1‑150 (MQ=255)
                                                  gaagaaGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATc                                                                                                    >  2:1935935/1‑150 (MQ=255)
                                                    agaagaCATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAg                                                                                                  <  2:2814603/150‑1 (MQ=255)
                                                    agaagaCATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAg                                                                                                  >  1:1884388/1‑150 (MQ=255)
                                                      aagaCATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGAt                                                                                                >  2:2449616/1‑150 (MQ=255)
                                                      aagaCATTTATCTGACTCACATCACACTTTTATCCCCTTGTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGAt                                                                                                >  2:3604944/1‑150 (MQ=255)
                                                       agaCATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGAtt                                                                                               <  2:713368/150‑1 (MQ=255)
                                                        gaCATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTg                                                                                              <  2:3448862/150‑1 (MQ=255)
                                                         aCATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTgg                                                                                             <  1:4423099/150‑1 (MQ=255)
                                                         aCATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTgg                                                                                             <  1:1040503/150‑1 (MQ=255)
                                                            tttATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCt                                                                                          <  2:3768264/150‑1 (MQ=255)
                                                             ttATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTg                                                                                         >  2:4171475/1‑150 (MQ=255)
                                                               aTCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGAc                                                                                       >  1:346627/1‑150 (MQ=255)
                                                                   gACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCt                                                                                   <  1:903454/150‑1 (MQ=255)
                                                                    aCTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTa                                                                                  <  1:2068771/150‑1 (MQ=255)
                                                                       cacatcacaCTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATtc                                                                               <  2:4473964/150‑1 (MQ=255)
                                                                       cacatcacaCTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATtc                                                                               >  1:1080844/1‑150 (MQ=255)
                                                                        acatcacaCTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATtca                                                                              <  1:3105459/150‑1 (MQ=255)
                                                                        acatcacaCTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATtca                                                                              <  1:3312773/150‑1 (MQ=255)
                                                                         catcacaCTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATtcat                                                                             <  1:2556555/150‑1 (MQ=255)
                                                                            cacaCTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGATTGGCCTGACGGCTATTCATCAg                                                                          >  2:1681415/1‑150 (MQ=255)
                                                                              cacTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATACCTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGca                                                                        <  1:356959/150‑1 (MQ=255)
                                                                              cacTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGca                                                                        <  1:3894831/150‑1 (MQ=255)
                                                                                     aTCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTc                                                                 >  1:3515338/1‑150 (MQ=255)
                                                                                      tCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCg                                                                >  1:3352739/1‑150 (MQ=255)
                                                                                      tCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCg                                                                <  1:2757625/150‑1 (MQ=255)
                                                                                          cTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACAt                                                            >  2:1996601/1‑150 (MQ=255)
                                                                                             ttGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTgg                                                         <  1:4334826/150‑1 (MQ=255)
                                                                                                  ggAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCa                                                    >  2:4401481/1‑150 (MQ=255)
                                                                                                   gAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCAt                                                   >  2:1123419/1‑150 (MQ=255)
                                                                                                      gCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCgg                                                >  1:1596153/1‑150 (MQ=255)
                                                                                                       cTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGt                                               <  2:3979323/150‑1 (MQ=255)
                                                                                                           aTTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCgg                                           >  2:2358054/1‑150 (MQ=255)
                                                                                                              ccAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCa                                        <  2:285290/150‑1 (MQ=255)
                                                                                                              ccAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCa                                        >  1:1821989/1‑150 (MQ=255)
                                                                                                                aGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATa                                      <  1:2382447/150‑1 (MQ=255)
                                                                                                                  gCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATaaa                                    >  2:738065/1‑150 (MQ=255)
                                                                                                                    tGGCGTAATAATAACCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAAt                                  <  2:665352/150‑1 (MQ=255)
                                                                                                                       cGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCt                               <  1:280661/150‑1 (MQ=255)
                                                                                                                          aataataaCCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCg                            >  2:243280/1‑150 (MQ=255)
                                                                                                                          aataataaCCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCg                            <  1:4456354/150‑1 (MQ=255)
                                                                                                                           ataataaCCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCgg                           >  2:3911489/1‑150 (MQ=255)
                                                                                                                             aataaCCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCt                         <  1:1916993/150‑1 (MQ=255)
                                                                                                                              ataaCCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTa                        >  2:2036377/1‑150 (MQ=255)
                                                                                                                               taaCCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTAc                       >  2:539669/1‑150 (MQ=255)
                                                                                                                                aaCCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACt                      >  1:4278431/1‑150 (MQ=255)
                                                                                                                                     tACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCcg                 >  1:3677098/1‑150 (MQ=255)
                                                                                                                                      aCAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCcgc                >  1:1182669/1‑150 (MQ=255)
                                                                                                                                      aCAATAACTGGAATAAATTGTCAACTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCcgc                >  1:3116356/1‑150 (MQ=255)
                                                                                                                                         aTAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGtt             <  1:3081816/150‑1 (MQ=255)
                                                                                                                                          tAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGTTa            >  1:4123109/1‑150 (MQ=255)
                                                                                                                                          tAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGTTa            >  2:1763972/1‑150 (MQ=255)
                                                                                                                                           aaCTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGTTaa           >  2:1011357/1‑150 (MQ=255)
                                                                                                                                                 aaTAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGTTAAtcctcc     >  2:372244/1‑150 (MQ=255)
                                                                                                                                                 aaTAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGTTAAtcctcc     >  1:4053341/1‑150 (MQ=255)
                                                                                                                                                  aTAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGTTAATCCTccc    <  1:3928187/150‑1 (MQ=255)
                                                                                                                                                   tAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGTTAATCCTCCCg   <  1:2648178/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTTCAGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGTTAATCCTCCCGC  >  CP009273/3067181‑3067478

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: