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,8470AT97.7% 716.1 / 4.0 219I331I (ATT→ATA) ybdKweak gamma‑glutamyl:cysteine ligase
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base T (110/104);  total (112/107)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 221 total aligned reads displayed.

CTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGG  >  CP009273/601702‑601991
                                                                                                                                                 |                                                                                                                                                
ctGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTAt                                                                                                                                                <  2:2929602/148‑1 (MQ=255)
  gTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATcg                                                                                                                                              >  2:1503493/1‑148 (MQ=255)
    aaGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATcgcg                                                                                                                                            <  2:529294/148‑1 (MQ=255)
     aGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCgcgc                                                                                                                                           <  1:2890856/148‑1 (MQ=255)
     aGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCgcgc                                                                                                                                           >  2:1141820/1‑148 (MQ=255)
        cAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATAGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCg                                                                                                                                        >  1:2670246/1‑148 (MQ=255)
          aTTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGAt                                                                                                                                      >  2:3532535/1‑148 (MQ=255)
            tCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGc                                                                                                                                    >  2:795098/1‑148 (MQ=255)
            tCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGATGCGCTTCATTCAGACCGCTGACGACCTGGAGATGCAGGGCCTCTATCGCGCTCGATGc                                                                                                                                    <  2:32570/148‑1 (MQ=255)
               cGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCAcc                                                                                                                                 >  2:3178451/1‑148 (MQ=255)
                ggCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCa                                                                                                                                <  2:2949067/148‑1 (MQ=255)
                     aGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAAtttt                                                                                                                           >  1:1590609/1‑148 (MQ=255)
                     aGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAAtttt                                                                                                                           >  1:3432289/1‑148 (MQ=255)
                       tCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAt                                                                                                                         <  1:3299527/148‑1 (MQ=255)
                        cACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAtg                                                                                                                        >  1:1368981/1‑148 (MQ=255)
                         aCCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAtgt                                                                                                                       >  2:1960621/1‑148 (MQ=255)
                          ccGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAtgtg                                                                                                                      <  1:3184581/148‑1 (MQ=255)
                          ccGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAtgtg                                                                                                                      >  1:663464/1‑148 (MQ=255)
                               ccAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGaa                                                                                                                 >  2:2574884/1‑148 (MQ=255)
                                cAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAg                                                                                                                <  1:3452382/148‑1 (MQ=255)
                                  gATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAggg                                                                                                              >  1:1947247/1‑148 (MQ=255)
                                  gATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGACGCAGGGCCTCTATCGCGCTCGATGAACCAATTTTATGTGCGGAAggg                                                                                                              >  2:1643170/1‑148 (MQ=255)
                                   aTCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAgggg                                                                                                             >  1:2214633/1‑148 (MQ=255)
                                   aTCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAgggg                                                                                                             >  1:700733/1‑148 (MQ=255)
                                   aTCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAgggg                                                                                                             >  1:1645317/1‑148 (MQ=255)
                                    tctcACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGACGGGGc                                                                                                            >  1:3191355/1‑148 (MQ=255)
                                        acaATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAtt                                                                                                        <  2:1910495/148‑1 (MQ=255)
                                        acaATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAtt                                                                                                        >  1:3242777/1‑148 (MQ=255)
                                        acaATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAtt                                                                                                        >  2:1188404/1‑148 (MQ=255)
                                        acaATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAtt                                                                                                        >  2:1464198/1‑148 (MQ=255)
                                             gCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTc                                                                                                   >  2:3064641/1‑148 (MQ=255)
                                             gCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTc                                                                                                   >  2:393860/1‑148 (MQ=255)
                                             gCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTc                                                                                                   >  2:393979/1‑148 (MQ=255)
                                             gCTTGTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTc                                                                                                   <  1:2727014/148‑1 (MQ=255)
                                              cTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTcc                                                                                                  <  2:3294351/148‑1 (MQ=255)
                                                ttttCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAg                                                                                                >  2:2409632/1‑148 (MQ=255)
                                                ttttCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAg                                                                                                >  1:538703/1‑148 (MQ=255)
                                                   tCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCaa                                                                                             <  1:1234974/148‑1 (MQ=255)
                                                    cACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAAt                                                                                            <  2:3240277/148‑1 (MQ=255)
                                                        aGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCa                                                                                        <  2:1840059/148‑1 (MQ=255)
                                                        aGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCa                                                                                        <  2:553185/148‑1 (MQ=255)
                                                         gCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCaa                                                                                       >  2:2180956/1‑148 (MQ=255)
                                                             aaTCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTa                                                                                   <  1:3117149/148‑1 (MQ=255)
                                                             aaTCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTa                                                                                   >  1:2379450/1‑148 (MQ=255)
                                                              aTCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTAt                                                                                  <  1:1435957/148‑1 (MQ=255)
                                                                 aGCGAGCCGCCATCGGCGACGAAATCGCGCATGAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCtt                                                                               <  2:888498/148‑1 (MQ=255)
                                                                   cGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTcc                                                                             >  1:1939704/1‑148 (MQ=255)
                                                                   cGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTcc                                                                             >  2:680215/1‑148 (MQ=255)
                                                                   cGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTcc                                                                             >  1:1733123/1‑148 (MQ=255)
                                                                   cGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTcc                                                                             >  1:1732972/1‑148 (MQ=255)
                                                                    gAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCg                                                                            <  1:630389/148‑1 (MQ=255)
                                                                     aGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGt                                                                           >  1:1172020/1‑148 (MQ=255)
                                                                       ccgccATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTa                                                                         <  2:717458/148‑1 (MQ=255)
                                                                       ccgccATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTa                                                                         >  2:3120451/1‑148 (MQ=255)
                                                                          ccATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCg                                                                      >  1:2022967/1‑148 (MQ=255)
                                                                           cATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCgg                                                                     >  1:1056838/1‑148 (MQ=255)
                                                                           cATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCgg                                                                     <  2:694736/148‑1 (MQ=255)
                                                                            aTCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGAAATCGCAAGGTATCTTCCGTTAGCGGt                                                                    <  1:377652/148‑1 (MQ=255)
                                                                             tCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTc                                                                   >  1:508035/1‑148 (MQ=255)
                                                                             tCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTc                                                                   >  2:2335752/1‑148 (MQ=255)
                                                                             tCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTc                                                                   >  2:278796/1‑148 (MQ=255)
                                                                                gCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTcgac                                                                >  2:1655997/1‑148 (MQ=255)
                                                                                   acgaAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGAt                                                             >  2:1727048/1‑148 (MQ=255)
                                                                                   acgaAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGAt                                                             <  2:758186/148‑1 (MQ=255)
                                                                                   acgaAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGAt                                                             >  1:2730417/1‑148 (MQ=255)
                                                                                    cgaAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGAtc                                                            <  2:912031/148‑1 (MQ=255)
                                                                                    cgaAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGAtc                                                            >  2:870816/1‑148 (MQ=255)
                                                                                     gaAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGAtct                                                           <  2:3331575/148‑1 (MQ=255)
                                                                                        aTCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCa                                                        <  2:2310074/148‑1 (MQ=255)
                                                                                              cATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGc                                                  >  2:3423262/1‑148 (MQ=255)
                                                                                               aTTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCg                                                 >  1:2356416/1‑148 (MQ=255)
                                                                                               aTTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCg                                                 >  1:2428796/1‑148 (MQ=255)
                                                                                                 tAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGa                                               <  2:2735433/148‑1 (MQ=255)
                                                                                                      gcgcTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATATTACGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGt                                          >  1:759053/1‑148 (MQ=255)
                                                                                                      gcgcTTCATTCAGAACGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGt                                          >  1:3405465/1‑148 (MQ=255)
                                                                                                       cgcTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtg                                         >  2:1563829/1‑148 (MQ=255)
                                                                                                        gcTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtga                                        <  2:3023088/148‑1 (MQ=255)
                                                                                                         cTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtgat                                       <  2:2842982/148‑1 (MQ=255)
                                                                                                           tcattcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtgatga                                     >  2:2276497/1‑148 (MQ=255)
                                                                                                           tcattcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtgatga                                     <  1:1930643/148‑1 (MQ=255)
                                                                                                              ttcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGc                                  <  1:1750467/148‑1 (MQ=255)
                                                                                                              ttcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGc                                  <  1:683068/148‑1 (MQ=255)
                                                                                                              ttcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGc                                  <  2:664001/148‑1 (MQ=255)
                                                                                                              ttcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGc                                  <  1:2162205/148‑1 (MQ=255)
                                                                                                              ttcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTCTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGc                                  <  1:2246805/148‑1 (MQ=255)
                                                                                                                        tGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGccc                        <  1:897361/148‑1 (MQ=255)
                                                                                                                         gacgacCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCa                       >  1:3371150/1‑148 (MQ=255)
                                                                                                                         gacgacCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCa                       <  2:2454915/148‑1 (MQ=255)
                                                                                                                                 ggCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACa               <  2:3399273/148‑1 (MQ=255)
                                                                                                                                  gCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAg              <  2:3386423/148‑1 (MQ=255)
                                                                                                                                   cGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAgg             <  1:2745029/148‑1 (MQ=255)
                                                                                                                                    gATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGc            <  1:2518992/148‑1 (MQ=255)
                                                                                                                                    gATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGc            <  1:2519229/148‑1 (MQ=255)
                                                                                                                                    gATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGc            <  2:1579151/148‑1 (MQ=255)
                                                                                                                                    gATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGc            <  2:2484718/148‑1 (MQ=255)
                                                                                                                                     aTGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGAATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTTATGACGCCTTCAAGCCCATATCGACATGcc           >  1:825658/1‑148 (MQ=255)
                                                                                                                                         aGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGa       >  2:3351559/1‑148 (MQ=255)
                                                                                                                                         aGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGa       <  1:2817734/148‑1 (MQ=255)
                                                                                                                                             ccTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACg   >  2:2099418/1‑148 (MQ=255)
                                                                                                                                              ctctATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACgg  <  1:2289062/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: ‑