Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,768,689 C→A 8.7% S215* (TCA→TAA)  ydiB → quinate/shikimate 5‑dehydrogenase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,768,6890CA 8.7% 1024.1 / 45.6 274S215* (TCA→TAA) ydiBquinate/shikimate 5‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (115/135);  new base A (16/8);  total (131/143)
Fisher's exact test for biased strand distribution p-value = 5.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 276 total aligned reads displayed.

CGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGGCTG  >  CP009273/1768540‑1768838
                                                                                                                                                     |                                                                                                                                                     
    ttCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTg                                                                                                                                                   <  2:663207/150‑1 (MQ=255)
      cgcgCATCTGGTTTATTAAAACACCGTTTGTTTCGTCACGGTCACCGATATCGCCGTTCAGCAAGCCTTTGCTGAAGACATTGCTTCCTCCGACATTTTAACCAATGGCAAAAAAGTGGGTATGTAACCCCTTGAGAATGAATCATTTGt                                                                                                                                                 <  2:240918/150‑1 (MQ=255)
      cgcgCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGt                                                                                                                                                 <  1:1948379/150‑1 (MQ=255)
        cgcAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTa                                                                                                                                               >  2:3265183/1‑150 (MQ=255)
        cgcAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTa                                                                                                                                               <  2:3874084/150‑1 (MQ=255)
             cGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGat                                                                                                                                          <  1:4175787/150‑1 (MQ=255)
             cGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGat                                                                                                                                          <  1:4175775/150‑1 (MQ=255)
                gTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATc                                                                                                                                       >  1:1443459/1‑150 (MQ=255)
                  tAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAg                                                                                                                                     >  1:2455884/1‑150 (MQ=255)
                   aaTGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGt                                                                                                                                    >  2:2834679/1‑150 (MQ=255)
                   aaTGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGt                                                                                                                                    >  2:1511997/1‑150 (MQ=255)
                     tGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCt                                                                                                                                  <  1:3252738/150‑1 (MQ=255)
                     tGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCt                                                                                                                                  >  1:579228/1‑150 (MQ=255)
                      gAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTg                                                                                                                                 >  1:2722911/1‑150 (MQ=255)
                       aaaaCACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGt                                                                                                                                >  1:9768/1‑150 (MQ=255)
                       aaaaCACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGt                                                                                                                                <  2:2152227/150‑1 (MQ=255)
                        aaaCACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAACAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGtt                                                                                                                               <  2:3641247/150‑1 (MQ=255)
                         aaCACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTa                                                                                                                              <  2:2983197/150‑1 (MQ=255)
                            acCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACAt                                                                                                                           >  1:2519844/1‑150 (MQ=255)
                            acCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACAt                                                                                                                           <  2:1170092/150‑1 (MQ=255)
                            acCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACAt                                                                                                                           >  2:219793/1‑150 (MQ=255)
                              cGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATcc                                                                                                                         >  2:4620736/1‑150 (MQ=255)
                              cGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATcc                                                                                                                         >  1:1601654/1‑150 (MQ=255)
                                 ttGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCggg                                                                                                                      <  1:1701911/150‑1 (MQ=255)
                                   gtgtCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCAAAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGAc                                                                                                                    <  1:1022639/150‑1 (MQ=255)
                                       cgtcACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCt                                                                                                                >  1:4406413/1‑150 (MQ=255)
                                           aCGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTc                                                                                                            <  2:39284/150‑1 (MQ=255)
                                           aCGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTc                                                                                                            >  2:1382182/1‑150 (MQ=255)
                                           aCGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTc                                                                                                            >  1:852744/1‑150 (MQ=255)
                                            cGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCa                                                                                                           <  1:2449347/150‑1 (MQ=255)
                                             ggTCACCGATCTCGCCGATCAGCGAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCAc                                                                                                          >  2:2634675/1‑150 (MQ=255)
                                             ggTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCAc                                                                                                          <  2:4221525/150‑1 (MQ=255)
                                              gTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACt                                                                                                         <  1:2182929/150‑1 (MQ=255)
                                                  ccGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAAt                                                                                                     >  1:259724/1‑150 (MQ=255)
                                                   cGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATg                                                                                                    <  1:3182488/150‑1 (MQ=255)
                                                   cGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATg                                                                                                    >  2:557052/1‑150 (MQ=255)
                                                   cGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATg                                                                                                    <  1:3182513/150‑1 (MQ=255)
                                                       ctcGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCgtg                                                                                                >  2:2240776/1‑150 (MQ=255)
                                                            cGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATaa                                                                                           <  2:127192/150‑1 (MQ=255)
                                                            cGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATaa                                                                                           <  1:1881524/150‑1 (MQ=255)
                                                            cGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATaa                                                                                           <  2:519990/150‑1 (MQ=255)
                                                            cGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATaa                                                                                           <  2:1209146/150‑1 (MQ=255)
                                                             gATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAAc                                                                                          <  2:2062509/150‑1 (MQ=255)
                                                             gATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAAc                                                                                          >  1:1358159/1‑150 (MQ=255)
                                                             gATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAAc                                                                                          >  1:1357556/1‑150 (MQ=255)
                                                               tCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAAccc                                                                                        <  1:3784545/150‑1 (MQ=255)
                                                                   cAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCat                                                                                    >  1:873527/1‑150 (MQ=255)
                                                                       ccTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGa                                                                                >  2:3105103/1‑150 (MQ=255)
                                                                        cTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGAc                                                                               <  1:360847/150‑1 (MQ=255)
                                                                        cTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGAc                                                                               <  1:2932109/150‑1 (MQ=255)
                                                                                 aGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTAtt                                                                      <  2:4573804/150‑1 (MQ=255)
                                                                                   cccTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgc                                                                    <  1:4083200/150‑1 (MQ=255)
                                                                                    ccTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgca                                                                   <  1:4079013/150‑1 (MQ=255)
                                                                                    ccTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgca                                                                   >  2:3236390/1‑150 (MQ=255)
                                                                                     cTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgcag                                                                  >  1:3936100/1‑150 (MQ=255)
                                                                                      tGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgcagc                                                                 <  1:553953/150‑1 (MQ=255)
                                                                                       ggCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgcagca                                                                >  2:3543563/1‑150 (MQ=255)
                                                                                       ggCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgcagca                                                                <  2:896581/150‑1 (MQ=255)
                                                                                        gCTTCCGCCGACATTTTAACCAATGGCACCAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGcagcag                                                               <  1:4045047/150‑1 (MQ=255)
                                                                                        gCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGcagcag                                                               >  1:3688887/1‑150 (MQ=255)
                                                                                         cTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAgg                                                              >  1:822350/1‑150 (MQ=255)
                                                                                          ttCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGgc                                                             >  2:821678/1‑150 (MQ=255)
                                                                                           tCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGgcg                                                            <  2:3732089/150‑1 (MQ=255)
                                                                                             cgccgACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGca                                                          <  1:4684/150‑1 (MQ=255)
                                                                                              gccgACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGcaa                                                         >  2:1653284/1‑150 (MQ=255)
                                                                                                cgACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGcaaca                                                       >  1:302464/1‑150 (MQ=255)
                                                                                                cgACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGcaaca                                                       >  1:2615988/1‑150 (MQ=255)
                                                                                                  aCATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAg                                                     >  2:2268362/1‑150 (MQ=255)
                                                                                                    aTTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCt                                                   >  2:2401247/1‑150 (MQ=255)
                                                                                                     ttttAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTg                                                  >  1:109469/1‑150 (MQ=255)
                                                                                                        tAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGtt                                               <  1:3998720/150‑1 (MQ=255)
                                                                                                          aCCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGc                                             >  2:1783167/1‑150 (MQ=255)
                                                                                                           ccAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCa                                            <  2:242849/150‑1 (MQ=255)
                                                                                                           ccAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCa                                            <  2:1601654/150‑1 (MQ=255)
                                                                                                           ccAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCa                                            <  2:2873430/150‑1 (MQ=255)
                                                                                                             aaTGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCaaa                                          >  1:656645/1‑150 (MQ=255)
                                                                                                              aTGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCaaaa                                         >  1:4193619/1‑150 (MQ=255)
                                                                                                                  cacaAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGAt                                     >  2:1821575/1‑150 (MQ=255)
                                                                                                                  cacaAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGAt                                     >  2:3059896/1‑150 (MQ=255)
                                                                                                                    caAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTg                                   <  1:1047053/150‑1 (MQ=255)
                                                                                                                    caAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTg                                   >  2:3363539/1‑150 (MQ=255)
                                                                                                                        aGTGGGTATGAAACCCCTTGAGATTGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATgg                               <  1:2387127/150‑1 (MQ=255)
                                                                                                                           gggTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATa                            <  1:3656211/150‑1 (MQ=255)
                                                                                                                           gggTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATa                            >  2:672660/1‑150 (MQ=255)
                                                                                                                              tATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACgg                         <  2:1915807/150‑1 (MQ=255)
                                                                                                                               aTGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGc                        >  1:4289626/1‑150 (MQ=255)
                                                                                                                                 gAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAt                      <  2:4136976/150‑1 (MQ=255)
                                                                                                                                  aaaCCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtg                     >  1:2072096/1‑150 (MQ=255)
                                                                                                                                    aCCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgtt                   >  2:4075954/1‑150 (MQ=255)
                                                                                                                                     ccccTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttg                  >  1:1215684/1‑150 (MQ=255)
                                                                                                                                       ccTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttgtt                <  1:2211613/150‑1 (MQ=255)
                                                                                                                                        cTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttgttg               <  1:4574598/150‑1 (MQ=255)
                                                                                                                                         ttGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAAAAAGCTGGTTGCAAAAAGATTTATGTATACGGCAtgttgttgt              >  1:1828796/1‑150 (MQ=255)
                                                                                                                                           gagaATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTgg            >  2:40730/1‑150 (MQ=255)
                                                                                                                                           gagaATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTgg            >  2:1616328/1‑150 (MQ=255)
                                                                                                                                             gaatgaatCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCa          <  1:3093168/150‑1 (MQ=255)
                                                                                                                                              aatgaatCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGAATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCaa         >  1:4479447/1‑150 (MQ=255)
                                                                                                                                                 gaatAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCATTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAggg      <  2:807454/150‑1 (MQ=255)
                                                                                                                                                   atCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGGc    <  2:15552/150‑1 (MQ=255)
                                                                                                                                                    tCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGGCt   <  2:259724/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGGCTG  >  CP009273/1768540‑1768838

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: