Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,898,349 G→A 5.1% F203F (TTC→TTT bglF ← fused beta‑glucoside‑specific PTS enzymes: IIA component/IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,898,3490GA 5.1% 1270.0 / 37.4 336F203F (TTC→TTTbglFfused beta‑glucoside‑specific PTS enzymes: IIA component/IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base G (182/137);  new base A (7/10);  total (189/147)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 337 total aligned reads displayed.

TCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGA  >  CP009273/3898201‑3898497
                                                                                                                                                    |                                                                                                                                                    
 cAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGaa                                                                                                                                                    >  2:1624743/1‑150 (MQ=255)
 cAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGaa                                                                                                                                                    <  1:3523365/150‑1 (MQ=255)
   gACAGAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAAt                                                                                                                                                  <  2:3511497/150‑1 (MQ=255)
     cATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATcc                                                                                                                                                <  2:1739614/150‑1 (MQ=255)
      aTAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCa                                                                                                                                               >  2:1311736/1‑150 (MQ=255)
        aGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGc                                                                                                                                             >  1:568596/1‑150 (MQ=255)
        aGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGc                                                                                                                                             >  1:1466925/1‑150 (MQ=255)
           aaTGGTGTGAAGAAATTTGTGATGGCCGACGGTGACCACGCATTAAGTCGGCGGTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAAGTCAACAATGTGACCGGAATACCCAGGAAATCCAGcccc                                                                                                                                          <  2:1026466/150‑1 (MQ=255)
              ggTGTGAAGAAATTTTTGCTTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAgc                                                                                                                                       <  2:3892224/150‑1 (MQ=255)
                 gtgAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCa                                                                                                                                    >  1:3807289/1‑150 (MQ=255)
                    aagaaATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGAAAATCCAGCCCCAGCGCATcc                                                                                                                                 <  2:524734/150‑1 (MQ=255)
                     agaaATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCg                                                                                                                                >  2:3604518/1‑150 (MQ=255)
                       aaaTTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCt                                                                                                                              <  1:4487014/150‑1 (MQ=255)
                       aaaTTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCt                                                                                                                              <  1:4075373/150‑1 (MQ=255)
                        aaTTTTTGATTGCCGACGGTAACAACGCATTAAGTCTGCGTTACAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCtt                                                                                                                             <  1:3344505/150‑1 (MQ=255)
                         aTTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCttt                                                                                                                            <  2:4549426/150‑1 (MQ=255)
                         aTTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCttt                                                                                                                            <  1:1504142/150‑1 (MQ=255)
                            tttGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGAAAATCCAGCCCCAGCGCATCCGCTTTTTg                                                                                                                         <  1:2232622/150‑1 (MQ=255)
                              tGATTGCCGACGGTAAACACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACTAGTAATTCAACAATTTGACCGGAATACACAGGAAATCCAGCCCCAGCGCATCCGCTATTTGcc                                                                                                                       >  1:9554/1‑150 (MQ=255)
                               gATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGccc                                                                                                                      >  2:724850/1‑150 (MQ=255)
                                aTTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCg                                                                                                                     >  2:4299675/1‑150 (MQ=255)
                                   gCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGAAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTtc                                                                                                                  <  2:3250779/150‑1 (MQ=255)
                                     cGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTtctc                                                                                                                >  1:1538248/1‑150 (MQ=255)
                                      gACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCg                                                                                                               >  1:89893/1‑150 (MQ=255)
                                       aCGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGAAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGa                                                                                                              >  2:399894/1‑150 (MQ=255)
                                        cGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGaa                                                                                                             >  1:4190293/1‑150 (MQ=255)
                                        cGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGaa                                                                                                             >  2:1211046/1‑150 (MQ=255)
                                          gTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAg                                                                                                           <  2:212255/150‑1 (MQ=255)
                                           tAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGc                                                                                                          >  1:2605784/1‑150 (MQ=255)
                                           tAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGc                                                                                                          >  1:2164893/1‑150 (MQ=255)
                                             aCCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAg                                                                                                        >  1:5001250/1‑150 (MQ=255)
                                                aCGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCa                                                                                                     >  1:3586490/1‑150 (MQ=255)
                                                 cGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAg                                                                                                    >  1:4859787/1‑150 (MQ=255)
                                                  gCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGa                                                                                                   <  1:316946/150‑1 (MQ=255)
                                                  gCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGa                                                                                                   <  1:511561/150‑1 (MQ=255)
                                                   cATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGaa                                                                                                  <  1:1221845/150‑1 (MQ=255)
                                                        aGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAAt                                                                                             <  1:1920782/150‑1 (MQ=255)
                                                         gTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATg                                                                                            >  2:2215084/1‑150 (MQ=255)
                                                          tCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATgg                                                                                           <  2:880447/150‑1 (MQ=255)
                                                           cGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGa                                                                                          >  1:1185898/1‑150 (MQ=255)
                                                            ggCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGAt                                                                                         >  1:2423210/1‑150 (MQ=255)
                                                            ggCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGAt                                                                                         >  1:2095913/1‑150 (MQ=255)
                                                            ggCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGAt                                                                                         >  2:4313964/1‑150 (MQ=255)
                                                               gTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCa                                                                                      >  2:4611387/1‑150 (MQ=255)
                                                                 tCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACt                                                                                    >  2:4894620/1‑150 (MQ=255)
                                                                   cAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTTATGGATGCACTaa                                                                                  <  1:759380/150‑1 (MQ=255)
                                                                   cAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTaa                                                                                  >  2:3785465/1‑150 (MQ=255)
                                                                    aGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAAc                                                                                 >  2:1082310/1‑150 (MQ=255)
                                                                      aaTGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGc                                                                               >  2:3162419/1‑150 (MQ=255)
                                                                       aTGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCt                                                                              <  2:2804087/150‑1 (MQ=255)
                                                                       aTGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCt                                                                              <  2:537986/150‑1 (MQ=255)
                                                                        tgctgcACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTc                                                                             >  1:2722183/1‑150 (MQ=255)
                                                                           tgcACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTccac                                                                          >  2:440266/1‑150 (MQ=255)
                                                                             cacaACCAGGCAGAAAAAATATTGTGAATAACCGATGACGAGTAATTCATCAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGATAGCAGTCAGAATTAATGGATGCACTAACGCTccacca                                                                        <  1:218210/150‑1 (MQ=255)
                                                                             cacaACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTccacca                                                                        <  1:4393910/150‑1 (MQ=255)
                                                                                aaCCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAAta                                                                     >  1:474522/1‑149 (MQ=255)
                                                                                      gCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGAAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATg                                                               >  2:1890338/1‑150 (MQ=255)
                                                                                        aGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGc                                                             >  1:248275/1‑150 (MQ=255)
                                                                                         gAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCa                                                            >  1:1025595/1‑150 (MQ=255)
                                                                                         gAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCa                                                            <  2:2038646/150‑1 (MQ=255)
                                                                                           aaaaaTAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGt                                                          <  1:710177/150‑1 (MQ=255)
                                                                                           aaaaaTAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGt                                                          >  2:3844160/1‑150 (MQ=255)
                                                                                             aaaTAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTaa                                                        >  2:2423729/1‑150 (MQ=255)
                                                                                             aaaTAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTaa                                                        >  1:2900535/1‑150 (MQ=255)
                                                                                              aataatGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTaaa                                                       >  1:3724705/1‑150 (MQ=255)
                                                                                              aataatGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGAAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTaaa                                                       <  1:3689793/150‑1 (MQ=255)
                                                                                                 aatGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATgg                                                    <  1:1836904/150‑1 (MQ=255)
                                                                                                    gggAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGAtt                                                 <  2:2956281/150‑1 (MQ=255)
                                                                                                    gggAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGAtt                                                 >  2:3188607/1‑150 (MQ=255)
                                                                                                        aTAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTAccg                                             <  2:1491507/150‑1 (MQ=255)
                                                                                                        aTAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTAccg                                             >  2:1761232/1‑150 (MQ=255)
                                                                                                        aTAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGAAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTAccg                                             <  1:4683660/150‑1 (MQ=255)
                                                                                                          aaCCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTAccgcc                                           <  1:4963329/150‑1 (MQ=255)
                                                                                                            ccGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGa                                         >  1:2961728/1‑150 (MQ=255)
                                                                                                             cGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGaa                                        >  1:3878771/1‑150 (MQ=255)
                                                                                                             cGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGaa                                        >  2:3899550/1‑150 (MQ=255)
                                                                                                             cGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCAGGGCAGTAAATGGATTACCGCCGaa                                        <  2:4958267/150‑1 (MQ=255)
                                                                                                               aTGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGc                                      >  2:2352698/1‑150 (MQ=255)
                                                                                                                 gACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGt                                    >  2:1972849/1‑150 (MQ=255)
                                                                                                                 gACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGt                                    <  1:4406169/150‑1 (MQ=255)
                                                                                                                   cGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGttt                                  >  2:1397258/1‑150 (MQ=255)
                                                                                                                    gAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTc                                 >  1:4791124/1‑150 (MQ=255)
                                                                                                                       tAATTCAACAATGTGACCGGAATACCCAGAAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTcccc                              <  2:4094454/150‑1 (MQ=255)
                                                                                                                        aaTTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCcg                             <  2:3391874/150‑1 (MQ=255)
                                                                                                                            caacaaTGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGgt                         >  2:1795171/1‑150 (MQ=255)
                                                                                                                               caaTGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTa                      >  1:339720/1‑150 (MQ=255)
                                                                                                                               caaTGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTa                      >  1:167206/1‑150 (MQ=255)
                                                                                                                               caaTGTGACCGGAATACCCAGAAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTa                      <  1:4363670/150‑1 (MQ=255)
                                                                                                                                 aTGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATc                    <  1:1624743/150‑1 (MQ=255)
                                                                                                                                  tgtgACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATcc                   <  2:3050675/150‑1 (MQ=255)
                                                                                                                                   gtgACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATccc                  <  2:1216022/150‑1 (MQ=255)
                                                                                                                                     gACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAg                >  2:4788137/1‑150 (MQ=255)
                                                                                                                                       ccGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGa              <  2:2635779/150‑1 (MQ=255)
                                                                                                                                          gAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTa           <  1:3242774/150‑1 (MQ=255)
                                                                                                                                          gAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTa           <  1:4551645/150‑1 (MQ=255)
                                                                                                                                             tACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCg        >  1:978172/1‑150 (MQ=255)
                                                                                                                                              aCCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCgg       <  1:824087/150‑1 (MQ=255)
                                                                                                                                               cccAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTCACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCggg      <  1:406112/150‑1 (MQ=255)
                                                                                                                                               cccAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCggg      >  2:1014153/1‑150 (MQ=255)
                                                                                                                                                 cAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGgaa    <  2:5100838/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGA  >  CP009273/3898201‑3898497

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: