Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 796,379 A→T 5.0% P55P (CCA→CCT ybhH → putative PrpF family isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273796,3790AT 5.0% 1013.7 / 28.3 289P55P (CCA→CCTybhHputative PrpF family isomerase
Reads supporting (aligned to +/- strand):  ref base A (130/143);  new base T (7/8);  total (137/152)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 292 total aligned reads displayed.

CGTGATGATGCGAGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGCATTCGCCAT  >  CP009273/796232‑796528
                                                                                                                                                   |                                                                                                                                                     
cGTGATGATGCGAGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACt                                                                                                                                                     >  2:2965619/1‑150 (MQ=255)
 gTGATGATGCGAGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTg                                                                                                                                                    >  2:2604877/1‑150 (MQ=255)
  tgatgaTGCGAGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCTCTGa                                                                                                                                                   <  2:2240460/150‑1 (MQ=255)
  tgatgaTGCGAGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGa                                                                                                                                                   <  2:2106345/150‑1 (MQ=255)
   gatgatGCGAGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGAc                                                                                                                                                  <  1:4598926/150‑1 (MQ=255)
       atgCGAGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGt                                                                                                                                              <  1:415122/150‑1 (MQ=255)
          cGAGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTaaa                                                                                                                                           <  2:4387992/150‑1 (MQ=255)
           gAGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAg                                                                                                                                          <  1:208085/150‑1 (MQ=255)
            aGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGt                                                                                                                                         >  1:2197573/1‑150 (MQ=255)
             ggTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTc                                                                                                                                        <  1:3068437/150‑1 (MQ=255)
                ggAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGACGcc                                                                                                                                     >  2:1148682/1‑150 (MQ=255)
                  aaCCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCat                                                                                                                                   >  1:4787276/1‑150 (MQ=255)
                    ccTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCatta                                                                                                                                 >  2:4859960/1‑150 (MQ=255)
                    ccTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCatta                                                                                                                                 <  2:3237091/150‑1 (MQ=255)
                     cTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCattat                                                                                                                                <  2:1395885/150‑1 (MQ=255)
                      tCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCattatt                                                                                                                               >  2:538926/1‑150 (MQ=255)
                       cGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCAttatta                                                                                                                              <  2:94716/150‑1 (MQ=255)
                         aGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGc                                                                                                                            <  1:2051257/150‑1 (MQ=255)
                          ggggCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGcc                                                                                                                           <  1:2313722/150‑1 (MQ=255)
                          ggggCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGcc                                                                                                                           <  1:3804229/150‑1 (MQ=255)
                           gggCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCg                                                                                                                          <  2:3405601/150‑1 (MQ=255)
                                cgTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCa                                                                                                                     <  2:2317515/150‑1 (MQ=255)
                                cgTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCa                                                                                                                     <  1:3109104/150‑1 (MQ=255)
                                cgTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCa                                                                                                                     <  2:5047266/150‑1 (MQ=255)
                                   tCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCTCTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCg                                                                                                                  >  2:45019/1‑150 (MQ=255)
                                   tCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCg                                                                                                                  <  2:923193/150‑1 (MQ=255)
                                       gTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATcc                                                                                                              >  2:3146633/1‑150 (MQ=255)
                                         tAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCcgc                                                                                                            >  1:3434944/1‑150 (MQ=255)
                                         tAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCcgc                                                                                                            <  1:2856726/150‑1 (MQ=255)
                                          aGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCcgcg                                                                                                           <  1:999869/150‑1 (MQ=255)
                                           gCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGt                                                                                                          <  1:5130051/150‑1 (MQ=255)
                                             ggAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGc                                                                                                        >  2:4595399/1‑150 (MQ=255)
                                                aCATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGa                                                                                                     >  1:3201770/1‑150 (MQ=255)
                                                  aTTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATg                                                                                                   <  2:2849165/150‑1 (MQ=255)
                                                    ttACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGGAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTc                                                                                                 >  1:1069362/1‑150 (MQ=255)
                                                     tACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCg                                                                                                <  2:2834/150‑1 (MQ=255)
                                                     tACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCg                                                                                                <  2:3302569/150‑1 (MQ=255)
                                                     tACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCg                                                                                                <  1:3638130/150‑1 (MQ=255)
                                                       cccGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGAt                                                                                              >  1:4096780/1‑150 (MQ=255)
                                                        ccGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGAtt                                                                                             <  2:3103419/150‑1 (MQ=255)
                                                             gATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTg                                                                                        <  2:4437146/150‑1 (MQ=255)
                                                                cAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGttt                                                                                     <  1:1223160/150‑1 (MQ=255)
                                                                 aaaCGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTg                                                                                    <  2:4321071/150‑1 (MQ=255)
                                                                    cGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTgcgc                                                                                 >  2:406233/1‑150 (MQ=255)
                                                                        gcgcGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGt                                                                             >  1:1418771/1‑150 (MQ=255)
                                                                            gATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATc                                                                         >  2:4556205/1‑150 (MQ=255)
                                                                                aaaTATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTc                                                                     <  1:3597622/150‑1 (MQ=255)
                                                                                aaaTATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTc                                                                     >  2:3828471/1‑150 (MQ=255)
                                                                                aaaTATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTc                                                                     <  2:855200/150‑1 (MQ=255)
                                                                                  atatTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCAt                                                                   >  1:356743/1‑150 (MQ=255)
                                                                                    atTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGa                                                                 >  1:1928627/1‑150 (MQ=255)
                                                                                    atTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGa                                                                 <  2:4900997/150‑1 (MQ=255)
                                                                                     ttGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAg                                                                <  1:3657340/150‑1 (MQ=255)
                                                                                     ttGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAg                                                                >  2:1749683/1‑150 (MQ=255)
                                                                                       gATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCa                                                              <  1:368779/150‑1 (MQ=255)
                                                                                        aTGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCaa                                                             <  2:2427269/150‑1 (MQ=255)
                                                                                          ggCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACg                                                           <  1:1009302/150‑1 (MQ=255)
                                                                                            cAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACgtg                                                         <  1:2022692/150‑1 (MQ=255)
                                                                                             aaTTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACgtgt                                                        >  2:4152042/1‑150 (MQ=255)
                                                                                               ttATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCg                                                      <  1:4297879/150‑1 (MQ=255)
                                                                                                 aTGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGAt                                                    <  2:421385/150‑1 (MQ=255)
                                                                                                 aTGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGAt                                                    <  2:4699201/150‑1 (MQ=255)
                                                                                                   gggTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAc                                                  >  2:3319952/1‑150 (MQ=255)
                                                                                                    ggTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAcc                                                 <  1:1674375/150‑1 (MQ=255)
                                                                                                     gTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCTCTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCa                                                >  1:356125/1‑150 (MQ=255)
                                                                                                      ttCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTATTCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCAc                                               <  2:2873341/150‑1 (MQ=255)
                                                                                                          ggTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCTCTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCt                                           >  1:3232368/1‑150 (MQ=255)
                                                                                                          ggTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCt                                           >  1:1711044/1‑150 (MQ=255)
                                                                                                          ggTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCt                                           <  2:4959812/150‑1 (MQ=255)
                                                                                                           gTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTa                                          >  2:3007861/1‑150 (MQ=255)
                                                                                                           gTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTa                                          <  1:4626636/150‑1 (MQ=255)
                                                                                                           gTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTAACTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTa                                          >  2:3610587/1‑150 (MQ=255)
                                                                                                            tAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTaa                                         <  1:3343919/150‑1 (MQ=255)
                                                                                                             aaCGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCTCTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAAc                                        <  2:2698475/150‑1 (MQ=255)
                                                                                                             aaCGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAAc                                        <  1:26774/150‑1 (MQ=255)
                                                                                                             aaCGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAAc                                        >  1:651120/1‑150 (MQ=255)
                                                                                                             aaCGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAAc                                        >  2:3020090/1‑150 (MQ=255)
                                                                                                                gATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGc                                     <  1:875992/150‑1 (MQ=255)
                                                                                                                  tCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCgg                                   >  1:4813631/1‑150 (MQ=255)
                                                                                                                  tCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCgg                                   >  2:3659943/1‑150 (MQ=255)
                                                                                                                  tCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCgg                                   <  1:1915340/150‑1 (MQ=255)
                                                                                                                   cTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGc                                  <  2:662104/150‑1 (MQ=255)
                                                                                                                   cTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGc                                  >  2:1812087/1‑150 (MQ=255)
                                                                                                                     ggAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCaa                                >  1:4474744/1‑150 (MQ=255)
                                                                                                                      gAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAat                               >  1:824937/1‑150 (MQ=255)
                                                                                                                      gAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAat                               >  2:3658743/1‑150 (MQ=255)
                                                                                                                         aTTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGGTTGCGCAGGTACTCGTTCATGAGCAACGAGTCGATACCACGCCTAACTGCGGCAATATg                            >  2:201711/1‑150 (MQ=255)
                                                                                                                          ttGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGc                           <  1:4610313/150‑1 (MQ=255)
                                                                                                                            gACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTg                         <  2:1910982/150‑1 (MQ=255)
                                                                                                                               ggTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCt                      <  2:3722368/150‑1 (MQ=255)
                                                                                                                               ggTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCt                      <  1:2989400/150‑1 (MQ=255)
                                                                                                                                 tATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTgg                    >  2:3597677/1‑150 (MQ=255)
                                                                                                                                 tATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTgg                    <  1:4173324/150‑1 (MQ=255)
                                                                                                                                      gcggcgGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTg               <  2:2054856/150‑1 (MQ=255)
                                                                                                                                       cggcggTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCCACGTGTCGATACCACGCCTAACTGCGGCAATATGCGGTCTGGTGTTgg              >  2:1674502/1‑150 (MQ=255)
                                                                                                                                       cggcggTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTgg              >  2:2158481/1‑150 (MQ=255)
                                                                                                                                         gcggTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTgggg            <  1:1800825/150‑1 (MQ=255)
                                                                                                                                           ggTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGCa          >  2:3915323/1‑150 (MQ=255)
                                                                                                                                              aaTCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGCATTc       <  1:1351852/150‑1 (MQ=255)
                                                                                                                                                   aCTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGCATTCGCCAt  <  1:4983292/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGTGATGATGCGAGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGCATTCGCCAT  >  CP009273/796232‑796528

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: