Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,969,299 A→G 6.1% Y265C (TAC→TGC)  wecE → TDP‑4‑oxo‑6‑deoxy‑D‑glucose transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,969,2990AG 6.1% 536.7 / 22.2 179Y265C (TAC→TGC) wecETDP‑4‑oxo‑6‑deoxy‑D‑glucose transaminase
Reads supporting (aligned to +/- strand):  ref base A (82/86);  new base G (7/4);  total (89/90)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-01

Only 100 of 184 total aligned reads displayed.

TTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGAAAGAAGCGGA  >  CP009273/3969154‑3969446
                                                                                                                                                 |                                                                                                                                                   
ttcCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACg                                                                                                                                                   <  2:1978053/148‑1 (MQ=255)
 tcCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGa                                                                                                                                                  <  2:2214585/148‑1 (MQ=255)
 tcCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGa                                                                                                                                                  <  2:2214018/148‑1 (MQ=255)
  ccGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGAt                                                                                                                                                 <  2:523350/148‑1 (MQ=255)
  ccGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGAt                                                                                                                                                 <  1:1881089/148‑1 (MQ=255)
     tGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCg                                                                                                                                              >  2:2963077/1‑148 (MQ=255)
         cAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAg                                                                                                                                          <  1:567380/148‑1 (MQ=255)
         cAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAg                                                                                                                                          >  1:2170669/1‑148 (MQ=255)
            gTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAgcgc                                                                                                                                       >  1:2402392/1‑148 (MQ=255)
            gTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAgcgc                                                                                                                                       >  2:149793/1‑148 (MQ=255)
             tctatatatatACGTGGCGCGATATTTGCTCCAGATATTTGATGTCCGATCTGCTAGCTGCATTCCTTTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTGTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGcc                                                                                                                                      <  2:3394803/145‑1 (MQ=255)
               gATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCctc                                                                                                                                    <  1:2873050/148‑1 (MQ=255)
                aTAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCctct                                                                                                                                   >  2:1665959/1‑148 (MQ=255)
                  aaaTATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTgg                                                                                                                                 <  2:2030975/148‑1 (MQ=255)
                     tataCCTTGCGCTATATTGGCTCCAGATATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTATGGCGCTGTGGCAAAACTAATACGATGCTTTAGCGCCTCTGGCGa                                                                                                                              <  2:1035128/148‑1 (MQ=255)
                     tataCCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCTTTAGCGCCTCTGGCGa                                                                                                                              <  2:1035448/148‑1 (MQ=255)
                      ataCCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGaa                                                                                                                             >  1:1769901/1‑148 (MQ=255)
                      ataCCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGaa                                                                                                                             >  1:2409508/1‑148 (MQ=255)
                      ataCCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGaa                                                                                                                             >  1:2782299/1‑148 (MQ=255)
                       taCCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGaaa                                                                                                                            >  2:1850321/1‑148 (MQ=255)
                        aCCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAg                                                                                                                           <  2:3500240/148‑1 (MQ=255)
                         ccTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGc                                                                                                                          <  1:2949209/148‑1 (MQ=255)
                          cTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGcc                                                                                                                         <  2:72566/148‑1 (MQ=255)
                            ggCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCgg                                                                                                                       <  2:3193889/148‑1 (MQ=255)
                                cgATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGt                                                                                                                   <  2:208075/148‑1 (MQ=255)
                                   tatTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTCCTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATc                                                                                                                >  2:2311972/1‑148 (MQ=255)
                                   tatTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATc                                                                                                                >  1:1097839/1‑148 (MQ=255)
                                   tatTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATc                                                                                                                >  1:3236848/1‑148 (MQ=255)
                                   tatTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATc                                                                                                                >  1:3216708/1‑148 (MQ=255)
                                    atTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTGCGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCg                                                                                                               >  1:432539/1‑148 (MQ=255)
                                                 aTTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGAt                                                                                                  <  1:2135711/148‑1 (MQ=255)
                                                 aTTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGAt                                                                                                  <  1:49311/148‑1 (MQ=255)
                                                  tttGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGAtt                                                                                                 <  1:2189266/148‑1 (MQ=255)
                                                  tttGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGAtt                                                                                                 >  2:3221821/1‑148 (MQ=255)
                                                   ttGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTc                                                                                                >  2:2496705/1‑148 (MQ=255)
                                                   ttGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTc                                                                                                <  1:2344839/148‑1 (MQ=255)
                                                   ttGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAAAGTCTGGAGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTc                                                                                                >  2:1791393/1‑148 (MQ=255)
                                                    tGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTcc                                                                                               <  1:1522511/148‑1 (MQ=255)
                                                        gTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGAt                                                                                           <  1:1572234/148‑1 (MQ=255)
                                                           cGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGc                                                                                        <  2:2345929/148‑1 (MQ=255)
                                                           cGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGc                                                                                        <  2:762515/148‑1 (MQ=255)
                                                           cGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGc                                                                                        <  1:3058550/148‑1 (MQ=255)
                                                             aTCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTGCGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTg                                                                                      <  1:1335057/148‑1 (MQ=255)
                                                               cTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCg                                                                                    <  1:383134/148‑1 (MQ=255)
                                                               cTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCg                                                                                    <  1:25914/148‑1 (MQ=255)
                                                                 gCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTg                                                                                  >  1:1309038/1‑148 (MQ=255)
                                                                 gCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTg                                                                                  >  2:1755864/1‑148 (MQ=255)
                                                                 gCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTg                                                                                  >  1:3083400/1‑148 (MQ=255)
                                                                 gCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTg                                                                                  >  2:1307070/1‑148 (MQ=255)
                                                                     gCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGa                                                                              >  1:1549961/1‑148 (MQ=255)
                                                                     gCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGa                                                                              >  2:95172/1‑148 (MQ=255)
                                                                      cTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTGCGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGaa                                                                             >  1:3314649/1‑148 (MQ=255)
                                                                        gCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAAcg                                                                           <  1:78572/148‑1 (MQ=255)
                                                                          aTACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAAcgcg                                                                         <  2:2670619/148‑1 (MQ=255)
                                                                            aCCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCa                                                                       >  2:2507544/1‑148 (MQ=255)
                                                                            aCCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCa                                                                       >  1:50904/1‑148 (MQ=255)
                                                                             ccTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTGCGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGATGCGTGCAGAACGCGCat                                                                      <  1:2027213/148‑1 (MQ=255)
                                                                             ccTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCat                                                                      >  2:3190287/1‑148 (MQ=255)
                                                                             ccTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCat                                                                      <  2:890160/148‑1 (MQ=255)
                                                                              cTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCata                                                                     >  1:3280287/1‑148 (MQ=255)
                                                                               tgtgGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCatat                                                                    <  2:3253583/148‑1 (MQ=255)
                                                                                gtgGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATg                                                                   <  2:2542632/148‑1 (MQ=255)
                                                                                 tgGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGt                                                                  <  2:590231/148‑1 (MQ=255)
                                                                                 tgGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGt                                                                  >  1:1685142/1‑148 (MQ=255)
                                                                                      gcAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTAc                                                             <  2:1511516/148‑1 (MQ=255)
                                                                                         aCTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACAtt                                                          >  2:680146/1‑148 (MQ=255)
                                                                                           tGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTaa                                                        <  1:168708/148‑1 (MQ=255)
                                                                                            ggAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCTGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTAAATTaaa                                                       >  2:110232/1‑148 (MQ=255)
                                                                                            ggAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTaaa                                                       >  1:1097462/1‑148 (MQ=255)
                                                                                            ggAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTaaa                                                       >  1:1533624/1‑148 (MQ=255)
                                                                                            ggAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTaaa                                                       <  1:1481400/148‑1 (MQ=255)
                                                                                            ggAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTaaa                                                       <  2:2391140/148‑1 (MQ=255)
                                                                                                gcagcGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGc                                                   >  1:266247/1‑148 (MQ=255)
                                                                                                  agcGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCgg                                                 >  1:786021/1‑148 (MQ=255)
                                                                                                  agcGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCgg                                                 <  2:2103395/148‑1 (MQ=255)
                                                                                                     ggATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGat                                              <  2:3015078/148‑1 (MQ=255)
                                                                                                      gATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTGCGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGata                                             <  2:2565649/148‑1 (MQ=255)
                                                                                                       aTCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGatat                                            >  1:423482/1‑148 (MQ=255)
                                                                                                       aTCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGatat                                            >  2:2773895/1‑148 (MQ=255)
                                                                                                          gTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATtga                                         <  2:2409447/148‑1 (MQ=255)
                                                                                                           tATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATtgat                                        >  2:1581476/1‑148 (MQ=255)
                                                                                                             tCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATtgatga                                      <  1:3067009/148‑1 (MQ=255)
                                                                                                               aaCCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGAcc                                    <  1:279854/148‑1 (MQ=255)
                                                                                                                aCCAGCAACGTCTGGCACTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCg                                   >  1:2501677/1‑148 (MQ=255)
                                                                                                                       aCGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGcgcg                            >  1:1178274/1‑148 (MQ=255)
                                                                                                                       aCGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGcgcg                            <  2:3527719/148‑1 (MQ=255)
                                                                                                                        cGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGt                           >  1:3468249/1‑148 (MQ=255)
                                                                                                                           cTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGa                        >  2:1104646/1‑148 (MQ=255)
                                                                                                                           cTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGa                        <  1:680146/148‑1 (MQ=255)
                                                                                                                            tGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGAt                       <  2:1547428/148‑1 (MQ=255)
                                                                                                                                 cTGTGGCAAAACTACTACGATGCGTGAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACt                  <  1:715994/148‑1 (MQ=255)
                                                                                                                                      gCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTg             >  2:2882014/1‑148 (MQ=255)
                                                                                                                                      gCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTg             >  2:1206837/1‑148 (MQ=255)
                                                                                                                                        aaaaCTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGaa           <  1:3162195/148‑1 (MQ=255)
                                                                                                                                          aaCTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGAaag         <  1:143167/148‑1 (MQ=255)
                                                                                                                                           actactGCGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGACCGCGCATAGGTTCTACATTAAACTGCGGGATCGTGAGGACCGGAGCGCGTTGATGAAGTGTCTGCaaga        >  2:2869733/1‑148 (MQ=255)
                                                                                                                                           actactACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGAaaga        >  2:2814491/1‑148 (MQ=255)
                                                                                                                                            ctactaCGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGAaagaa       <  1:1474760/148‑1 (MQ=255)
                                                                                                                                               ctacGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGAAAGAAGCg    >  1:1721642/1‑148 (MQ=255)
                                                                                                                                                 acGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGAAAGAAGCGGa  >  2:382828/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGAAAGAAGCGGA  >  CP009273/3969154‑3969446

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: