Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 601,847 A→T 94.5% 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,8470AT94.5% 568.5 / 18.1 181I331I (ATT→ATAybdKweak gamma‑glutamyl:cysteine ligase
Reads supporting (aligned to +/- strand):  ref base A (6/4);  new base T (67/104);  total (73/108)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 182 total aligned reads displayed.

CTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGG  >  CP009273/601702‑601991
                                                                                                                                                 |                                                                                                                                                
     aGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCgcgc                                                                                                                                           >  1:460719/1‑148 (MQ=255)
      ggCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCt                                                                                                                                          <  1:640548/148‑1 (MQ=255)
          gttCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGAt                                                                                                                                      >  2:1789388/2‑148 (MQ=255)
          aTTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGAt                                                                                                                                      <  2:932162/148‑1 (MQ=255)
          aTTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGAt                                                                                                                                      >  2:716731/1‑148 (MQ=255)
            tCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGc                                                                                                                                    <  1:1002908/148‑1 (MQ=255)
            tCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGc                                                                                                                                    >  2:2051030/1‑148 (MQ=255)
            tCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGc                                                                                                                                    >  2:1253704/1‑148 (MQ=255)
            tCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGc                                                                                                                                    <  1:2113687/148‑1 (MQ=255)
            tCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGc                                                                                                                                    <  1:1002907/148‑1 (MQ=255)
                ggCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCa                                                                                                                                >  2:428585/1‑148 (MQ=255)
                 gCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCaa                                                                                                                               >  2:3074085/1‑148 (MQ=255)
                     aGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAAtttt                                                                                                                           <  2:2708738/148‑1 (MQ=255)
                      gTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTa                                                                                                                          <  2:347956/148‑1 (MQ=255)
                       tCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAt                                                                                                                         <  1:1572928/148‑1 (MQ=255)
                        cACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAtg                                                                                                                        <  2:2450978/148‑1 (MQ=255)
                          ccGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAtgtg                                                                                                                      >  2:2914359/1‑148 (MQ=255)
                          ccGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAtgtg                                                                                                                      <  1:2776094/148‑1 (MQ=255)
                          ccGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAtgtg                                                                                                                      >  2:2090969/1‑148 (MQ=255)
                           cGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGc                                                                                                                     >  1:286116/1‑148 (MQ=255)
                              cccAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGa                                                                                                                  <  2:3116885/148‑1 (MQ=255)
                               ccAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGaa                                                                                                                 <  2:637614/148‑1 (MQ=255)
                               ccAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGaa                                                                                                                 <  1:866289/148‑1 (MQ=255)
                               ccAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGaa                                                                                                                 <  1:860720/148‑1 (MQ=255)
                                cAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAg                                                                                                                <  1:1527503/148‑1 (MQ=255)
                                cAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAg                                                                                                                >  2:422038/1‑148 (MQ=255)
                                cAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAg                                                                                                                >  2:726473/1‑148 (MQ=255)
                                cAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAg                                                                                                                <  1:1527473/148‑1 (MQ=255)
                                 aGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGaacg                                                                                                               <  2:2844653/148‑3 (MQ=255)
                                 aGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAgg                                                                                                               <  2:1256245/148‑1 (MQ=255)
                                 aGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAgg                                                                                                               <  2:885892/148‑1 (MQ=255)
                                  gATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAggg                                                                                                              >  2:3154847/1‑148 (MQ=255)
                                  gATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAggg                                                                                                              >  2:1600472/1‑148 (MQ=255)
                                     ctcACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCg                                                                                                           >  1:2373067/1‑148 (MQ=255)
                                     ctcACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCg                                                                                                           >  1:2373148/1‑148 (MQ=255)
                                     ctcACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCg                                                                                                           <  2:1551028/148‑1 (MQ=255)
                                       cacaATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAt                                                                                                         >  1:2988922/1‑148 (MQ=255)
                                       cacaATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAt                                                                                                         >  1:3106498/1‑148 (MQ=255)
                                        acaATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAtt                                                                                                        >  1:20085/1‑148 (MQ=255)
                                         caATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAttt                                                                                                       <  1:1604426/148‑1 (MQ=255)
                                         caATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAttt                                                                                                       <  1:1604519/148‑1 (MQ=255)
                                          aaTGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAtttt                                                                                                      >  1:1179585/1‑148 (MQ=255)
                                           aTGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCCCCAATTTTCTGTGCGGAAGGGGCGAttttt                                                                                                     <  2:1125778/148‑1 (MQ=255)
                                           aTGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAttttt                                                                                                     >  2:724659/1‑148 (MQ=255)
                                            tGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGGAGGGGCGAtttttt                                                                                                    <  2:764741/148‑1 (MQ=255)
                                             gCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTc                                                                                                   <  2:384211/148‑1 (MQ=255)
                                              cTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTcc                                                                                                  <  2:825119/148‑1 (MQ=255)
                                              cTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTcc                                                                                                  <  2:2523283/148‑1 (MQ=255)
                                              cTTGTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTCTGTGCGGAAGGGGCGATTTTTTcc                                                                                                  <  2:2483872/148‑1 (MQ=255)
                                                  ttCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCa                                                                                              <  1:243759/148‑1 (MQ=255)
                                                       cAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGc                                                                                         >  1:2134331/1‑148 (MQ=255)
                                                         gCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCaa                                                                                       >  2:1394158/1‑148 (MQ=255)
                                                          cccAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAg                                                                                      >  1:56406/1‑148 (MQ=255)
                                                          cccAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAg                                                                                      >  1:583577/1‑148 (MQ=255)
                                                            cAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGt                                                                                    >  1:1818836/1‑148 (MQ=255)
                                                            cAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGt                                                                                    <  1:3182505/148‑1 (MQ=255)
                                                            cAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGt                                                                                    <  1:3301979/148‑1 (MQ=255)
                                                             aaTCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTa                                                                                   >  2:2638335/1‑148 (MQ=255)
                                                              aTCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTAt                                                                                  <  2:2387770/148‑1 (MQ=255)
                                                               tCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATc                                                                                 <  2:643650/148‑1 (MQ=255)
                                                               tCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATc                                                                                 <  2:3000125/148‑1 (MQ=255)
                                                               tCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATc                                                                                 <  1:2774688/148‑1 (MQ=255)
                                                                    gAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCg                                                                            <  2:1366744/148‑1 (MQ=255)
                                                                    gAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCg                                                                            <  2:2707579/148‑1 (MQ=255)
                                                                      gccgccATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGtt                                                                          >  1:2458561/1‑148 (MQ=255)
                                                                       ccgccATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTa                                                                         >  1:2164979/1‑148 (MQ=255)
                                                                       ccgccATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTa                                                                         >  1:2645507/1‑148 (MQ=255)
                                                                           cATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCgg                                                                     <  1:1618356/148‑1 (MQ=255)
                                                                            aTCGTCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGt                                                                    <  1:900019/148‑1 (MQ=255)
                                                                             tCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTc                                                                   >  2:2258306/1‑148 (MQ=255)
                                                                             tCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTc                                                                   >  2:596584/1‑148 (MQ=255)
                                                                                   acgaAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGAt                                                             <  2:2006403/148‑1 (MQ=255)
                                                                                      aaaTCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGAtctc                                                          <  2:2804799/148‑1 (MQ=255)
                                                                                           gcgcATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAgtg                                                     <  2:2134331/148‑1 (MQ=255)
                                                                                            cgcATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAgtgt                                                    <  1:226489/148‑1 (MQ=255)
                                                                                            cgcATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAgtgt                                                    <  1:226525/148‑1 (MQ=255)
                                                                                              cATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGc                                                  <  1:1521970/148‑1 (MQ=255)
                                                                                              cATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGc                                                  <  1:1521350/148‑1 (MQ=255)
                                                                                                 tAGCTGCGCTTCTTTCAGACCGCTGACGCGCTGTCGAGGCAGGGGCTGTATCGCGCTCGATGCACCAATTTGATNTGCTGAAGGGGGGATTTTTTCCAGCAATCGGAAGGTATCTTCCGTTAGAGGTCGACGATCTCCAGTGTGCGGa                                               <  1:1679200/148‑1 (MQ=255)
                                                                                                  aGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGAt                                              >  2:818676/1‑148 (MQ=255)
                                                                                                      gcgcTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGt                                          <  1:1507696/148‑1 (MQ=255)
                                                                                                       cgcTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtg                                         <  2:330260/148‑1 (MQ=255)
                                                                                                       cgcTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtg                                         <  2:1847772/148‑1 (MQ=255)
                                                                                                       cgcTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtg                                         <  1:170864/148‑1 (MQ=255)
                                                                                                               tcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGcc                                 <  1:1253704/148‑1 (MQ=255)
                                                                                                                  gACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTc                              >  2:3305914/1‑148 (MQ=255)
                                                                                                                  gACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGGATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTc                              >  1:2516551/1‑148 (MQ=255)
                                                                                                                    ccGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCaa                            <  1:422038/148‑1 (MQ=255)
                                                                                                                    ccGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCaa                            >  2:1863568/1‑148 (MQ=255)
                                                                                                                               cTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGa                 <  2:514772/148‑1 (MQ=255)
                                                                                                                               cTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGa                 <  2:1362759/148‑1 (MQ=255)
                                                                                                                                tGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGAc                <  1:52988/148‑1 (MQ=255)
                                                                                                                                 ggCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACa               >  2:658515/1‑148 (MQ=255)
                                                                                                                                   cGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAgg             <  2:3201944/148‑1 (MQ=255)
                                                                                                                                    gATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCACTCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGc            <  1:2647552/148‑1 (MQ=255)
                                                                                                                                     aTGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGcc           <  1:684788/148‑1 (MQ=255)
                                                                                                                                           ggCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGaaa     >  2:30454/1‑148 (MQ=255)
                                                                                                                                             ccTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACg   <  2:2140849/148‑1 (MQ=255)
                                                                                                                                             ccTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACg   >  1:3053370/1‑148 (MQ=255)
                                                                                                                                              ctctATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACgg  <  2:2882578/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGG  >  CP009273/601702‑601991

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: