Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,563,702 T→A 10.1% E326D (GAA→GAT glgX ← glycogen debranching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,563,7020TA10.1% 578.0 / 47.3 197E326D (GAA→GATglgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base T (73/104);  new base A (7/13);  total (80/117)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 201 total aligned reads displayed.

ATATCCCACGGTTCAGCAATTAACTTCACCTGCGAGAGCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGCTACGGTTATCGATCCCACGCAGCGAG  >  CP009273/3563556‑3563849
                                                                                                                                                  |                                                                                                                                                   
   tCCCACGGTTCAGCAATTAACTTCACCTGCGAGAGCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTa                                                                                                                                                 >  1:123341/1‑148 (MQ=255)
   tCCCACGGTTCAGCAATTAACTTCACCTGCGAGAGCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTa                                                                                                                                                 <  1:1986807/148‑1 (MQ=255)
   tCCCACGGTTCAGCAATTAACTTCACCTGCGAGAGCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTa                                                                                                                                                 >  2:3222858/1‑148 (MQ=255)
              aGCAATTAACTTCACCTGCGAGAGCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACgc                                                                                                                                      <  2:3107357/148‑1 (MQ=255)
              aGCAATTAACTTCACCTGCGAGAGCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACgc                                                                                                                                      >  2:1250226/1‑148 (MQ=255)
               gCAATTAACTTCACCTGCGAGAGCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACgca                                                                                                                                     <  2:1172847/148‑1 (MQ=255)
                   ttAACTTCACCTGCGAGAGCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACgcaggca                                                                                                                                 <  2:3279616/148‑1 (MQ=255)
                    tAACTTCACCTGCGAGAGCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACgcaggcag                                                                                                                                <  1:1536756/148‑1 (MQ=255)
                    tAACTTCACCTGCGAGAGCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACgcaggcag                                                                                                                                <  1:3255867/148‑1 (MQ=255)
                      aCTTCACCTGCGAGAGCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCaggcaggc                                                                                                                              >  1:3302224/1‑148 (MQ=255)
                       cTTCACCTGCGAGAGCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGgcg                                                                                                                             <  2:1526495/148‑1 (MQ=255)
                           aCCTGCGAGAGCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTgg                                                                                                                         >  1:380280/1‑148 (MQ=255)
                             cTGCGAGAGCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCa                                                                                                                       <  1:937165/148‑1 (MQ=255)
                               gCGAGAGCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGAGGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATa                                                                                                                     <  2:153604/148‑1 (MQ=255)
                                  agagCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATc                                                                                                                  >  1:257210/1‑148 (MQ=255)
                                  agagCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTATCTACCCAATAACGCAGGCAGGCGCTGGCATAATc                                                                                                                  >  1:57389/1‑148 (MQ=255)
                                   gagCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATcc                                                                                                                 <  1:2811137/148‑1 (MQ=255)
                                   gagCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATcc                                                                                                                 <  1:2426298/148‑1 (MQ=255)
                                    agCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATcca                                                                                                                >  2:861265/1‑148 (MQ=255)
                                            gCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCAccgcc                                                                                                        >  2:390165/1‑148 (MQ=255)
                                                 tCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATg                                                                                                   <  1:598033/148‑1 (MQ=255)
                                                 tCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATg                                                                                                   >  2:583581/1‑148 (MQ=255)
                                                    ggATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACt                                                                                                >  2:3258348/1‑148 (MQ=255)
                                                    ggATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACt                                                                                                <  2:1590688/148‑1 (MQ=255)
                                                     gATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTc                                                                                               >  2:560822/1‑148 (MQ=255)
                                                        aGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCaaa                                                                                            >  1:3102462/1‑148 (MQ=255)
                                                        aGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCaaa                                                                                            >  2:1342427/1‑148 (MQ=255)
                                                         gCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAAt                                                                                           >  1:3079158/1‑148 (MQ=255)
                                                         gCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAAt                                                                                           <  1:3222858/148‑1 (MQ=255)
                                                         gCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAAt                                                                                           >  2:3171385/1‑148 (MQ=255)
                                                               aacaacGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTATCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCg                                                                                     >  2:1987644/1‑148 (MQ=255)
                                                                acaacGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCgt                                                                                    >  2:10062/1‑148 (MQ=255)
                                                                 caacGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTATCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCCCCGCCGGATGACTCAAATTGAGCgtg                                                                                   <  2:60252/148‑1 (MQ=255)
                                                                    cGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTg                                                                                >  2:2317591/1‑148 (MQ=255)
                                                                      gcgcATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTATCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGcc                                                                              >  2:1787973/1‑148 (MQ=255)
                                                                      gcgcATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTATCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGcc                                                                              >  2:2748840/1‑148 (MQ=255)
                                                                            ccTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGACATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAAcc                                                                        >  2:182849/1‑148 (MQ=255)
                                                                            ccTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTATCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAAcc                                                                        <  1:2183387/148‑1 (MQ=255)
                                                                             cTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCg                                                                       <  1:808258/148‑1 (MQ=255)
                                                                              tGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCgg                                                                      >  1:630930/1‑148 (MQ=255)
                                                                                aCGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTc                                                                    <  1:1904249/148‑1 (MQ=255)
                                                                                   gAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAg                                                                 >  1:847005/1‑148 (MQ=255)
                                                                                    aaCTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGt                                                                >  2:2860608/1‑148 (MQ=255)
                                                                                    aaCTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGt                                                                <  2:1071003/148‑1 (MQ=255)
                                                                                       tctGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTtgt                                                             <  2:3282192/148‑1 (MQ=255)
                                                                                        ctGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTtgtg                                                            >  2:568203/1‑148 (MQ=255)
                                                                                        ctGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTtgtg                                                            <  1:2487654/148‑1 (MQ=255)
                                                                                           gCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATa                                                         <  2:3170465/148‑1 (MQ=255)
                                                                                           gCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATa                                                         <  1:1726403/148‑1 (MQ=255)
                                                                                           gCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATa                                                         <  2:3170258/148‑1 (MQ=255)
                                                                                            cGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATaa                                                        <  2:667166/148‑1 (MQ=255)
                                                                                            cGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATaa                                                        <  2:2672234/148‑1 (MQ=255)
                                                                                             gTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAAt                                                       <  2:270285/148‑1 (MQ=255)
                                                                                                   cccATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGt                                                 <  1:703124/148‑1 (MQ=255)
                                                                                                    ccATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGtc                                                <  1:937081/148‑1 (MQ=255)
                                                                                                    ccATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGtc                                                <  2:2884888/148‑1 (MQ=255)
                                                                                                    ccATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGtc                                                <  2:1441686/148‑1 (MQ=255)
                                                                                                    ccATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTATCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGtc                                                <  1:565703/148‑1 (MQ=255)
                                                                                                    ccATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTATCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGtc                                                <  1:927689/148‑1 (MQ=255)
                                                                                                     cATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGtct                                               <  2:3146908/148‑1 (MQ=255)
                                                                                                      aTGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGtctt                                              <  2:1036384/148‑1 (MQ=255)
                                                                                                      aTGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGtctt                                              >  1:493420/1‑148 (MQ=255)
                                                                                                       tGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGtcttc                                             <  2:239028/148‑1 (MQ=255)
                                                                                                         aCTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTtctc                                           <  2:1144086/148‑1 (MQ=255)
                                                                                                          cTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTctct                                          <  1:10062/148‑1 (MQ=255)
                                                                                                           tGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCtt                                         >  1:1668668/1‑148 (MQ=255)
                                                                                                           tGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCtt                                         <  2:1344751/148‑1 (MQ=255)
                                                                                                            gccgccAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTa                                        <  1:1365527/148‑1 (MQ=255)
                                                                                                                 cAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCaa                                   <  2:2483780/148‑1 (MQ=255)
                                                                                                                  aGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCaat                                  <  1:3260903/148‑1 (MQ=255)
                                                                                                                   gATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCaata                                 <  1:2330390/148‑1 (MQ=255)
                                                                                                                    aTCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCaataa                                <  2:2036998/148‑1 (MQ=255)
                                                                                                                    aTCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCaataa                                <  1:2597145/148‑1 (MQ=255)
                                                                                                                    aTCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCaataa                                >  2:1399633/1‑148 (MQ=255)
                                                                                                                     tCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCaataat                               >  2:1628991/1‑148 (MQ=255)
                                                                                                                     tCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCaataat                               >  2:2980689/1‑148 (MQ=255)
                                                                                                                      cAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAataata                              <  1:2461204/148‑1 (MQ=255)
                                                                                                                      cAAAGCGGAAACCATCGACGTGGCAGGTATCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAataata                              <  2:1914840/148‑1 (MQ=255)
                                                                                                                       aaaGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAg                             <  2:1952594/148‑1 (MQ=255)
                                                                                                                       aaaGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAg                             <  1:2987403/148‑1 (MQ=255)
                                                                                                                        aaGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTtctc                                           <  2:187270/133‑1 (MQ=255)
                                                                                                                        aaGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTtctc                                           >  1:1858012/1‑133 (MQ=255)
                                                                                                                        aaGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGc                            <  1:2014881/148‑1 (MQ=255)
                                                                                                                          gCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGCTa                          <  2:3140331/148‑1 (MQ=255)
                                                                                                                             gAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGCTACgg                       >  2:283216/1‑148 (MQ=255)
                                                                                                                             gAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGCTACgg                       <  1:1790858/148‑1 (MQ=255)
                                                                                                                             gAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGACATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGCTACgg                       <  1:1250578/148‑1 (MQ=255)
                                                                                                                              aaaCCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGCTACGGt                      <  1:486246/148‑1 (MQ=255)
                                                                                                                              aaaCCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGCTACGGt                      <  1:2238122/148‑1 (MQ=255)
                                                                                                                                aCCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGCTACGGTTa                    >  1:1488889/1‑148 (MQ=255)
                                                                                                                                aCCATCGACGTGGCAGGTATCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGCTACGGTTa                    >  1:2163487/1‑148 (MQ=255)
                                                                                                                                    tCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGCTACGGTTATCGa                >  2:2103800/1‑148 (MQ=255)
                                                                                                                                       aCGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGCTACGGTTATCGATcc             <  1:2081395/148‑1 (MQ=255)
                                                                                                                                          tGGCAGGTATCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGCTACGGTTATCGATCCCAc          >  1:3132408/1‑148 (MQ=255)
                                                                                                                                             cAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGCTACGGTTATCGATCCCACGCa       >  1:2535490/1‑148 (MQ=255)
                                                                                                                                               ggTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGCTACGGTTATCGATCCCACGCAGc     >  2:2380781/1‑148 (MQ=255)
                                                                                                                                                gattCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGCTAAGGTTATCGATCCCACGCAGCg    >  1:520165/3‑148 (MQ=255)
                                                                                                                                                gTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGCTACGGTTATCGATCCCACGCAGCg    >  2:3043751/1‑148 (MQ=255)
                                                                                                                                                gTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGCTACGGTTATCGATCCCACGCAGCg    >  1:2185971/1‑148 (MQ=255)
                                                                                                                                                  ttCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGCTACGGTTATCGATCCCACGCAGCgag  >  1:2647402/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATATCCCACGGTTCAGCAATTAACTTCACCTGCGAGAGCACCGGGCAGTTCTGGATAGCGGTAAACAACGGCGCATCCTGACGGAACTCTGGCGTACGGCCCATGACTGCCGCCAGATCAAAGCGGAAACCATCGACGTGGCAGGTTTCTACCCAATAACGCAGGCAGGCGCTGGCATAATCCACCACCGCCGGATGACTCAAATTGAGCGTGTTGCCGCAACCGGTCCAGTTGTGATAATCGCCGTCTTCTCTTATCCAATAATAGCTACGGTTATCGATCCCACGCAGCGAG  >  CP009273/3563556‑3563849

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: