Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,120,556 A→T 16.2% L39Q (CTG→CAG)  wcaG ← bifunctional GDP‑fucose synthetase: GDP‑4‑dehydro‑6‑deoxy‑D‑mannose epimerase/ GDP‑4‑dehydro‑6‑L‑deoxygalactose reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,120,5560AT16.2% 447.5 / 69.5 165L39Q (CTG→CAG) wcaGbifunctional GDP‑fucose synthetase: GDP‑4‑dehydro‑6‑deoxy‑D‑mannose epimerase/ GDP‑4‑dehydro‑6‑L‑deoxygalactose reductase
Reads supporting (aligned to +/- strand):  ref base A (66/72);  new base T (19/8);  total (85/80)
Fisher's exact test for biased strand distribution p-value = 3.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 167 total aligned reads displayed.

CTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGACTCCAGCGCGATCGCCACGTCGT  >  CP009273/2120409‑2120702
                                                                                                                                                   |                                                                                                                                                  
ctctCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCa                                                                                                                                                    >  2:52641/1‑148 (MQ=255)
ctctCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCa                                                                                                                                                    >  2:52342/1‑148 (MQ=255)
   tcAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCt                                                                                                                                                 <  1:930836/148‑1 (MQ=255)
      atcatcATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCtcgt                                                                                                                                              <  2:584157/148‑1 (MQ=255)
      atcatcATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCtcgt                                                                                                                                              <  2:270822/148‑1 (MQ=255)
        catcatGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCtcgtcg                                                                                                                                            >  2:2875747/1‑148 (MQ=255)
           catGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCgg                                                                                                                                         >  1:787958/1‑148 (MQ=255)
            atGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCggg                                                                                                                                        >  1:2858870/1‑148 (MQ=255)
            atGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCggg                                                                                                                                        >  1:313744/1‑148 (MQ=255)
             tGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCTGCTCGTCGCGGGt                                                                                                                                       >  1:1329597/1‑148 (MQ=255)
              gTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTg                                                                                                                                      <  1:3007333/148‑1 (MQ=255)
                         tGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAg                                                                                                                           >  2:1944125/1‑148 (MQ=255)
                         tGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAg                                                                                                                           >  2:348129/1‑148 (MQ=255)
                            aaTCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTc                                                                                                                        >  2:2949412/1‑148 (MQ=255)
                            aaTCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTc                                                                                                                        <  1:2056181/148‑1 (MQ=255)
                            aaTCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTc                                                                                                                        >  2:708858/1‑148 (MQ=255)
                             aTCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTcc                                                                                                                       <  2:1007173/148‑1 (MQ=255)
                               ccgccgGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGACAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCcac                                                                                                                     >  1:2439468/1‑148 (MQ=255)
                                 gccgGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACAt                                                                                                                   <  2:672822/148‑1 (MQ=255)
                                 gccgGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACAt                                                                                                                   <  2:1061053/148‑1 (MQ=255)
                                 gccgGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTATGCTCGTCGCGGGTGCGTAATACCAGTTCCACAt                                                                                                                   <  2:252670/148‑1 (MQ=255)
                                     gATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCTGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCAcc                                                                                                               >  1:1078018/1‑148 (MQ=255)
                                      aTAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACcg                                                                                                              >  2:1321195/1‑148 (MQ=255)
                                      aTAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACcg                                                                                                              >  1:1784050/1‑148 (MQ=255)
                                       tAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCTGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACcgc                                                                                                             >  2:2236773/1‑148 (MQ=255)
                                       tAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACcgc                                                                                                             >  1:306670/1‑148 (MQ=255)
                                        aGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACcgcg                                                                                                            <  2:877024/148‑1 (MQ=255)
                                        aGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACcgcg                                                                                                            <  2:549563/148‑1 (MQ=255)
                                          gtgtTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCt                                                                                                          <  1:1563602/148‑1 (MQ=255)
                                          gtgtTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCt                                                                                                          >  2:1499244/1‑148 (MQ=255)
                                           tgttgtTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCTGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTg                                                                                                         <  2:1112559/148‑1 (MQ=255)
                                           tgttgtTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTg                                                                                                         <  1:1714660/148‑1 (MQ=255)
                                           tgttgtTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTg                                                                                                         <  1:1714817/148‑1 (MQ=255)
                                            gttgttGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGt                                                                                                        <  2:1304474/148‑1 (MQ=255)
                                            gttgttGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGt                                                                                                        <  1:687949/148‑1 (MQ=255)
                                              tgttgGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGGAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTc                                                                                                      >  1:2245263/1‑148 (MQ=255)
                                                 tgGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAg                                                                                                   <  2:1308872/148‑1 (MQ=255)
                                                          tGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCt                                                                                          >  1:1025905/1‑148 (MQ=255)
                                                          tGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCt                                                                                          >  1:1717354/1‑148 (MQ=255)
                                                               cccACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATgg                                                                                     <  1:3263362/148‑1 (MQ=255)
                                                                  aCTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCTGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCgg                                                                                  >  2:1908607/1‑148 (MQ=255)
                                                                   cTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGa                                                                                 <  1:1441609/148‑1 (MQ=255)
                                                                    tttCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGaa                                                                                <  2:3233735/148‑1 (MQ=255)
                                                                    tttCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGaa                                                                                <  2:1598867/148‑1 (MQ=255)
                                                                    tttCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGaa                                                                                <  2:1376113/148‑1 (MQ=255)
                                                                    tttCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGCTTCTGCTCGTCGCGGGTGCGTAAGACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGaa                                                                                <  2:2066887/148‑1 (MQ=255)
                                                                     ttCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCTGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAAc                                                                               >  2:792731/1‑148 (MQ=255)
                                                                     ttCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCTGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAAc                                                                               >  2:748391/1‑148 (MQ=255)
                                                                     ttCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCTGCGCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAAc                                                                               >  2:792738/1‑148 (MQ=255)
                                                                     ttCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAAc                                                                               >  2:586224/1‑148 (MQ=255)
                                                                         ccgccgccAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGAc                                                                           >  2:269289/1‑148 (MQ=255)
                                                                             cgccAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCTGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATc                                                                       <  2:321875/148‑1 (MQ=255)
                                                                             cgccAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCTGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATc                                                                       <  2:1411919/148‑1 (MQ=255)
                                                                                cAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCcg                                                                    <  1:2019190/148‑1 (MQ=255)
                                                                                  gATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCcgcg                                                                  <  2:3173587/148‑1 (MQ=255)
                                                                                      gACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGa                                                              >  2:479911/1‑148 (MQ=255)
                                                                                       aCCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGAc                                                             <  1:3176092/148‑1 (MQ=255)
                                                                                        ccTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCTGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGAcc                                                            <  1:1911165/148‑1 (MQ=255)
                                                                                            gTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCTGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCa                                                        >  1:2669291/1‑148 (MQ=255)
                                                                                            gTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCa                                                        >  2:1207672/1‑148 (MQ=255)
                                                                                            gTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCa                                                        <  2:2334276/148‑1 (MQ=255)
                                                                                              cAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAAt                                                      <  2:852669/148‑1 (MQ=255)
                                                                                              cAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAAt                                                      <  2:1486657/148‑1 (MQ=255)
                                                                                                aTACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATaa                                                    >  1:66255/1‑148 (MQ=255)
                                                                                                aTACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATaa                                                    >  2:2444461/1‑148 (MQ=255)
                                                                                                aTACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATaa                                                    <  2:352995/148‑1 (MQ=255)
                                                                                                aTACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATaa                                                    >  1:66217/1‑148 (MQ=255)
                                                                                                 tACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATaaa                                                   >  1:988709/1‑148 (MQ=255)
                                                                                                 tACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATaaa                                                   >  1:988650/1‑148 (MQ=255)
                                                                                                 tACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATaaa                                                   >  1:2746224/1‑148 (MQ=255)
                                                                                                 tACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATaaa                                                   >  2:1180883/1‑148 (MQ=255)
                                                                                                  aCGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATaaaa                                                  >  1:3053377/1‑148 (MQ=255)
                                                                                                       cGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCg                                             <  1:215874/148‑1 (MQ=255)
                                                                                                            gCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATACAAACTCgttgtt                                        <  2:1873721/148‑1 (MQ=255)
                                                                                                            gCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCgttgtt                                        <  2:2825301/148‑1 (MQ=255)
                                                                                                            gCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCgttgtt                                        <  2:2499112/148‑1 (MQ=255)
                                                                                                               aaGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTAc                                     >  1:1206419/1‑148 (MQ=255)
                                                                                                                aGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACt                                    >  2:2943041/1‑148 (MQ=255)
                                                                                                                      cATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCt                              <  1:269289/148‑1 (MQ=255)
                                                                                                                      cATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCt                              >  2:3308222/1‑148 (MQ=255)
                                                                                                                        tGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTa                            >  1:1506900/1‑148 (MQ=255)
                                                                                                                           aCGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGa                         <  2:609929/148‑1 (MQ=255)
                                                                                                                              gcgcGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGACTc                      <  1:816407/148‑1 (MQ=255)
                                                                                                                              gcgcGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGACTc                      <  1:2474100/148‑1 (MQ=255)
                                                                                                                              gcgcGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGACTc                      <  1:548087/148‑1 (MQ=255)
                                                                                                                              gcgcGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGACTc                      <  1:547817/148‑1 (MQ=255)
                                                                                                                               cgcgGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGACTcc                     >  2:1247797/1‑148 (MQ=255)
                                                                                                                                    cTGTCCAGCAGGTTCTGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGACTCCAgcgc                >  1:1670156/1‑148 (MQ=255)
                                                                                                                                    cTGTCCAGCAGGTTCTGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGACTCCAgcgc                >  2:222723/1‑148 (MQ=255)
                                                                                                                                    cTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGACTCCAgcgc                >  2:3007251/1‑148 (MQ=255)
                                                                                                                                      gTCCAGCAGGTTCTGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGACTCCAGCGCGa              <  1:2753337/148‑1 (MQ=255)
                                                                                                                                      gTCCAGCAGGTTCAGCTCGTCGCGGGTTCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCTCCTGATGGCGGAACCGACCATCCCGCGATGACCAGAAATAAAAACTCGTTGTTTACTCATGCTTATGACTCCAGCGCGa              <  1:1486522/148‑1 (MQ=255)
                                                                                                                                      gTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAATTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGACTCCAGCGCGa              <  2:1699853/148‑1 (MQ=255)
                                                                                                                                        ccAGCAGGTTCTGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGACTCCAGCGCGATc            >  1:250309/1‑148 (MQ=255)
                                                                                                                                        ccAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGACTCCAGCGCGATc            >  2:1457851/1‑148 (MQ=255)
                                                                                                                                          agcagGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGACTCCAGCGCGATCGc          >  2:1768396/1‑148 (MQ=255)
                                                                                                                                          agcagGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGACTCCAGCGCGATCGc          <  2:1807161/148‑1 (MQ=255)
                                                                                                                                                ttCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGACTCCAGCGCGATCGCCAcgtc    <  2:1025905/148‑1 (MQ=255)
                                                                                                                                                 tCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGACTCCAGCGCGATCGCCAcgtcg   <  2:2892442/148‑1 (MQ=255)
                                                                                                                                                  cAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGACTCCAGCGCGATCGCCAcgtcgt  <  1:2949412/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCATGCACGGCGCGGCTGTCCAGCAGGTTCAGCTCGTCGCGGGTGCGTAATACCAGTTCCACATCACCGCGCTGTTCGAGCTGCCGCCTGATGGCGGAACCGACCATCCCGCGATGACCAGCAATAAAAACTCGTTGTTTACTCATGCTTATGACTCCAGCGCGATCGCCACGTCGT  >  CP009273/2120409‑2120702

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: