Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 97,521 C→A 93.8% A90A (GCC→GCA murC → UDP‑N‑acetylmuramate:L‑alanine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927397,5210CA93.8% 1148.4 / 58.0 365A90A (GCC→GCAmurCUDP‑N‑acetylmuramate:L‑alanine ligase
Reads supporting (aligned to +/- strand):  ref base C (9/14);  new base A (170/172);  total (179/186)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 365 total aligned reads displayed.

GAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAG  >  CP009273/97375‑97667
                                                                                                                                                  |                                                                                                                                                  
 aaGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATa                                                                                                                                                >  1:3857527/1‑150 (MQ=255)
      tATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCg                                                                                                                                           <  1:22111/150‑1 (MQ=255)
         cagatcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGaa                                                                                                                                        >  2:3089366/1‑150 (MQ=255)
         cagatcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGaa                                                                                                                                        >  2:4008736/1‑150 (MQ=255)
          agatcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGaca                                                                                                                                       >  2:3236880/1‑148 (MQ=255)
           gatcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAAt                                                                                                                                      >  2:4752648/1‑150 (MQ=255)
           gatcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAAt                                                                                                                                      <  1:4423929/150‑1 (MQ=255)
              cagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGACATTGt                                                                                                                                   >  1:522405/1‑150 (MQ=255)
                  ggTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTcgcc                                                                                                                               <  1:3463232/150‑1 (MQ=255)
                  ggTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTcgcc                                                                                                                               <  2:3140685/150‑1 (MQ=255)
                  ggTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTcgcc                                                                                                                               >  1:3743166/1‑150 (MQ=255)
                    ttCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTcgccgc                                                                                                                             >  2:2930555/1‑150 (MQ=255)
                      ccGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTc                                                                                                                           >  2:3451778/1‑150 (MQ=255)
                        gATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCAt                                                                                                                         >  2:3052659/1‑150 (MQ=255)
                           ttAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGaa                                                                                                                      <  1:5036445/150‑1 (MQ=255)
                               cgcCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAgcgc                                                                                                                  <  2:4588914/150‑1 (MQ=255)
                                gcCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGcgcg                                                                                                                 >  1:4327650/1‑150 (MQ=255)
                                    aaTCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTAtt                                                                                                             <  2:421244/150‑1 (MQ=255)
                                        cGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCgg                                                                                                         <  1:964884/150‑1 (MQ=255)
                                         ggTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGt                                                                                                        <  1:4618772/150‑1 (MQ=255)
                                          gTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTg                                                                                                       <  1:4074080/150‑1 (MQ=255)
                                           tCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGATAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAATTTGTCGCCGCTCATGAAGCGCGTATTCCGGTGa                                                                                                      <  1:1739061/150‑1 (MQ=255)
                                             aCGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATc                                                                                                    >  1:2657615/1‑150 (MQ=255)
                                             aCGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATc                                                                                                    <  2:765235/150‑1 (MQ=255)
                                              cGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATcc                                                                                                   <  1:835845/150‑1 (MQ=255)
                                              cGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATcc                                                                                                   <  2:494415/150‑1 (MQ=255)
                                               gcagcaGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcg                                                                                                  >  2:689802/1‑150 (MQ=255)
                                                cagcagTTAATGAATCTGGGTGCGTAGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgt                                                                                                 <  1:760750/150‑1 (MQ=255)
                                                cagcagTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgt                                                                                                 <  1:655034/150‑1 (MQ=255)
                                                cagcagTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgt                                                                                                 >  2:137822/1‑150 (MQ=255)
                                                  gcagTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgtcg                                                                                               >  1:1060189/1‑150 (MQ=255)
                                                  gcagTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgtcg                                                                                               >  1:1438635/1‑150 (MQ=255)
                                                   cagTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgtcgt                                                                                              >  2:2906417/1‑150 (MQ=255)
                                                   cagTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgtcgt                                                                                              >  1:4567604/1‑150 (MQ=255)
                                                    agTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTg                                                                                             >  1:3611805/1‑150 (MQ=255)
                                                       tAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCg                                                                                          <  2:2216214/150‑1 (MQ=255)
                                                            aaTCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATg                                                                                     <  2:1811600/150‑1 (MQ=255)
                                                             aTCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgc                                                                                    >  2:199757/1‑150 (MQ=255)
                                                              tCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgct                                                                                   <  1:2459736/150‑1 (MQ=255)
                                                                   gTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgctggctg                                                                              >  1:291096/1‑150 (MQ=255)
                                                                       gacgaTTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt                                                                          >  1:3214973/1‑150 (MQ=255)
                                                                          gaTTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAAt                                                                       <  1:4001784/150‑1 (MQ=255)
                                                                          gaTTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTAGTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGGTAAt                                                                       >  2:906744/1‑150 (MQ=255)
                                                                            tttatttCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGc                                                                     <  2:3721532/150‑1 (MQ=255)
                                                                              tatttCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGt                                                                   <  2:2988826/150‑1 (MQ=255)
                                                                               atttCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGtt                                                                  <  1:3391930/150‑1 (MQ=255)
                                                                                  tCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTc                                                               <  2:3902192/150‑1 (MQ=255)
                                                                                    aaCCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCGGAGTTAATGCGTTTTCGt                                                             >  1:1778089/1‑150 (MQ=255)
                                                                                      ccATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCa                                                           >  1:4361567/1‑150 (MQ=255)
                                                                                      ccATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCa                                                           <  1:2505037/150‑1 (MQ=255)
                                                                                      ccATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGAAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGAATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAAAGCGTTTTCGTCa                                                           >  1:1805605/1‑150 (MQ=255)
                                                                                       cATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCAt                                                          <  2:1572970/150‑1 (MQ=255)
                                                                                        aTCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATg                                                         <  1:789634/150‑1 (MQ=255)
                                                                                        aTCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATg                                                         <  2:2318577/150‑1 (MQ=255)
                                                                                         tCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATgg                                                        >  1:4847310/1‑150 (MQ=255)
                                                                                         tCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATgg                                                        <  1:4918466/150‑1 (MQ=255)
                                                                                          cGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGc                                                       <  1:710132/150‑1 (MQ=255)
                                                                                           gCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCa                                                      >  2:4803848/1‑150 (MQ=255)
                                                                                             ccGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATc                                                    >  1:4907717/1‑150 (MQ=255)
                                                                                             ccGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATc                                                    >  1:1193142/1‑150 (MQ=255)
                                                                                              cGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCg                                                   >  2:1965975/1‑150 (MQ=255)
                                                                                                   aaCGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCAtt                                              >  2:2772455/1‑150 (MQ=255)
                                                                                                   aaCGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCAtt                                              <  1:1887288/150‑1 (MQ=255)
                                                                                                     cgtacgtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGc                                            >  2:1266125/1‑150 (MQ=255)
                                                                                                      gtacgtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGcc                                           <  1:1407225/150‑1 (MQ=255)
                                                                                                      gtacgtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGcc                                           <  1:4158879/150‑1 (MQ=255)
                                                                                                      gtacgtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGGTAATGCGTTTTCGTCATGGCATCGCCATTGcc                                           >  2:4323965/1‑150 (MQ=255)
                                                                                                       tacgtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCg                                          <  2:951502/150‑1 (MQ=255)
                                                                                                        acgtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCgg                                         >  1:2713446/1‑150 (MQ=255)
                                                                                                         cgtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGa                                        >  2:2577440/1‑150 (MQ=255)
                                                                                                           tGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAAc                                      <  1:473083/150‑1 (MQ=255)
                                                                                                           tGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAAc                                      >  2:4092027/1‑150 (MQ=255)
                                                                                                           tGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAAc                                      <  2:3542719/150‑1 (MQ=255)
                                                                                                           tGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAAc                                      <  2:2174410/150‑1 (MQ=255)
                                                                                                           tGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAAc                                      <  1:198446/150‑1 (MQ=255)
                                                                                                            gATGCCAGCGTGGTCGTTGTTTCCAGCGCGAGTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACg                                     <  2:2653966/150‑1 (MQ=255)
                                                                                                              tGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCa                                   <  2:31753/150‑1 (MQ=255)
                                                                                                                ccAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACg                                 <  2:3345692/150‑1 (MQ=255)
                                                                                                                  aGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGc                               >  2:5133273/1‑150 (MQ=255)
                                                                                                                  aGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGc                               >  1:1000319/1‑150 (MQ=255)
                                                                                                                   gCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCa                              <  2:3782052/150‑1 (MQ=255)
                                                                                                                     gTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCaaa                            >  1:989372/1‑150 (MQ=255)
                                                                                                                          cGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGAc                       <  2:673538/150‑1 (MQ=255)
                                                                                                                          cGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGAc                       >  1:3114587/1‑150 (MQ=255)
                                                                                                                          cGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGAc                       <  1:3052659/150‑1 (MQ=255)
                                                                                                                          cGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGAc                       >  2:4391654/1‑150 (MQ=255)
                                                                                                                           gttgttTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACa                      <  2:3857527/150‑1 (MQ=255)
                                                                                                                                ttCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAaccac                 <  2:4708762/150‑1 (MQ=255)
                                                                                                                                 tCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGAATCGCCATTGCCGGAACGCACGGCAAAACGACAaccaca                >  1:1540483/1‑149 (MQ=255)
                                                                                                                                  ccAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACcg               <  1:2386922/150‑1 (MQ=255)
                                                                                                                                  ccAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACcg               <  1:4955219/150‑1 (MQ=255)
                                                                                                                                   cAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACcgc              >  2:1200096/1‑150 (MQ=255)
                                                                                                                                     gcgcGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTACTGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGa            <  2:108211/150‑1 (MQ=255)
                                                                                                                                        cgATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATgg         >  2:3004444/1‑150 (MQ=255)
                                                                                                                                         gATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGt        >  2:148735/1‑150 (MQ=255)
                                                                                                                                          aTTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGtt       >  1:1279776/1‑150 (MQ=255)
                                                                                                                                           tttCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGttt      <  1:2787003/150‑1 (MQ=255)
                                                                                                                                           tttCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGttt      <  1:689847/150‑1 (MQ=255)
                                                                                                                                             tCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTcc    >  1:3153544/1‑150 (MQ=255)
                                                                                                                                              cTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCa   <  1:1581359/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAG  >  CP009273/97375‑97667

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: