Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 784,094 G→T 86.6% A67E (GCG→GAG)  galM ← aldose 1‑epimerase; type‑1 mutarotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273784,0940GT86.6% 452.4 / 73.7 194A67E (GCG→GAG) galMaldose 1‑epimerase; type‑1 mutarotase
Reads supporting (aligned to +/- strand):  ref base G (6/20);  new base T (94/74);  total (100/94)
Fisher's exact test for biased strand distribution p-value = 2.62e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 195 total aligned reads displayed.

AGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGCGTTGTTACGCAAAG  >  CP009273/783951‑784238
                                                                                                                                               |                                                                                                                                                
 gCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGcc                                                                                                                                             <  1:735736/148‑1 (MQ=255)
 gCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGcc                                                                                                                                             <  2:1095320/148‑1 (MQ=255)
    aCGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGcctga                                                                                                                                          >  1:1748863/1‑148 (MQ=255)
    aCGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGcctga                                                                                                                                          >  1:1196692/1‑148 (MQ=255)
    aCGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGcctga                                                                                                                                          <  2:1279548/148‑1 (MQ=255)
    aCGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGcctga                                                                                                                                          <  2:2004025/148‑1 (MQ=255)
     cGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGcctgat                                                                                                                                         >  2:2972401/1‑148 (MQ=255)
      gTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGcctgatc                                                                                                                                        <  1:2222173/148‑1 (MQ=255)
        ttGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGcctgatcct                                                                                                                                      >  2:60502/1‑148 (MQ=255)
        ttGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGcctgatcct                                                                                                                                      <  2:2764976/148‑1 (MQ=255)
         tGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGcctgatcctg                                                                                                                                     <  1:343397/148‑1 (MQ=255)
         tGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGcctgatcctg                                                                                                                                     <  2:2171086/148‑1 (MQ=255)
                cccTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCAt                                                                                                                              >  2:1037196/1‑148 (MQ=255)
                 ccTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCAtt                                                                                                                             >  1:2800329/1‑148 (MQ=255)
                      cGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCggg                                                                                                                        <  1:1909484/148‑1 (MQ=255)
                      cGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCggg                                                                                                                        <  1:1037196/148‑1 (MQ=255)
                       ggCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGc                                                                                                                       <  1:24090/148‑1 (MQ=255)
                        gCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCt                                                                                                                      <  2:2924953/148‑1 (MQ=255)
                          ccgccgTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTgg                                                                                                                    <  2:236260/148‑1 (MQ=255)
                          ccgccgTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTgg                                                                                                                    <  1:3246549/148‑1 (MQ=255)
                           cgccgTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGc                                                                                                                   >  1:2248666/1‑148 (MQ=255)
                             ccgTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGcac                                                                                                                 <  2:1533597/148‑1 (MQ=255)
                                tGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGc                                                                                                              <  2:1363518/148‑1 (MQ=255)
                                tGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGc                                                                                                              <  2:1472867/148‑1 (MQ=255)
                                tGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGc                                                                                                              >  2:972442/1‑148 (MQ=255)
                                tGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGc                                                                                                              >  2:973036/1‑148 (MQ=255)
                                     cTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGagc                                                                                                         >  1:2727532/1‑148 (MQ=255)
                                       ggTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGagcag                                                                                                       <  2:1339118/148‑1 (MQ=255)
                                           aaCGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGcc                                                                                                   <  1:1532682/148‑1 (MQ=255)
                                            aCGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCt                                                                                                  <  2:1201865/148‑1 (MQ=255)
                                            aCGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCt                                                                                                  >  1:1796425/1‑148 (MQ=255)
                                            aCGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCt                                                                                                  >  1:442678/1‑148 (MQ=255)
                                                 cTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCgg                                                                                             <  1:147149/148‑1 (MQ=255)
                                                  tGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGa                                                                                            >  1:2332421/1‑148 (MQ=255)
                                                  tGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGa                                                                                            <  1:2103035/148‑1 (MQ=255)
                                                  tGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGa                                                                                            <  1:2103921/148‑1 (MQ=255)
                                                   gACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGAc                                                                                           >  2:1923145/1‑148 (MQ=255)
                                                   gACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGAc                                                                                           >  2:2860814/1‑148 (MQ=255)
                                                    aCTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACg                                                                                          >  2:2935353/1‑148 (MQ=255)
                                                     cTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGc                                                                                         <  2:3207805/148‑1 (MQ=255)
                                                         gCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGcc                                                                                     >  1:1106300/1‑148 (MQ=255)
                                                           gAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCAt                                                                                   <  2:324137/148‑1 (MQ=255)
                                                              aGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCgg                                                                                <  2:178775/148‑1 (MQ=255)
                                                                cGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGaa                                                                              >  2:354854/1‑148 (MQ=255)
                                                                  tCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCCTCGGAAAg                                                                            >  1:2962629/1‑148 (MQ=255)
                                                                  tCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAg                                                                            >  2:705033/1‑148 (MQ=255)
                                                                  tCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAg                                                                            <  1:635488/148‑1 (MQ=255)
                                                                  tCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAg                                                                            <  2:1636345/148‑1 (MQ=255)
                                                                   cACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGc                                                                           >  1:2269297/1‑148 (MQ=255)
                                                                     cGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCgg                                                                         >  1:250633/1‑148 (MQ=255)
                                                                     cGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCgg                                                                         >  1:3280906/1‑148 (MQ=255)
                                                                      ggTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGa                                                                        <  1:1754074/148‑1 (MQ=255)
                                                                       gTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGAGAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGaa                                                                       <  1:12282/148‑1 (MQ=255)
                                                                        tttCACCGTCAAAGGTATACCGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAAt                                                                      <  2:2130362/148‑1 (MQ=255)
                                                                        tttCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAAt                                                                      <  2:2334085/148‑1 (MQ=255)
                                                                        tttCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAAt                                                                      <  2:877265/148‑1 (MQ=255)
                                                                         ttCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGTCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCTAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATa                                                                     >  1:2741247/1‑148 (MQ=255)
                                                                            aCCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACgg                                                                  >  2:2582216/1‑148 (MQ=255)
                                                                             ccGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATATCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACggg                                                                 <  1:1327047/148‑1 (MQ=255)
                                                                                tCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATAGGAAACCCGAAAAAGGGCgg                                                              >  1:1489545/1‑148 (MQ=255)
                                                                                 cAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGa                                                             >  2:1425464/1‑148 (MQ=255)
                                                                                   aaGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaa                                                           >  2:249957/1‑148 (MQ=255)
                                                                                   aaGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaa                                                           <  1:270461/148‑1 (MQ=255)
                                                                                    aGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaag                                                          >  2:854048/1‑148 (MQ=255)
                                                                                     ggTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagt                                                         >  1:1413385/1‑148 (MQ=255)
                                                                                     ggTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagt                                                         <  2:526820/148‑1 (MQ=255)
                                                                                     ggTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagt                                                         <  1:1011448/148‑1 (MQ=255)
                                                                                      gTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagta                                                        >  1:2035731/1‑148 (MQ=255)
                                                                                      gTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagta                                                        >  1:902227/1‑148 (MQ=255)
                                                                                       tataACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagtaa                                                       >  1:403037/1‑148 (MQ=255)
                                                                                       tataACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagtaa                                                       <  2:1413385/148‑1 (MQ=255)
                                                                                        ataACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagtaaa                                                      <  1:2563227/148‑1 (MQ=255)
                                                                                        ataACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagtaaa                                                      <  1:1689646/148‑1 (MQ=255)
                                                                                        ataACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagtaaa                                                      <  1:1689662/148‑1 (MQ=255)
                                                                                        ataACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCAGaaagtaaa                                                      <  1:60402/148‑1 (MQ=255)
                                                                                         taACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagtaaag                                                     >  2:2695295/1‑148 (MQ=255)
                                                                                         taACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagtaaag                                                     >  2:448664/1‑148 (MQ=255)
                                                                                         taACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagtaaag                                                     >  2:3030342/1‑148 (MQ=255)
                                                                                          aaCGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagtaaagt                                                    >  2:714353/1‑148 (MQ=255)
                                                                                                       gATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTcc                                       <  2:3121398/148‑1 (MQ=255)
                                                                                                        aTACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCa                                      <  2:3098735/148‑1 (MQ=255)
                                                                                                        aTACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCa                                      <  1:3111774/148‑1 (MQ=255)
                                                                                                         tACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCAt                                     <  2:295180/148‑1 (MQ=255)
                                                                                                         tACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCAt                                     <  2:1390919/148‑1 (MQ=255)
                                                                                                              ttGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCg                                <  2:240696/148‑1 (MQ=255)
                                                                                                              ttGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCg                                <  2:716285/148‑1 (MQ=255)
                                                                                                              ttGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCg                                <  2:2376444/148‑1 (MQ=255)
                                                                                                              ttGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCg                                <  1:1280633/148‑1 (MQ=255)
                                                                                                                ggCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTg                              <  2:1805678/148‑1 (MQ=255)
                                                                                                                ggCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTg                              <  2:873905/148‑1 (MQ=255)
                                                                                                                 gCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGa                             <  2:1483905/148‑1 (MQ=255)
                                                                                                                   aTAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACt                           <  2:1124996/148‑1 (MQ=255)
                                                                                                                   aTAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACt                           <  1:3213971/148‑1 (MQ=255)
                                                                                                                               agagGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGCg               <  1:2808740/148‑1 (MQ=255)
                                                                                                                                gagGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGCgt              >  1:1351434/1‑148 (MQ=255)
                                                                                                                                   gCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGCgttgt           >  2:1668714/1‑148 (MQ=255)
                                                                                                                                      cccAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCCTCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGCGTTGTTAc        >  2:469110/1‑148 (MQ=255)
                                                                                                                                        cAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCAGCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGCGTTGTTACGc      >  1:1720647/1‑148 (MQ=255)
                                                                                                                                          gAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGCGTTGTTACGCaa    <  2:504801/148‑1 (MQ=255)
                                                                                                                                            aaCGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGCGTTGTTACGCAAAg  <  1:2972401/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
AGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGCGTTGTTACGCAAAG  >  CP009273/783951‑784238

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: ‑