Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,245,682 T→G 20.1% T205P (ACC→CCC)  dhaK ← dihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,245,6820TG20.1% 607.0 / 200.7 299T205P (ACC→CCC) dhaKdihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit
Reads supporting (aligned to +/- strand):  ref base T (150/89);  new base G (40/20);  total (190/109)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

Only 100 of 306 total aligned reads displayed.

GAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTA  >  CP009273/1245533‑1245831
                                                                                                                                                     |                                                                                                                                                     
  aaCGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGtca                                                                                                                                                     >  2:1523982/1‑148 (MQ=255)
      cAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAaaga                                                                                                                                                 >  1:2093468/1‑150 (MQ=255)
      cAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAaaga                                                                                                                                                 >  1:1303215/1‑150 (MQ=255)
      cAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAaaga                                                                                                                                                 <  1:1881539/150‑1 (MQ=255)
       aaaGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAaagaa                                                                                                                                                >  1:1208257/1‑150 (MQ=255)
       aaaGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAaagaa                                                                                                                                                >  1:3619434/1‑150 (MQ=255)
        aaGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAaagaag                                                                                                                                               >  1:2326602/1‑150 (MQ=255)
        aaGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAaagaag                                                                                                                                               >  1:4409945/1‑150 (MQ=255)
        aaGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAaagaag                                                                                                                                               <  2:944623/150‑1 (MQ=255)
         aGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAgg                                                                                                                                              <  1:4218344/150‑1 (MQ=255)
           tGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGtt                                                                                                                                            <  1:895057/150‑1 (MQ=255)
             cGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTg                                                                                                                                          >  2:1002704/1‑150 (MQ=255)
              gATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGc                                                                                                                                         >  1:4195998/1‑150 (MQ=255)
               aTGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGcc                                                                                                                                        >  2:1571746/1‑150 (MQ=255)
                 ggTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCcg                                                                                                                                      >  2:1107166/1‑150 (MQ=255)
                    aTGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCgg                                                                                                                                   >  1:501293/1‑150 (MQ=255)
                      gAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCa                                                                                                                                 >  1:218936/1‑150 (MQ=255)
                      gAGCCATTTACCAGCAGGGAGTCGAACATTTCATCGACCGTTAGATCAAGGGAAGAGAACGGGCGGCGGTCAATACCCGGCACACCATGAATGCCGACGCCAAACTCCATCTCATTATACGCCAGGGTAAAAGAAGGTTTGCACACGGCa                                                                                                                                 >  1:4538182/1‑150 (MQ=255)
                       aGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAg                                                                                                                                <  1:2141021/150‑1 (MQ=255)
                       aGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAg                                                                                                                                >  2:2837148/1‑150 (MQ=255)
                       aGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAg                                                                                                                                <  2:1989751/150‑1 (MQ=255)
                        gCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAgg                                                                                                                               >  1:3506592/1‑150 (MQ=255)
                         ccATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGa                                                                                                                              <  2:3431666/150‑1 (MQ=255)
                         ccATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGa                                                                                                                              >  2:1399406/1‑150 (MQ=255)
                         ccATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGa                                                                                                                              >  2:2114772/1‑150 (MQ=255)
                         ccATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGa                                                                                                                              >  1:577972/1‑150 (MQ=255)
                           aTTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGTGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGAAc                                                                                                                            <  2:3698266/150‑1 (MQ=255)
                             ttACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGTGCGGCGGTCAATAACCGGCTCACCATGAATGCAGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGAACgg                                                                                                                          >  2:4741494/1‑150 (MQ=255)
                               aCCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGATGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTa                                                                                                                        <  2:1303215/150‑1 (MQ=255)
                               aCCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTa                                                                                                                        >  2:728913/1‑150 (MQ=255)
                                  agcagGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAg                                                                                                                     <  2:665539/150‑1 (MQ=255)
                                    cagGGTGTCGAACATTTAATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATATCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGc                                                                                                                   >  1:2177139/1‑150 (MQ=255)
                                         tgtCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCga                                                                                                              <  2:4404664/150‑1 (MQ=255)
                                               aCATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGa                                                                                                        >  2:3116225/1‑150 (MQ=255)
                                                cATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGAt                                                                                                       >  2:338592/1‑150 (MQ=255)
                                                  tttCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATAc                                                                                                     <  2:3786855/150‑1 (MQ=255)
                                                   ttCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATAcc                                                                                                    >  1:796996/1‑150 (MQ=255)
                                                       tCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTAtt                                                                                                <  2:2792983/150‑1 (MQ=255)
                                                        cGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTg                                                                                               <  2:2972629/150‑1 (MQ=255)
                                                        cGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTg                                                                                               >  2:406347/1‑150 (MQ=255)
                                                        cGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTg                                                                                               >  2:539577/1‑150 (MQ=255)
                                                         gACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGa                                                                                              >  1:3521400/1‑150 (MQ=255)
                                                         gACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGa                                                                                              >  1:1078196/1‑150 (MQ=255)
                                                               ttGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCGATTGAGTGGcc                                                                                        >  1:3858897/1‑150 (MQ=255)
                                                                 gATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCtt                                                                                      >  2:4624110/1‑150 (MQ=255)
                                                                   tCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGa                                                                                    <  2:1051067/150‑1 (MQ=255)
                                                                   tCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGa                                                                                    >  1:1018752/1‑150 (MQ=255)
                                                                    cAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGat                                                                                   <  2:4483562/150‑1 (MQ=255)
                                                                    cAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGat                                                                                   >  2:2601017/1‑150 (MQ=255)
                                                                     aaGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGatt                                                                                  >  2:2469378/1‑150 (MQ=255)
                                                                     aaGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGatt                                                                                  <  1:3362146/150‑1 (MQ=255)
                                                                     aaGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGatt                                                                                  <  1:3962486/150‑1 (MQ=255)
                                                                         gAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTc                                                                              >  2:3528147/1‑150 (MQ=255)
                                                                          aaGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCa                                                                             <  2:4355409/150‑1 (MQ=255)
                                                                          aaGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCa                                                                             >  1:3602154/1‑150 (MQ=255)
                                                                              gaAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCtt                                                                         >  1:4623258/1‑150 (MQ=255)
                                                                               aaGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTa                                                                        <  2:2996722/150‑1 (MQ=255)
                                                                               aaGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTa                                                                        >  2:1425758/1‑150 (MQ=255)
                                                                                 ggggCGGCGGTCAAGACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACg                                                                      <  1:1778277/150‑1 (MQ=255)
                                                                                     cggcggTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGcccc                                                                  <  1:4485899/150‑1 (MQ=255)
                                                                                     cggcggTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGcccc                                                                  <  1:2634734/150‑1 (MQ=255)
                                                                                        cggTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGt                                                               >  2:3049647/1‑150 (MQ=255)
                                                                                           tCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTcc                                                            >  1:4617572/1‑150 (MQ=255)
                                                                                           tCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTcc                                                            >  1:1322812/1‑150 (MQ=255)
                                                                                           tCAATACCCGGCTCACCATGAATGCCGAAGCCAAACACCATCTCATTATCCGCAAGGGTAAAAGAAGGTTTGCACGCGGCAGGAACTGTACAGGCACCGAGAGCGATACCTCTTTAGTGGCCTTGATTATTCAGCTTACGCCACAGATcc                                                            >  2:4207901/1‑150 (MQ=255)
                                                                                             aaTACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGc                                                          >  1:2058680/1‑150 (MQ=255)
                                                                                             aaTACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGc                                                          >  1:691554/1‑150 (MQ=255)
                                                                                                   cGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGc                                                    >  1:26792/1‑150 (MQ=255)
                                                                                                   cGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGc                                                    >  2:1010702/1‑150 (MQ=255)
                                                                                                   cGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGc                                                    >  1:2856087/1‑150 (MQ=255)
                                                                                                   cGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGc                                                    >  1:2938728/1‑150 (MQ=255)
                                                                                                    ggCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCg                                                   >  2:3736371/1‑150 (MQ=255)
                                                                                                      cTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTc                                                 >  2:3130907/1‑150 (MQ=255)
                                                                                                      cTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTc                                                 <  2:1208257/150‑1 (MQ=255)
                                                                                                       tCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTcc                                                >  1:3445702/1‑150 (MQ=255)
                                                                                                           cATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTg                                            <  1:3081764/150‑1 (MQ=255)
                                                                                                             tGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAg                                          >  1:823853/1‑150 (MQ=255)
                                                                                                                   ccGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCa                                    >  2:886516/1‑150 (MQ=255)
                                                                                                                    cGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGGTCCGCACAGGCGTCCAGTGAGTCGCCAc                                   <  2:4311038/150‑1 (MQ=255)
                                                                                                                       cGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCt                                <  1:2545245/150‑1 (MQ=255)
                                                                                                                        gCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTc                               >  2:3608656/1‑150 (MQ=255)
                                                                                                                        gCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTc                               >  1:4569749/1‑150 (MQ=255)
                                                                                                                          cAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTccg                             <  1:1714599/150‑1 (MQ=255)
                                                                                                                            aaCTCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTccgcc                           >  1:2638109/1‑150 (MQ=255)
                                                                                                                               tCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCt                        >  2:4002640/1‑150 (MQ=255)
                                                                                                                               tCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCt                        <  1:338592/150‑1 (MQ=255)
                                                                                                                               tCCATCTCATTATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCt                        <  1:3745454/150‑1 (MQ=255)
                                                                                                                                  aTCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTgcg                     >  2:643723/1‑150 (MQ=255)
                                                                                                                                    ctcATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGcc                   <  1:1576862/150‑1 (MQ=255)
                                                                                                                                        ttATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACg               <  1:1992105/150‑1 (MQ=255)
                                                                                                                                        ttATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACg               >  1:2986464/1‑150 (MQ=255)
                                                                                                                                        ttATCCGCCAGGGGAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACg               <  2:1246150/150‑1 (MQ=255)
                                                                                                                                          aTCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAg             <  2:4195998/150‑1 (MQ=255)
                                                                                                                                              gCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGtttt         >  1:1060231/1‑150 (MQ=255)
                                                                                                                                               ccAGGGGAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGttttt        <  1:32213/150‑1 (MQ=255)
                                                                                                                                                cAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTc       >  2:2148082/1‑150 (MQ=255)
                                                                                                                                                 aGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCa      >  1:4373196/1‑150 (MQ=255)
                                                                                                                                                    gTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAAtt   <  2:2326602/150‑1 (MQ=255)
                                                                                                                                                     tAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTa  >  2:391517/1‑150 (MQ=255)
                                                                                                                                                     tAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTa  >  2:951221/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTA  >  CP009273/1245533‑1245831

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: