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.0% 652.3 / 7.3 201I331I (ATT→ATA) ybdKweak gamma‑glutamyl:cysteine ligase
Reads supporting (aligned to +/- strand):  ref base A (2/4);  new base T (86/109);  total (88/113)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 202 total aligned reads displayed.

TGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGG  >  CP009273/601703‑601991
                                                                                                                                                |                                                                                                                                                
tGTAAGGCAAGTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATc                                                                                                                                               <  2:2282763/148‑1 (MQ=255)
 gTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATcg                                                                                                                                              >  2:2600158/1‑148 (MQ=255)
   aaGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATcgcg                                                                                                                                            >  2:40253/1‑148 (MQ=255)
     ggCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCt                                                                                                                                          <  1:1416275/148‑1 (MQ=255)
      gCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTc                                                                                                                                         >  2:1343570/1‑148 (MQ=255)
        aaTTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGa                                                                                                                                       <  1:2261213/148‑1 (MQ=255)
        aaTTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGa                                                                                                                                       <  1:2261166/148‑1 (MQ=255)
          ttCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATg                                                                                                                                     >  2:1227328/1‑148 (MQ=255)
               ggCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCa                                                                                                                                >  1:3164978/1‑148 (MQ=255)
               ggCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCa                                                                                                                                >  2:2904458/1‑148 (MQ=255)
                   tAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAAttt                                                                                                                            >  2:1256002/1‑148 (MQ=255)
                     gTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTa                                                                                                                          <  2:492342/148‑1 (MQ=255)
                      tCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAt                                                                                                                         >  2:1137358/1‑148 (MQ=255)
                       cACCGGCCCAGATCTCCCAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAtg                                                                                                                        <  1:1000756/148‑1 (MQ=255)
                       cACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAtg                                                                                                                        <  2:3369080/148‑1 (MQ=255)
                              ccAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGaa                                                                                                                 <  2:3795843/148‑1 (MQ=255)
                               cAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAg                                                                                                                <  2:1246221/148‑1 (MQ=255)
                               cAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAg                                                                                                                >  2:567025/1‑148 (MQ=255)
                               cAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAg                                                                                                                >  2:538204/1‑148 (MQ=255)
                                 gATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAggg                                                                                                              >  2:73815/1‑148 (MQ=255)
                                 gATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAggg                                                                                                              <  1:1885519/148‑1 (MQ=255)
                                  aTCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAgggg                                                                                                             <  1:43496/148‑1 (MQ=255)
                                  aTCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAgggg                                                                                                             <  1:3463480/148‑1 (MQ=255)
                                   tctcACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGc                                                                                                            >  1:1173595/1‑148 (MQ=255)
                                     tcACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGa                                                                                                          <  2:1988639/148‑1 (MQ=255)
                                     tcACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGa                                                                                                          <  2:2165697/148‑1 (MQ=255)
                                      cacaATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAt                                                                                                         <  2:651955/148‑1 (MQ=255)
                                      cacaATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAt                                                                                                         <  2:2743404/148‑1 (MQ=255)
                                        caATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAttt                                                                                                       <  1:3203477/148‑1 (MQ=255)
                                        caATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAttt                                                                                                       >  1:639928/1‑148 (MQ=255)
                                        caATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAttt                                                                                                       >  1:2463947/1‑148 (MQ=255)
                                        caATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAttt                                                                                                       >  1:1733204/1‑148 (MQ=255)
                                        caATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAttt                                                                                                       <  1:1815978/148‑1 (MQ=255)
                                          aTGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAttttt                                                                                                     >  1:3578035/1‑148 (MQ=255)
                                          aTGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAttttt                                                                                                     >  2:2934653/1‑148 (MQ=255)
                                           tGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAtttttt                                                                                                    <  1:2583130/148‑1 (MQ=255)
                                            gCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTc                                                                                                   <  1:2387379/148‑1 (MQ=255)
                                               ttttCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAg                                                                                                >  1:2595483/1‑148 (MQ=255)
                                               ttttCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAg                                                                                                <  1:3662539/148‑1 (MQ=255)
                                                 ttCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCa                                                                                              <  2:3100787/148‑1 (MQ=255)
                                                 ttCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCa                                                                                              <  1:2212070/148‑1 (MQ=255)
                                                  tCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCaa                                                                                             <  1:2529279/148‑1 (MQ=255)
                                                         cccAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAg                                                                                      >  2:265718/1‑148 (MQ=255)
                                                         cccAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAg                                                                                      >  2:266988/1‑148 (MQ=255)
                                                          ccAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAgg                                                                                     <  1:1781261/148‑1 (MQ=255)
                                                           cAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCACTCGCAAGGt                                                                                    >  1:1998601/1‑148 (MQ=255)
                                                           cAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGt                                                                                    >  2:3525242/1‑148 (MQ=255)
                                                             aTCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCCGCAATCGCAAGGTAt                                                                                  >  1:115215/1‑148 (MQ=255)
                                                             aTCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTAt                                                                                  <  1:3355058/148‑1 (MQ=255)
                                                              tccgCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATc                                                                                 <  1:1704491/145‑1 (MQ=255)
                                                              tCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATc                                                                                 >  2:1077958/1‑148 (MQ=255)
                                                              tCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATc                                                                                 >  2:1077907/1‑148 (MQ=255)
                                                              tCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATc                                                                                 >  2:224258/1‑148 (MQ=255)
                                                               cAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCt                                                                                <  1:2171535/148‑1 (MQ=255)
                                                                 gCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTc                                                                              <  1:378499/148‑1 (MQ=255)
                                                                 gCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTc                                                                              <  1:3557876/148‑1 (MQ=255)
                                                                   gAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCg                                                                            <  1:2110436/148‑1 (MQ=255)
                                                                   gAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCg                                                                            <  1:3378877/148‑1 (MQ=255)
                                                                   gAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCg                                                                            <  2:1399810/148‑1 (MQ=255)
                                                                       cgccATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAg                                                                        >  2:1733825/1‑148 (MQ=255)
                                                                         ccATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCg                                                                      >  2:3295076/1‑148 (MQ=255)
                                                                           aTCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGt                                                                    <  1:1092161/148‑1 (MQ=255)
                                                                                cgacgaAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTcgacg                                                               <  2:2449164/148‑1 (MQ=255)
                                                                                cgacgaAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTcgacg                                                               <  1:57286/148‑1 (MQ=255)
                                                                                  acgaAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGAt                                                             >  1:2694728/1‑148 (MQ=255)
                                                                                  acgaAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGAt                                                             >  1:778202/1‑148 (MQ=255)
                                                                                      aaTCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTcc                                                         >  2:2803809/1‑148 (MQ=255)
                                                                                       aTCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCa                                                        >  2:1534309/1‑148 (MQ=255)
                                                                                       aTCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCa                                                        <  1:1213748/148‑1 (MQ=255)
                                                                                       aTCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCa                                                        >  2:3025411/1‑148 (MQ=255)
                                                                                         cgcgCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAgt                                                      >  2:3637919/1‑148 (MQ=255)
                                                                                               ttAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCgg                                                >  1:526012/1‑148 (MQ=255)
                                                                                                     gcgcTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGt                                          <  1:2600158/148‑1 (MQ=255)
                                                                                                      cgcTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtg                                         <  1:463559/148‑1 (MQ=255)
                                                                                                      cgcTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtg                                         <  1:2656947/148‑1 (MQ=255)
                                                                                                      cgcTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtg                                         <  1:2803809/148‑1 (MQ=255)
                                                                                                      cgcTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtg                                         <  1:740740/148‑1 (MQ=255)
                                                                                                       gcTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtga                                        <  2:3584069/148‑1 (MQ=255)
                                                                                                        cTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtgat                                       <  2:1452993/148‑1 (MQ=255)
                                                                                                        cTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtgat                                       <  2:1453054/148‑1 (MQ=255)
                                                                                                              tcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGcc                                 >  2:1502639/1‑148 (MQ=255)
                                                                                                              tcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGcc                                 >  1:3174318/1‑148 (MQ=255)
                                                                                                                 gACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTc                              <  1:3261974/148‑1 (MQ=255)
                                                                                                                        gacgacCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCa                       <  1:1727625/148‑1 (MQ=255)
                                                                                                                        gacgacCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCa                       <  1:996418/148‑1 (MQ=255)
                                                                                                                         acgacCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCAt                      <  1:2039453/148‑1 (MQ=255)
                                                                                                                          cgacCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATa                     >  2:413142/1‑148 (MQ=255)
                                                                                                                           gacCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAg                    <  1:3043745/148‑1 (MQ=255)
                                                                                                                               tGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGAc                <  1:846296/148‑1 (MQ=255)
                                                                                                                               tGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGAc                >  1:1545619/1‑148 (MQ=255)
                                                                                                                               tGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGAc                >  1:1226876/1‑148 (MQ=255)
                                                                                                                               tGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGAc                <  1:887277/148‑1 (MQ=255)
                                                                                                                                   gATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGc            <  2:2213895/148‑1 (MQ=255)
                                                                                                                                   gATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGc            <  2:2385811/148‑1 (MQ=255)
                                                                                                                                    aTGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGcc           >  2:2261533/1‑148 (MQ=255)
                                                                                                                                    aTGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGcc           <  2:3237908/148‑1 (MQ=255)
                                                                                                                                     tGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCt          <  1:2004892/148‑1 (MQ=255)
                                                                                                                                         gggCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGaa      <  1:3025411/148‑1 (MQ=255)
                                                                                                                                          ggCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGaaa     <  1:1227328/148‑1 (MQ=255)
                                                                                                                                            ccTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACg   <  1:2079268/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGG  >  CP009273/601703‑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: ‑