Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,767,242 G→T 6.4% A534A (GCG→GCT mtlA → fused mannitol‑specific PTS enzymes: IIA components/IIB components/IIC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,767,2420GT 6.4% 1160.4 / 31.0 315A534A (GCG→GCTmtlAfused mannitol‑specific PTS enzymes: IIA components/IIB components/IIC components
Reads supporting (aligned to +/- strand):  ref base G (160/134);  new base T (7/13);  total (167/148)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 321 total aligned reads displayed.

GCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGA  >  CP009273/3767094‑3767391
                                                                                                                                                    |                                                                                                                                                     
  tttGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATg                                                                                                                                                    <  1:2075372/150‑1 (MQ=255)
    tGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCt                                                                                                                                                  <  2:1222908/150‑1 (MQ=255)
      aCGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTgg                                                                                                                                                >  2:496099/1‑150 (MQ=255)
        gATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGAt                                                                                                                                              >  1:5027649/1‑150 (MQ=255)
        gATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGAt                                                                                                                                              >  2:5081526/1‑150 (MQ=255)
        gATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGAt                                                                                                                                              >  1:1993693/1‑150 (MQ=255)
           tCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGt                                                                                                                                           >  2:225375/1‑150 (MQ=255)
            ccAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCTATGCTGGATCGTg                                                                                                                                          <  2:2845649/150‑1 (MQ=255)
             cAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGa                                                                                                                                         >  2:1806421/1‑150 (MQ=255)
              aGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGaa                                                                                                                                        >  1:2244916/1‑150 (MQ=255)
                tGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGACGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGaaaa                                                                                                                                      <  1:5086307/150‑1 (MQ=255)
                tGCTAACCTGTTAAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAATAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGaaaa                                                                                                                                      <  1:2335563/150‑1 (MQ=255)
                   tAACCTGTTCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACt                                                                                                                                   <  1:2989606/150‑1 (MQ=255)
                    aaCCTGTTCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTg                                                                                                                                  <  1:803666/150‑1 (MQ=255)
                    aaCCTGTTCAAGCTAGGCGCGGAGAACATCTTCCACGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTTCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAATCTg                                                                                                                                  >  1:1785469/1‑150 (MQ=255)
                     aCCTGTTCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAACGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGa                                                                                                                                 >  2:3101910/1‑150 (MQ=255)
                     aCCTGTTCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGa                                                                                                                                 >  2:830531/1‑150 (MQ=255)
                       cTGTTCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGAcc                                                                                                                               <  2:4855690/150‑1 (MQ=255)
                         gTTCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCTATGCTGGATCGTGAAAAACTGAcccc                                                                                                                             >  2:4448076/1‑150 (MQ=255)
                           tCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGa                                                                                                                           >  2:4668468/1‑150 (MQ=255)
                             aaGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACt                                                                                                                         <  2:1068438/150‑1 (MQ=255)
                             aaGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACt                                                                                                                         <  2:2130749/150‑1 (MQ=255)
                                 tAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCTATGCTGGATCGTGAAAAACTGACCCCGACTTATc                                                                                                                     <  1:440429/150‑1 (MQ=255)
                                  cggCGCCGAGACCATCGTCCCCCGCCGCACAGCGGCCACCCCAGAAGACGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCt                                                                                                                    <  1:4586056/149‑1 (MQ=255)
                                  aGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCt                                                                                                                    <  2:3230703/150‑1 (MQ=255)
                                    gcgcGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTgg                                                                                                                  <  1:1636270/150‑1 (MQ=255)
                                          agaACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGt                                                                                                            >  1:4824336/1‑150 (MQ=255)
                                           gaACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTc                                                                                                           <  1:4617037/150‑1 (MQ=255)
                                           gaACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTc                                                                                                           <  1:1846307/150‑1 (MQ=255)
                                           gaACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTc                                                                                                           >  2:3196608/1‑150 (MQ=255)
                                            aaCATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCTATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCt                                                                                                          <  1:3434321/150‑1 (MQ=255)
                                            aaCATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCt                                                                                                          >  2:2215857/1‑150 (MQ=255)
                                             aCATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTa                                                                                                         >  1:2636726/1‑150 (MQ=255)
                                             aCATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTa                                                                                                         >  1:4163887/1‑150 (MQ=255)
                                              cATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTAt                                                                                                        <  2:3356798/150‑1 (MQ=255)
                                               aTCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATc                                                                                                       >  2:4117890/1‑150 (MQ=255)
                                               aTCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATc                                                                                                       >  2:1299657/1‑150 (MQ=255)
                                                 cTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATcgc                                                                                                     <  1:1485488/150‑1 (MQ=255)
                                                 cTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATcgc                                                                                                     <  2:4488720/150‑1 (MQ=255)
                                                 cTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATcgc                                                                                                     >  2:3474473/1‑150 (MQ=255)
                                                    ccTCGGTCTCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGt                                                                                                  >  2:3610380/1‑150 (MQ=255)
                                                    ccTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGt                                                                                                  >  2:936201/1‑150 (MQ=255)
                                                            gCAAAGCGGCAACCATAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGTATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACacg                                                                                          <  2:2282642/150‑1 (MQ=255)
                                                              aaaGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGTCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACacggt                                                                                        >  1:4253649/1‑150 (MQ=255)
                                                              aaaGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCTATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACacggt                                                                                        >  2:4785362/1‑150 (MQ=255)
                                                              aaaGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACacggt                                                                                        >  2:3761335/1‑150 (MQ=255)
                                                               aaGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACacggta                                                                                       <  1:629318/150‑1 (MQ=255)
                                                                aGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACacggtac                                                                                      >  2:2520868/1‑150 (MQ=255)
                                                                          acagaagaagCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCg                                                                            <  2:4030462/148‑1 (MQ=255)
                                                                          aaaGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCg                                                                            <  1:5261978/150‑1 (MQ=255)
                                                                               agaagCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGa                                                                       <  2:1504340/150‑1 (MQ=255)
                                                                                gaagCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTATGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGAt                                                                      <  2:20139/150‑1 (MQ=255)
                                                                                 aagCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATc                                                                     <  2:856496/150‑1 (MQ=255)
                                                                                    cGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATcgcg                                                                  <  1:3097980/150‑1 (MQ=255)
                                                                                     gATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGt                                                                 >  1:4920862/1‑150 (MQ=255)
                                                                                     gATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGt                                                                 >  1:5191730/1‑150 (MQ=255)
                                                                                      aTTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTa                                                                >  1:48481/1‑150 (MQ=255)
                                                                                      aTTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTa                                                                >  2:4640356/1‑150 (MQ=255)
                                                                                      aTTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTa                                                                >  1:353142/1‑150 (MQ=255)
                                                                                        tCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCTATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACt                                                              >  2:3974825/1‑150 (MQ=255)
                                                                                         cGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTg                                                             >  1:1753892/1‑150 (MQ=255)
                                                                                         cGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTg                                                             >  2:2414273/1‑150 (MQ=255)
                                                                                         cGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTg                                                             >  2:707450/1‑150 (MQ=255)
                                                                                         cGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTg                                                             <  2:4891707/150‑1 (MQ=255)
                                                                                          gTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGa                                                            >  2:3667811/1‑150 (MQ=255)
                                                                                                 tGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAATGATCGCGTACTGAAAACggg                                                     >  1:4032877/1‑150 (MQ=255)
                                                                                                    cGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGcgt                                                  >  2:3064769/1‑150 (MQ=255)
                                                                                                      agcagcTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGcgtcg                                                <  1:995367/150‑1 (MQ=255)
                                                                                                       gcagcTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGcgtcgt                                               >  1:3344587/1‑150 (MQ=255)
                                                                                                        cagcTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCgtg                                              >  2:4222882/1‑150 (MQ=255)
                                                                                                          gcTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGtt                                            <  1:3904090/150‑1 (MQ=255)
                                                                                                          gcTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGtt                                            <  1:2954140/150‑1 (MQ=255)
                                                                                                          gcTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGtt                                            >  1:1698012/1‑150 (MQ=255)
                                                                                                           cTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTc                                           <  1:4089327/150‑1 (MQ=255)
                                                                                                           cTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTc                                           <  1:360134/150‑1 (MQ=255)
                                                                                                               tGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGcc                                       <  1:2748779/150‑1 (MQ=255)
                                                                                                                gAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCa                                      <  1:3846920/150‑1 (MQ=255)
                                                                                                                  aaGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGt                                    >  1:2211754/1‑150 (MQ=255)
                                                                                                                         ttACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGa                             >  1:1001416/1‑150 (MQ=255)
                                                                                                                          tACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCAGACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGaa                            >  1:2196919/1‑150 (MQ=255)
                                                                                                                            cGTTGAGCCGGAATTCGTTCAGGCGATGCTGGATCGTGAAAATCTGACCCCGACTTATCTGGGTGAGTATATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAgg                          <  2:2662845/150‑1 (MQ=255)
                                                                                                                            cGTTGAGCCGGAATACGTTCAGGCTATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAgg                          <  1:2474757/150‑1 (MQ=255)
                                                                                                                               tGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGt                       >  2:3838885/1‑150 (MQ=255)
                                                                                                                               tGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGt                       >  2:167849/1‑150 (MQ=255)
                                                                                                                                gAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTg                      <  1:2521024/150‑1 (MQ=255)
                                                                                                                                gAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTg                      <  1:3891284/150‑1 (MQ=255)
                                                                                                                                 aGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTgc                     >  1:4910571/1‑150 (MQ=255)
                                                                                                                                 aGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTgc                     >  2:1492363/1‑150 (MQ=255)
                                                                                                                                  gCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTgcg                    >  1:1096357/1‑150 (MQ=255)
                                                                                                                                   ccGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTgcgc                   >  2:428794/1‑150 (MQ=255)
                                                                                                                                     ggAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCtt                 >  1:1116522/1‑150 (MQ=255)
                                                                                                                                      gAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTc                >  2:4659613/1‑150 (MQ=255)
                                                                                                                                      gAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTc                >  2:4179889/1‑150 (MQ=255)
                                                                                                                                      gAATACGTTCAGGCGATGCTCGATCGTGAAAAACTGACCCCGACTTATCAGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTc                >  2:2932713/1‑150 (MQ=255)
                                                                                                                                       aaTACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCg               >  1:4055331/1‑150 (MQ=255)
                                                                                                                                           cGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTga           <  2:331722/150‑1 (MQ=255)
                                                                                                                                                aGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTgaagaag      <  1:3130795/150‑1 (MQ=255)
                                                                                                                                                 ggCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTgaagaaga     >  2:4408678/1‑150 (MQ=255)
                                                                                                                                                   cGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGaagaagaag   >  2:3901976/1‑150 (MQ=255)
                                                                                                                                                    gATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAagaagaaga  >  1:1396249/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTCCTCGGTCGCAAAGCGGCAACCAAAGAAGAAGCGATTCGTTTTGCTGGCGAGCAGCTGGTGAAAGGCGGTTACGTTGAGCCGGAATACGTTCAGGCGATGCTGGATCGTGAAAAACTGACCCCGACTTATCTGGGTGAGTCTATCGCGGTGCCACACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGA  >  CP009273/3767094‑3767391

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: