Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 97,521 C→A 100% 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,5210CA99.3% 911.3 / ‑2.0 269A90A (GCC→GCAmurCUDP‑N‑acetylmuramate:L‑alanine ligase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (124/143);  total (126/143)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 269 total aligned reads displayed.

AAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCC  >  CP009273/97376‑97665
                                                                                                                                                 |                                                                                                                                                
aaGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGAt                                                                                                                                               <  1:305494/149‑1 (MQ=255)
 aGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATa                                                                                                                                              <  1:2834614/149‑1 (MQ=255)
  ggTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATaa                                                                                                                                             <  2:3788572/149‑1 (MQ=255)
  ggTTATCAGATCAGTGGTTCCGATTGAGCGCCAAATCCGGTCACGCAGCAGTTACTGACTCTGGGTGCGACGATTTATTTCCACCATCGCCCGGAAACCGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATaa                                                                                                                                             <  2:1969202/149‑1 (MQ=255)
   gTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTACAACGCCATTTCTGCAGATAAc                                                                                                                                            >  1:1631228/1‑149 (MQ=255)
         agatcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGaa                                                                                                                                      >  1:2469235/1‑149 (MQ=255)
         agatcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGaa                                                                                                                                      >  2:3310615/1‑149 (MQ=255)
          gatcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAACACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGaaa                                                                                                                                     >  2:3329124/1‑149 (MQ=255)
          gatcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGaaa                                                                                                                                     >  2:1785902/1‑149 (MQ=255)
             cagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTg                                                                                                                                  <  1:3030780/149‑1 (MQ=255)
               gTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGCTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTc                                                                                                                                <  1:2750874/149‑1 (MQ=255)
                tGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTcg                                                                                                                               >  2:2345805/1‑149 (MQ=255)
                 ggTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTcgc                                                                                                                              <  1:2860885/149‑1 (MQ=255)
                 ggTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTcgc                                                                                                                              <  2:4723/149‑1 (MQ=255)
                 ggTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTcgc                                                                                                                              <  1:3858830/149‑1 (MQ=255)
                   ttCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTcgccg                                                                                                                            >  2:1390607/1‑149 (MQ=255)
                    tCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTcgccgc                                                                                                                           <  2:3540785/149‑1 (MQ=255)
                       gATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCa                                                                                                                        >  2:4464926/1‑149 (MQ=255)
                        aTTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGACATTGTCGCCGCTCAt                                                                                                                       >  2:1639687/1‑149 (MQ=255)
                        aTTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCAt                                                                                                                       <  2:3617005/149‑1 (MQ=255)
                        aTTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCAt                                                                                                                       <  2:1280344/149‑1 (MQ=255)
                         tttAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATg                                                                                                                      >  1:1163298/1‑149 (MQ=255)
                         tttAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATg                                                                                                                      <  2:2740368/149‑1 (MQ=255)
                                  aaaTCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCGCATGAAGcgcgga                                                                                                             >  1:4409447/1‑147 (MQ=255)
                                  aaaTCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTGGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTa                                                                                                             >  1:1642827/1‑149 (MQ=255)
                                     tCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTc                                                                                                          >  2:4142410/1‑149 (MQ=255)
                                     tCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAAGTGTCGCCGCTCATGAAGCGCGTATTc                                                                                                          >  1:4586337/1‑149 (MQ=255)
                                       cGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCg                                                                                                        <  1:491353/149‑1 (MQ=255)
                                        ggTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCgg                                                                                                       <  1:474807/149‑1 (MQ=255)
                                         gTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGt                                                                                                      <  1:3732922/149‑1 (MQ=255)
                                           cACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGa                                                                                                    >  2:1840111/1‑149 (MQ=255)
                                            aCGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGAt                                                                                                   >  2:2698685/1‑149 (MQ=255)
                                            aCGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGAt                                                                                                   <  1:2431513/149‑1 (MQ=255)
                                             cGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATc                                                                                                  <  1:4295531/149‑1 (MQ=255)
                                              gcagcaGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATcc                                                                                                 >  2:1367562/1‑149 (MQ=255)
                                                agcagTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgt                                                                                               >  1:1918648/1‑149 (MQ=255)
                                                   agTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgtcgt                                                                                            <  2:1451858/149‑1 (MQ=255)
                                                   agTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgtcgt                                                                                            <  2:1452004/149‑1 (MQ=255)
                                                   agTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgtcgt                                                                                            >  2:3541996/1‑149 (MQ=255)
                                                    gTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTg                                                                                           >  1:3740435/1‑149 (MQ=255)
                                                    gTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCAATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTg                                                                                           >  2:2740159/1‑149 (MQ=255)
                                                          gAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGaaa                                                                                     <  2:3950864/149‑1 (MQ=255)
                                                            aTCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATg                                                                                   >  1:919378/1‑149 (MQ=255)
                                                             tCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTGGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGCAATgc                                                                                  >  2:500497/1‑149 (MQ=255)
                                                              cTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgct                                                                                 <  2:2152902/149‑1 (MQ=255)
                                                                gggTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgctgg                                                                               <  1:2228363/149‑1 (MQ=255)
                                                                 ggTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgctggc                                                                              <  2:4766091/149‑1 (MQ=255)
                                                                 ggTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgctggc                                                                              <  1:1347492/149‑1 (MQ=255)
                                                                  gTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGGTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgctggct                                                                             <  2:1141177/149‑1 (MQ=255)
                                                                    gCGACGATTTATTTCAACCATCGCCCGGAAATCGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTTATCCGTCGTGCCGAAATGCTGGCTGa                                                                           >  2:1598457/1‑149 (MQ=255)
                                                                    gCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGa                                                                           <  2:69818/149‑1 (MQ=255)
                                                                     cgacgaTTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAg                                                                          <  2:3081930/149‑1 (MQ=255)
                                                                     cgacgaTTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAg                                                                          <  1:784071/149‑1 (MQ=255)
                                                                     cgacgaTTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAg                                                                          <  1:2345805/149‑1 (MQ=255)
                                                                      gacgaTTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGt                                                                         >  2:830325/1‑149 (MQ=255)
                                                                       acgaTTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt                                                                        >  2:1871324/1‑149 (MQ=255)
                                                                       acgaTTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt                                                                        >  2:420314/1‑149 (MQ=255)
                                                                       acgaTTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt                                                                        >  1:2639940/1‑149 (MQ=255)
                                                                         gaTTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTaa                                                                      <  1:2650422/149‑1 (MQ=255)
                                                                          atttatttCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAAt                                                                     <  2:254322/149‑1 (MQ=255)
                                                                          atttatttCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAAt                                                                     <  1:4148480/149‑1 (MQ=255)
                                                                              atttCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTGGTTTCCAGCGCGCTTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGt                                                                 >  2:2220745/1‑149 (MQ=255)
                                                                                ttCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGttt                                                               <  1:1446385/149‑1 (MQ=255)
                                                                                 tCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGtttt                                                              >  1:4617505/1‑149 (MQ=255)
                                                                                  cAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTc                                                             <  2:812027/149‑1 (MQ=255)
                                                                                      cATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTATGCAGATAAACCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAAATCTGGCTGAGTTAATTCGTTTTCGTCa                                                         >  1:2015485/1‑149 (MQ=255)
                                                                                      cATAGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCa                                                         >  2:1781040/1‑149 (MQ=255)
                                                                                       aTCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCAt                                                        >  2:143064/1‑149 (MQ=255)
                                                                                       aTCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCAt                                                        <  2:872095/149‑1 (MQ=255)
                                                                                         cGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAAACCGGAAATTGTCGCCGCGCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATTCGTTTTAGTCATgg                                                      >  1:1609497/1‑149 (MQ=255)
                                                                                         cGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTGGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCAGGAAGCGCGTATTCCCGGGATCCGGCGTGCCGCAATGCTGGCTGAGTTAATGCGTTTTCGTCAAgg                                                      >  1:3996435/1‑149 (MQ=255)
                                                                                             cGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTATGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATc                                                  >  1:1496089/1‑149 (MQ=255)
                                                                                                 aaaCGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCa                                              <  1:1148777/149‑1 (MQ=255)
                                                                                                  aaCGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCAt                                             >  1:3586683/1‑149 (MQ=255)
                                                                                                     gtacgtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAACTTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGc                                          >  1:1053065/1‑149 (MQ=255)
                                                                                                     gtacgtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTACGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATTCGTTTTCGTCATTGCATCGCCAtttc                                          >  1:3442323/1‑147 (MQ=255)
                                                                                                         gtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGa                                      <  1:2298276/149‑1 (MQ=255)
                                                                                                         gtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGa                                      >  2:4169333/1‑149 (MQ=255)
                                                                                                           gATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAAc                                    >  2:3915241/1‑149 (MQ=255)
                                                                                                            aTGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATTAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACg                                   >  1:2136290/1‑149 (MQ=255)
                                                                                                             tGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGc                                  <  1:32113/149‑1 (MQ=255)
                                                                                                                    gTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCaa                           >  2:299360/1‑149 (MQ=255)
                                                                                                                      ggTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCaaaa                         <  1:3397661/149‑1 (MQ=255)
                                                                                                                       gTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAAc                        <  2:73476/149‑1 (MQ=255)
                                                                                                                       gTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAAc                        >  2:1528748/1‑149 (MQ=255)
                                                                                                                       gTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAAc                        >  2:1430070/1‑149 (MQ=255)
                                                                                                                       gTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAAc                        <  1:182709/149‑1 (MQ=255)
                                                                                                                       gTCGTTGTTTCCAGCGCGATTTCAGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTACGGTTATACGTCGTGCCGAAATGCTGGCTGAGTTAATTCGTTTTCGTCATGGCATCGCCATTGCCGGAAAGCACGGCAAAAc                        >  2:2042956/1‑149 (MQ=255)
                                                                                                                           ttgttTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACa                    >  1:2566677/1‑149 (MQ=255)
                                                                                                                            tgttTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCTAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACaa                   <  2:4046266/149‑1 (MQ=255)
                                                                                                                              tttCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAACACGACAacc                 >  2:2868428/1‑149 (MQ=255)
                                                                                                                               ttCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAacca                <  1:969948/149‑1 (MQ=255)
                                                                                                                                tCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAaccac               <  2:3441289/149‑1 (MQ=255)
                                                                                                                                 ccAGCGCGATTTCTGCAGATAACCCGGACATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAaccacc              <  2:1556421/149‑1 (MQ=255)
                                                                                                                                 ccAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAaccacc              >  1:3450172/1‑149 (MQ=255)
                                                                                                                                   aGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACcgc            >  2:3369736/1‑149 (MQ=255)
                                                                                                                                      gcgATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGAt         <  2:1694146/149‑1 (MQ=255)
                                                                                                                                        gATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATgg       <  2:574805/149‑1 (MQ=255)
                                                                                                                                           ttCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGttt    <  1:3307485/149‑1 (MQ=255)
                                                                                                                                            tCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTc   <  2:3117837/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCC  >  CP009273/97376‑97665

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: