Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 97,521 C→A 82.0% 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,5210CA82.0% 704.3 / 176.8 315A90A (GCC→GCAmurCUDP‑N‑acetylmuramate:L‑alanine ligase
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base A (147/111);  total (174/141)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 315 total aligned reads displayed.

ATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCA  >  CP009273/97373‑97666
                                                                                                                                                    |                                                                                                                                                 
aTGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCg                                                                                                                                                  >  1:2360703/1‑150 (MQ=255)
  gAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGAt                                                                                                                                                <  1:2277170/150‑1 (MQ=255)
  gAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGAt                                                                                                                                                >  2:3266885/1‑150 (MQ=255)
     ggTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAAc                                                                                                                                             >  2:615464/1‑150 (MQ=255)
        tATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCg                                                                                                                                          >  1:630174/1‑150 (MQ=255)
         atcagatcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCgg                                                                                                                                         <  2:3046363/150‑1 (MQ=255)
         atcagatcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCATATAACCCgg                                                                                                                                         >  1:1906704/1‑150 (MQ=255)
         atcagatcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCgg                                                                                                                                         >  2:3472444/1‑150 (MQ=255)
           cagatcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGaa                                                                                                                                       >  2:593132/1‑150 (MQ=255)
           cagatcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGaa                                                                                                                                       >  2:4054835/1‑150 (MQ=255)
             gatcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAAt                                                                                                                                     >  2:2408432/1‑150 (MQ=255)
              atcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTGGTTTCCAGCGCGATTTCTTCAGATAACCCCGAAAtt                                                                                                                                    >  1:776308/1‑150 (MQ=255)
                  gTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTGGTTTCCAGCGCGATTTCTGCAGATAACCCCGAAAATTTcg                                                                                                                                >  1:792822/1‑150 (MQ=255)
                      ttCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTcgccgc                                                                                                                            >  1:4085055/1‑150 (MQ=255)
                         cGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCa                                                                                                                         >  2:2565250/1‑150 (MQ=255)
                         cGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAAACCGGAAATTGTCGCCGCTCa                                                                                                                         >  1:1124848/1‑150 (MQ=255)
                          gATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCAt                                                                                                                        >  2:1007833/1‑150 (MQ=255)
                           aTTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATg                                                                                                                       >  2:1671971/1‑150 (MQ=255)
                               aGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAgc                                                                                                                   >  1:3223371/1‑150 (MQ=255)
                               aGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAgc                                                                                                                   >  2:1416240/1‑150 (MQ=255)
                                gcgcCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAgcg                                                                                                                  <  2:3600008/150‑1 (MQ=255)
                                  gcCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGcgcg                                                                                                                <  2:3179398/150‑1 (MQ=255)
                                   ccAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGt                                                                                                               >  1:1012152/1‑150 (MQ=255)
                                     aaaTCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTAt                                                                                                             <  1:3338068/150‑1 (MQ=255)
                                        tCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTcc                                                                                                          <  1:1308263/150‑1 (MQ=255)
                                            gTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTg                                                                                                      >  2:2198989/1‑150 (MQ=255)
                                              cACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGAt                                                                                                    <  1:3249598/150‑1 (MQ=255)
                                                cGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATcc                                                                                                  <  1:2733753/150‑1 (MQ=255)
                                                 gcagcaGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcg                                                                                                 <  1:1065666/150‑1 (MQ=255)
                                                   agcagTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgtc                                                                                               >  1:2699728/1‑150 (MQ=255)
                                                   agcagTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTAGCCGCTCATGAAGCGCGTATTCCGGTGATCcgtc                                                                                               >  1:2335434/1‑150 (MQ=255)
                                                    gcagTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgtcg                                                                                              <  1:1685578/150‑1 (MQ=255)
                                                     cagTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgtcgt                                                                                             >  1:2119849/1‑150 (MQ=255)
                                                     cagTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgtcgt                                                                                             <  2:4492372/150‑1 (MQ=255)
                                                      agTTTATGAATCTTTGTGCGACGATTTATTTCAACGTTCGCACGGTATAAGTACGTTATTGCAGCGTGGTCGTTGTTTCCAGAGGTTTTTCTGCAGTTAACCATGATTTTGTAGCCTCTCATGAAGCTCTTATTCAGGTTTTCCGTAGTg                                                                                            <  2:655476/150‑1 (MQ=255)
                                                      agTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTg                                                                                            >  2:1664563/1‑150 (MQ=255)
                                                        ttAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGcc                                                                                          >  2:1857539/1‑150 (MQ=255)
                                                        ttAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGcc                                                                                          >  2:4132151/1‑150 (MQ=255)
                                                          aaTGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGa                                                                                        >  2:1468769/1‑150 (MQ=255)
                                                          aaTGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGa                                                                                        <  1:4408455/150‑1 (MQ=255)
                                                             gAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAAt                                                                                     <  2:2896493/150‑1 (MQ=255)
                                                             gAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAAt                                                                                     <  1:1231949/150‑1 (MQ=255)
                                                              aaTCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATg                                                                                    <  2:1140165/150‑1 (MQ=255)
                                                               aTCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgc                                                                                   <  2:2516402/150‑1 (MQ=255)
                                                                tCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgct                                                                                  <  1:3066699/150‑1 (MQ=255)
                                                                tCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgct                                                                                  <  2:2337333/150‑1 (MQ=255)
                                                                  tGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgctgg                                                                                <  1:4471404/150‑1 (MQ=255)
                                                                  tGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgctgg                                                                                >  1:3956499/1‑150 (MQ=255)
                                                                     gTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgctggctg                                                                             >  1:2039510/1‑150 (MQ=255)
                                                                       gCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTAATCAGATAAACCCGAAAATTTAGCCGCGCAAGAAGCGCGTATTACCGTGATCCGTCGTGCCGCAATTCTGGCGGAg                                                                           >  1:824345/1‑150 (MQ=255)
                                                                         gacgaTTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt                                                                         >  2:268068/1‑150 (MQ=255)
                                                                            gaTTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAAt                                                                      <  2:4002201/150‑1 (MQ=255)
                                                                             atttatttCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATg                                                                     <  1:1118130/150‑1 (MQ=255)
                                                                             atttatttCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATg                                                                     <  1:425334/150‑1 (MQ=255)
                                                                             atttatttCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATg                                                                     >  2:2206615/1‑150 (MQ=255)
                                                                              tttatttCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGc                                                                    <  2:4451797/150‑1 (MQ=255)
                                                                                tatttCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGTAAATTTCGCCCCGAATGGAGCGCGTATTACGGTTGTACGTCGTTCCGAAATGCGGGGGGGGTTAATGCgg                                                                  >  2:4340758/1‑149 (MQ=255)
                                                                                   ttCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGtttt                                                               <  2:3900878/150‑1 (MQ=255)
                                                                                   ttCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGtttt                                                               <  1:2976852/150‑1 (MQ=255)
                                                                                     cAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCg                                                             >  1:3651759/1‑150 (MQ=255)
                                                                                     cAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCg                                                             >  2:4497806/1‑150 (MQ=255)
                                                                                     cAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCg                                                             >  2:3745266/1‑150 (MQ=255)
                                                                                     cAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCg                                                             <  1:9520/150‑1 (MQ=255)
                                                                                      aaCCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAAAGCGTTTTCGt                                                            >  1:1439093/1‑150 (MQ=255)
                                                                                       aCCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTc                                                           >  2:3215872/1‑150 (MQ=255)
                                                                                       aCCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGGTAATGCGTGTTCGTc                                                           >  2:3215180/1‑150 (MQ=255)
                                                                                       aCCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTc                                                           >  2:1183444/1‑150 (MQ=255)
                                                                                         cATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCAt                                                         >  2:2029693/1‑150 (MQ=255)
                                                                                         cATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCAt                                                         >  1:2241220/1‑150 (MQ=255)
                                                                                         cATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCAt                                                         >  1:1747529/1‑150 (MQ=255)
                                                                                          aTCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATg                                                        <  2:896204/150‑1 (MQ=255)
                                                                                          aTCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATg                                                        >  1:4531365/1‑150 (MQ=255)
                                                                                             gCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAAACCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCa                                                     >  1:2573530/1‑150 (MQ=255)
                                                                                               ccGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATc                                                   <  1:4054821/150‑1 (MQ=255)
                                                                                                 ggAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGc                                                 <  1:85672/150‑1 (MQ=255)
                                                                                                   aaaaCGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCa                                               <  2:580282/150‑1 (MQ=255)
                                                                                                   aaaaCGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCa                                               <  2:2454015/150‑1 (MQ=255)
                                                                                                   aaaaCGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTACGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCa                                               >  1:3914010/1‑150 (MQ=255)
                                                                                                    aaaCGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAATCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCGGAGTTAATGCGTTTTCGTAATGGCATAGccct                                              >  1:2818201/1‑148 (MQ=255)
                                                                                                    aaaCGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCAt                                              <  1:3590356/150‑1 (MQ=255)
                                                                                                       cgtacgtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGc                                           <  1:1921728/150‑1 (MQ=255)
                                                                                                        gtacgtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGcc                                          >  2:3725150/1‑150 (MQ=255)
                                                                                                        gtacgtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGcc                                          <  2:262561/150‑1 (MQ=255)
                                                                                                        gtacgtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGcc                                          <  2:2525956/150‑1 (MQ=255)
                                                                                                         tacgtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCg                                         <  1:1497919/150‑1 (MQ=255)
                                                                                                         tacgtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCg                                         >  2:271859/1‑150 (MQ=255)
                                                                                                             tGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAAc                                     <  1:3711065/150‑1 (MQ=255)
                                                                                                             tGATGCCAGCGTGGTCGTTGTGTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAAc                                     <  1:3813867/150‑1 (MQ=255)
                                                                                                                   cAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACgg                               >  2:2155890/1‑150 (MQ=255)
                                                                                                                     gCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCa                             <  1:2498229/150‑1 (MQ=255)
                                                                                                                     gCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCa                             <  2:1300718/150‑1 (MQ=255)
                                                                                                                      cGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCaa                            <  1:725959/150‑1 (MQ=255)
                                                                                                                               tgttTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAac                   >  2:3808432/1‑150 (MQ=255)
                                                                                                                                 tttCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAacca                 <  2:1166244/150‑1 (MQ=255)
                                                                                                                                   tCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAaccacc               >  2:4393839/1‑150 (MQ=255)
                                                                                                                                      aGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACcgcg            >  2:2915024/1‑150 (MQ=255)
                                                                                                                                      aGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACcgcg            <  2:792850/150‑1 (MQ=255)
                                                                                                                                         gcgATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATg         <  1:593132/150‑1 (MQ=255)
                                                                                                                                           gATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGt       >  1:13221/1‑150 (MQ=255)
                                                                                                                                             tttCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGttt     <  1:1947372/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCA  >  CP009273/97373‑97666

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: