Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 601,847 A→T 100% I331I (ATT→ATA)  ybdK ← weak gamma‑glutamyl:cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273601,8470AT98.6% 932.5 / 1.1 281I331I (ATT→ATA) ybdKweak gamma‑glutamyl:cysteine ligase
Reads supporting (aligned to +/- strand):  ref base A (0/4);  new base T (139/138);  total (139/142)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-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 285 total aligned reads displayed.

TCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGT  >  CP009273/601701‑601992
                                                                                                                                                  |                                                                                                                                                 
  tGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATcg                                                                                                                                               <  1:1687549/149‑1 (MQ=255)
   gTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATcgc                                                                                                                                              >  2:1564938/1‑149 (MQ=255)
     aaGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCgcgc                                                                                                                                            >  2:2515436/1‑149 (MQ=255)
     aaGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCgcgc                                                                                                                                            >  2:2456416/1‑149 (MQ=255)
     aaGGCAATTCACGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCgcgc                                                                                                                                            >  1:2069630/1‑149 (MQ=255)
      aGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCt                                                                                                                                           >  1:1793646/1‑149 (MQ=255)
        gCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCg                                                                                                                                         >  1:1959162/1‑149 (MQ=255)
          aaTTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGAt                                                                                                                                       <  2:2974848/149‑1 (MQ=255)
          aaTTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGAt                                                                                                                                       <  1:605992/149‑1 (MQ=255)
             tCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCa                                                                                                                                    <  1:2868929/149‑1 (MQ=255)
                  gCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAAt                                                                                                                               <  2:684108/149‑1 (MQ=255)
                    ttAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAAttt                                                                                                                             >  1:1492458/1‑149 (MQ=255)
                      aGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTa                                                                                                                           <  1:1327047/149‑1 (MQ=255)
                       gTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAt                                                                                                                          <  1:1513002/149‑1 (MQ=255)
                       gTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAt                                                                                                                          <  1:2359525/149‑1 (MQ=255)
                        tCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATAGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAtg                                                                                                                         <  2:3687496/149‑1 (MQ=255)
                         cACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAtgt                                                                                                                        <  2:947432/149‑1 (MQ=255)
                           ccGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGc                                                                                                                      >  1:39139/1‑149 (MQ=255)
                            cGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCg                                                                                                                     <  2:3478278/149‑1 (MQ=255)
                             ggCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCgg                                                                                                                    >  2:1196127/1‑149 (MQ=255)
                             ggCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCgg                                                                                                                    >  2:143177/1‑149 (MQ=255)
                                ccAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAg                                                                                                                 >  1:3992927/1‑149 (MQ=255)
                                 cAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAgg                                                                                                                >  1:752445/1‑149 (MQ=255)
                                 cAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAgg                                                                                                                >  2:2855529/1‑149 (MQ=255)
                                     tctcACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCg                                                                                                            >  1:853457/1‑149 (MQ=255)
                                      ctcACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGa                                                                                                           <  1:3120055/149‑1 (MQ=255)
                                         acaATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAttt                                                                                                        <  1:37113/149‑1 (MQ=255)
                                          caATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAtttt                                                                                                       <  2:1695528/149‑1 (MQ=255)
                                             tGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTc                                                                                                    <  2:3545314/149‑1 (MQ=255)
                                              gCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTcc                                                                                                   <  1:2162186/149‑1 (MQ=255)
                                              gCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTcc                                                                                                   <  1:9860/149‑1 (MQ=255)
                                                tttttCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAg                                                                                                 >  1:1632947/1‑149 (MQ=255)
                                                 ttttCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGc                                                                                                >  2:1262574/1‑149 (MQ=255)
                                                 ttttCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGc                                                                                                >  1:225379/1‑149 (MQ=255)
                                                   ttCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCaa                                                                                              <  1:1207495/149‑1 (MQ=255)
                                                    tCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAAt                                                                                             <  1:3868785/149‑1 (MQ=255)
                                                     cACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATTGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATc                                                                                            <  2:3418262/149‑1 (MQ=255)
                                                          gCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTCTCGCGCTCGATGCACCAATGTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAg                                                                                       >  2:2836537/1‑149 (MQ=255)
                                                          gCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAg                                                                                       >  2:776278/1‑149 (MQ=255)
                                                           cccAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAgg                                                                                      >  2:2671235/1‑149 (MQ=255)
                                                             cAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTa                                                                                    <  1:845432/149‑1 (MQ=255)
                                                             cAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTa                                                                                    >  1:821709/1‑149 (MQ=255)
                                                               aTCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATc                                                                                  <  2:1477488/149‑1 (MQ=255)
                                                                 cAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCtt                                                                                <  2:3436170/149‑1 (MQ=255)
                                                                  aGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTc                                                                               <  2:2519005/149‑1 (MQ=255)
                                                                   gCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTcc                                                                              >  1:1003692/1‑149 (MQ=255)
                                                                   gCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTcc                                                                              >  1:584452/1‑149 (MQ=255)
                                                                    cGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCg                                                                             <  1:654822/149‑1 (MQ=255)
                                                                     gAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGt                                                                            <  1:3564208/149‑1 (MQ=255)
                                                                       gccgccATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTa                                                                          >  1:2608502/1‑149 (MQ=255)
                                                                       gccgccATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTa                                                                          >  2:2141210/1‑149 (MQ=255)
                                                                       gccgccATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTa                                                                          >  2:3774325/1‑149 (MQ=255)
                                                                       gccgccATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTa                                                                          >  1:1532598/1‑149 (MQ=255)
                                                                        ccgccATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAg                                                                         >  2:1191673/1‑149 (MQ=255)
                                                                           ccATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGAgg                                                                      >  2:3796095/1‑149 (MQ=255)
                                                                            cATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGt                                                                     >  1:3956149/1‑149 (MQ=255)
                                                                             aTCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTc                                                                    >  1:2531418/1‑149 (MQ=255)
                                                                             aTCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTc                                                                    >  1:1097178/1‑149 (MQ=255)
                                                                              tCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTcg                                                                   >  2:2197995/1‑149 (MQ=255)
                                                                               cGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTcga                                                                  >  1:4064168/1‑149 (MQ=255)
                                                                                   gacgaAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGAt                                                              >  1:475937/1‑149 (MQ=255)
                                                                                   gacgaAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGAt                                                              >  2:2783693/1‑149 (MQ=255)
                                                                                       aaaTCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTcc                                                          >  1:3521389/1‑149 (MQ=255)
                                                                                       aaaTCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTcc                                                          <  1:3598924/149‑1 (MQ=255)
                                                                                       aaaTCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTACAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTcc                                                          >  1:22091/1‑149 (MQ=255)
                                                                                         aTCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAg                                                        <  2:3300817/149‑1 (MQ=255)
                                                                                                aTTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCgg                                                 >  2:459941/1‑149 (MQ=255)
                                                                                                 ttAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTACAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGa                                                >  1:154264/1‑149 (MQ=255)
                                                                                                 ttAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGa                                                <  2:1027149/149‑1 (MQ=255)
                                                                                                  tAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGAt                                               <  1:3625485/149‑1 (MQ=255)
                                                                                                    gCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCg                                             >  2:2947888/1‑149 (MQ=255)
                                                                                                      tGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGt                                           <  1:116269/149‑1 (MQ=255)
                                                                                                       gcgcTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtg                                          >  2:175926/1‑149 (MQ=255)
                                                                                                           ttcattcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtgatga                                      <  1:1014219/149‑1 (MQ=255)
                                                                                                           ttcattcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTATATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtgatga                                      <  2:2069630/149‑1 (MQ=255)
                                                                                                            tcattcaGACAGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGATCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTACAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGAc                                     >  2:3165934/1‑149 (MQ=255)
                                                                                                             cattcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACg                                    >  2:630423/1‑149 (MQ=255)
                                                                                                             cattcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACg                                    >  1:399521/1‑149 (MQ=255)
                                                                                                             cattcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTAAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACg                                    >  2:2250670/1‑149 (MQ=255)
                                                                                                                  aGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTc                               >  1:2994663/1‑149 (MQ=255)
                                                                                                                         tGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCa                        <  1:1414682/149‑1 (MQ=255)
                                                                                                                              acCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCg                   >  1:230366/1‑149 (MQ=255)
                                                                                                                               ccTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGa                  <  1:2456416/149‑1 (MQ=255)
                                                                                                                               ccTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGa                  <  1:2515436/149‑1 (MQ=255)
                                                                                                                                cTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGAc                 <  2:2307397/149‑1 (MQ=255)
                                                                                                                                 tGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACa                <  2:3992927/149‑1 (MQ=255)
                                                                                                                                 tGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTACAGTGTGCGGATAGGTGATGACGCCTTCAAGCCCATAGCGCCa                >  1:187414/1‑149 (MQ=255)
                                                                                                                                  ggCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTTTGCGGATCGGTGATGACGCCTTCAAGCCCATAGAGACAg               >  1:3377123/1‑149 (MQ=255)
                                                                                                                                  ggCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAg               >  1:3431858/1‑149 (MQ=255)
                                                                                                                                   gCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAgg              >  1:2477140/1‑149 (MQ=255)
                                                                                                                                     gATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGcc            >  2:1928724/1‑149 (MQ=255)
                                                                                                                                     gATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGcc            <  1:3814159/149‑1 (MQ=255)
                                                                                                                                        gCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTgg         <  2:2608502/149‑1 (MQ=255)
                                                                                                                                           gggCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGaaa      <  1:2783693/149‑1 (MQ=255)
                                                                                                                                           gggCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGaaa      >  2:168233/1‑149 (MQ=255)
                                                                                                                                             gCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACg    <  2:1176831/149‑1 (MQ=255)
                                                                                                                                              ccTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACgg   <  1:3091017/149‑1 (MQ=255)
                                                                                                                                              ccTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACgg   >  1:3755340/1‑149 (MQ=255)
                                                                                                                                               ctctATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGt  >  2:3922372/1‑149 (MQ=255)
                                                                                                                                               ctctATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGt  >  1:215149/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGT  >  CP009273/601701‑601992

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: ‑