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% 749.6 / 9.1 232I331I (ATT→ATA) ybdKweak gamma‑glutamyl:cysteine ligase
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base T (104/121);  total (107/125)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 232 total aligned reads displayed.

TAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACG  >  CP009273/601705‑601990
                                                                                                                                              |                                                                                                                                               
tAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATcgc                                                                                                                                            >  1:1924131/1‑148 (MQ=255)
 aaGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATcgcg                                                                                                                                           <  2:2654284/148‑1 (MQ=255)
  aGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCgcgc                                                                                                                                          <  1:309962/148‑1 (MQ=255)
  aGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCgcgc                                                                                                                                          <  2:196449/148‑1 (MQ=255)
  aGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCgcgc                                                                                                                                          >  2:2864371/1‑148 (MQ=255)
   ggCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCt                                                                                                                                         <  1:2672038/148‑1 (MQ=255)
   ggCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCt                                                                                                                                         <  1:199015/148‑1 (MQ=255)
     cAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCg                                                                                                                                       >  1:3436820/1‑148 (MQ=255)
     cAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCg                                                                                                                                       <  1:3090706/148‑1 (MQ=255)
        ttCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATg                                                                                                                                    <  2:3442634/148‑1 (MQ=255)
         tCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGc                                                                                                                                   >  2:124753/1‑148 (MQ=255)
          cccGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCa                                                                                                                                  >  1:1354522/1‑148 (MQ=255)
           ccGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCAc                                                                                                                                 <  1:751137/148‑1 (MQ=255)
           ccGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCAc                                                                                                                                 >  2:636083/1‑148 (MQ=255)
             ggCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCa                                                                                                                               >  1:775638/1‑148 (MQ=255)
             ggCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCa                                                                                                                               <  1:3159865/148‑1 (MQ=255)
                  aGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAAtttt                                                                                                                          <  1:3017648/148‑1 (MQ=255)
                   gTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTa                                                                                                                         >  2:1086287/1‑148 (MQ=255)
                   gTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTa                                                                                                                         <  2:1162134/148‑1 (MQ=255)
                   gTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTa                                                                                                                         <  2:2837760/148‑1 (MQ=255)
                     cACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAtg                                                                                                                       <  1:97506/148‑1 (MQ=255)
                     cACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAtg                                                                                                                       >  1:3374219/1‑148 (MQ=255)
                       ccGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTAtgtg                                                                                                                     <  2:223148/148‑1 (MQ=255)
                        cGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGc                                                                                                                    >  2:1036157/1‑148 (MQ=255)
                               gATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAggg                                                                                                             <  2:3514463/148‑1 (MQ=255)
                               gATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAggg                                                                                                             >  2:2056738/1‑148 (MQ=255)
                               gATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAggg                                                                                                             <  2:1111564/148‑1 (MQ=255)
                               gATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAggg                                                                                                             <  2:1895945/148‑1 (MQ=255)
                               gATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAggg                                                                                                             <  2:1895401/148‑1 (MQ=255)
                                aTCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAgggg                                                                                                            >  2:3008187/1‑148 (MQ=255)
                                 tctcACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGc                                                                                                           <  2:548984/148‑1 (MQ=255)
                                      caATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAttt                                                                                                      <  1:639405/148‑1 (MQ=255)
                                         tGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGAtttttt                                                                                                   >  1:1218135/1‑148 (MQ=255)
                                           cTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTcc                                                                                                 <  1:3034131/148‑1 (MQ=255)
                                                tCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCaa                                                                                            <  2:918420/148‑1 (MQ=255)
                                                 cACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAAt                                                                                           <  1:1617456/148‑1 (MQ=255)
                                                    cAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGc                                                                                        <  1:2626327/148‑1 (MQ=255)
                                                      gCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCaa                                                                                      >  2:1283927/1‑148 (MQ=255)
                                                          aaTCAGCGAGCCGCCATCGGCGACGAAATAGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTa                                                                                  >  1:1191289/1‑148 (MQ=255)
                                                           aTCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTAt                                                                                 >  1:2613182/1‑148 (MQ=255)
                                                            tCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATc                                                                                >  1:356101/1‑148 (MQ=255)
                                                            tCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATc                                                                                <  1:389024/148‑1 (MQ=255)
                                                              aGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCtt                                                                              <  1:2472659/148‑1 (MQ=255)
                                                               gCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTc                                                                             <  2:995059/148‑1 (MQ=255)
                                                                cGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTcc                                                                            >  1:3274563/1‑148 (MQ=255)
                                                                cGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTcc                                                                            <  1:1148982/148‑1 (MQ=255)
                                                                 gAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCg                                                                           <  2:2892766/148‑1 (MQ=255)
                                                                   gccgccATCGGCGACGAAATCGCGAATTAGCTGCGCTTCATTCATACCGCTGACGAACTGGCTATGCAGGGCCTCTATCTAGCTCGATGAACCAATTTTATGTGCGGAAGGTGCGATTTTTTCCAGCAAACGCAAGGTATCTTCCgat                                                                         >  2:3192517/1‑146 (MQ=255)
                                                                    ccgccATCGGCGACGAAATCGCGCATGAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTa                                                                        <  2:1805809/148‑1 (MQ=255)
                                                                       ccATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCg                                                                     >  2:439840/1‑148 (MQ=255)
                                                                       ccATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCg                                                                     >  2:1876590/1‑148 (MQ=255)
                                                                        cATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCgg                                                                    >  2:2838596/1‑148 (MQ=255)
                                                                            ggCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTcga                                                                >  1:724667/1‑148 (MQ=255)
                                                                              cgacgaAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTcgacg                                                              <  2:2694747/148‑1 (MQ=255)
                                                                                acgaAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGAt                                                            <  1:496555/148‑1 (MQ=255)
                                                                                acgaAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGAt                                                            >  2:214065/1‑148 (MQ=255)
                                                                                 cgaAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGAtc                                                           <  2:2668602/148‑1 (MQ=255)
                                                                                   aaaTCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGAtctc                                                         <  1:1619852/148‑1 (MQ=255)
                                                                                   aaaTCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGAtctc                                                         >  1:770064/1‑148 (MQ=255)
                                                                                    aaTCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTcc                                                        <  1:3008187/148‑1 (MQ=255)
                                                                                     aTCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCa                                                       <  1:1002541/148‑1 (MQ=255)
                                                                                      tCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAg                                                      >  2:3024362/1‑148 (MQ=255)
                                                                                       cgcgCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAgt                                                     >  2:821845/1‑148 (MQ=255)
                                                                                        gcgcATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAgtg                                                    <  2:1763806/148‑1 (MQ=255)
                                                                                          gcATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGtgtg                                                  <  1:2864371/148‑1 (MQ=255)
                                                                                             ttAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCACGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCgg                                               >  2:2388682/1‑148 (MQ=255)
                                                                                             ttAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCgg                                               >  1:631909/1‑148 (MQ=255)
                                                                                              tAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCTCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGa                                              <  2:2319245/148‑1 (MQ=255)
                                                                                                 cTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCg                                           >  2:958602/1‑148 (MQ=255)
                                                                                                 cTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCg                                           >  1:3440594/1‑148 (MQ=255)
                                                                                                 cTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCAGTTAGCGGTCGACGATCTCCAGTGTGCGGATCg                                           >  2:1602620/1‑148 (MQ=255)
                                                                                                   gcgcTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGt                                         <  2:2420595/148‑1 (MQ=255)
                                                                                                   gcgcTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGt                                         <  1:1034004/148‑1 (MQ=255)
                                                                                                      cTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtgat                                      >  2:525801/1‑148 (MQ=255)
                                                                                                      cTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGtgat                                      >  2:2289881/1‑148 (MQ=255)
                                                                                                       ttcattcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGTGGAGATTTTTTCCAGCAATCGCAATGTATCTTGCGTTAGCGGTCGACTATCTCCAGTGTGCGGATCGGtgatg                                     <  1:188151/148‑1 (MQ=255)
                                                                                                            tcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGcc                                <  1:1082117/148‑1 (MQ=255)
                                                                                                            tcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGcc                                >  2:1233643/1‑148 (MQ=255)
                                                                                                            tcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGcc                                <  2:394108/148‑1 (MQ=255)
                                                                                                            tcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGcc                                <  2:150201/148‑1 (MQ=255)
                                                                                                            tcaGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGcc                                >  1:540930/1‑148 (MQ=255)
                                                                                                             caGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCt                               >  1:3241139/1‑148 (MQ=255)
                                                                                                              aGACCGCTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCtt                              <  1:3052057/148‑1 (MQ=255)
                                                                                                                    cTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGcc                        >  2:2117271/1‑148 (MQ=255)
                                                                                                                    cTGACGACCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGcc                        <  2:1753124/148‑1 (MQ=255)
                                                                                                                       acgacCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCAt                     <  2:775349/148‑1 (MQ=255)
                                                                                                                         gacCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAg                   <  2:1256673/148‑1 (MQ=255)
                                                                                                                         gacCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAg                   <  1:2437857/148‑1 (MQ=255)
                                                                                                                         gacCTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAg                   <  1:2804337/148‑1 (MQ=255)
                                                                                                                            cTGGCGATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGa                >  2:3395779/1‑148 (MQ=255)
                                                                                                                                 gATGCAGGGCCTCTCTCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGc           <  2:208823/148‑1 (MQ=255)
                                                                                                                                 gATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGc           <  2:76610/148‑1 (MQ=255)
                                                                                                                                 gATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGc           <  1:883338/148‑1 (MQ=255)
                                                                                                                                 gATGCAGGGCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGc           <  2:2444100/148‑1 (MQ=255)
                                                                                                                                       gggCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGaa     >  1:2525541/1‑148 (MQ=255)
                                                                                                                                        ggCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGaaa    <  1:1284487/148‑1 (MQ=255)
                                                                                                                                        ggCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGaaa    <  1:439840/148‑1 (MQ=255)
                                                                                                                                        ggCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGaaa    >  1:2886103/1‑148 (MQ=255)
                                                                                                                                        ggCCTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGaaa    <  1:3095306/148‑1 (MQ=255)
                                                                                                                                          ccTCTATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACg  <  2:2512066/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
TAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCGCTCGATGCACCAATTTTATGTGCGGAAGGGGCGATTTTTTCCAGCAATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACG  >  CP009273/601705‑601990

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