Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,394,875 G→T 94.4% A124A (GCG→GCT purA → adenylosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,394,8750GT94.4% 752.7 / 32.0 258A124A (GCG→GCTpurAadenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base G (5/10);  new base T (108/135);  total (113/145)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 258 total aligned reads displayed.

CTCCGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTCC  >  CP009273/4394730‑4395017
                                                                                                                                                 |                                                                                                                                              
ctcCGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTaa                                                                                                                                              <  1:3138525/148‑1 (MQ=255)
     gCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGa                                                                                                                                         <  1:2619828/148‑1 (MQ=255)
     gCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGa                                                                                                                                         <  2:373342/148‑1 (MQ=255)
      ccGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGAt                                                                                                                                        >  2:3244358/1‑148 (MQ=255)
              atgaaagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGcacca                                                                                                                                >  1:3311558/1‑148 (MQ=255)
              atgaaagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGcacca                                                                                                                                <  2:384120/148‑1 (MQ=255)
               tgaaagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGcaccac                                                                                                                               >  2:687481/1‑148 (MQ=255)
                gaaagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCaccacc                                                                                                                              >  1:2702585/1‑148 (MQ=255)
                 aaagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCg                                                                                                                             >  1:1946145/1‑148 (MQ=255)
                   agagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGt                                                                                                                           >  2:1310515/1‑148 (MQ=255)
                      gatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGt                                                                                                                        >  1:2977633/1‑148 (MQ=255)
                       atgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTg                                                                                                                       >  1:2506109/1‑148 (MQ=255)
                       atgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTg                                                                                                                       <  1:589739/148‑1 (MQ=255)
                        tgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTgg                                                                                                                      >  1:3172151/1‑148 (MQ=255)
                        tgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTgg                                                                                                                      <  1:1675353/148‑1 (MQ=255)
                        tgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTgg                                                                                                                      <  1:1980389/148‑1 (MQ=255)
                        tgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTgg                                                                                                                      >  1:2338110/1‑148 (MQ=255)
                           aagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTAt                                                                                                                   <  2:1383844/148‑1 (MQ=255)
                            agaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATc                                                                                                                  >  2:1786275/1‑148 (MQ=255)
                               aCTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCggg                                                                                                               <  2:2393148/148‑1 (MQ=255)
                                 tGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGcc                                                                                                             >  2:494344/1‑148 (MQ=255)
                                  ggAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCt                                                                                                            >  2:1420806/1‑148 (MQ=255)
                                   gAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTg                                                                                                           >  2:3124383/1‑148 (MQ=255)
                                   gAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTg                                                                                                           >  1:2053167/1‑148 (MQ=255)
                                   gAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTg                                                                                                           >  1:2029940/1‑148 (MQ=255)
                                      gACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCtt                                                                                                        <  1:485342/148‑1 (MQ=255)
                                      gACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCtt                                                                                                        >  1:2219213/1‑148 (MQ=255)
                                       aCCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTa                                                                                                       <  2:1668252/148‑1 (MQ=255)
                                            ggCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAg                                                                                                  <  1:2648553/148‑1 (MQ=255)
                                              cATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGAt                                                                                                >  2:1609261/1‑148 (MQ=255)
                                              cATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGAt                                                                                                >  1:3402615/1‑148 (MQ=255)
                                                 ccccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATaaa                                                                                             <  1:2390349/148‑1 (MQ=255)
                                                 ccccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATaaa                                                                                             <  1:1879112/148‑1 (MQ=255)
                                                 ccccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATaaa                                                                                             >  2:2135557/1‑148 (MQ=255)
                                                  cccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAg                                                                                            <  2:2774178/148‑1 (MQ=255)
                                                  cccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAg                                                                                            <  2:3112313/148‑1 (MQ=255)
                                                    cGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTa                                                                                          <  1:2903828/148‑1 (MQ=255)
                                                       tCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCa                                                                                       <  2:603151/148‑1 (MQ=255)
                                                       tCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCa                                                                                       <  2:2913916/148‑1 (MQ=255)
                                                       tCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGAGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCa                                                                                       >  1:3552637/1‑148 (MQ=255)
                                                         gTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACg                                                                                     <  2:2290529/148‑1 (MQ=255)
                                                         gTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACg                                                                                     <  2:61020/148‑1 (MQ=255)
                                                         gTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACg                                                                                     <  2:1260413/148‑1 (MQ=255)
                                                          tGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGt                                                                                    >  2:2333893/1‑148 (MQ=255)
                                                            aGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTcg                                                                                  >  1:1295796/1‑148 (MQ=255)
                                                              cGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTcgcg                                                                                >  1:1661967/1‑148 (MQ=255)
                                                                tCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGt                                                                              >  1:98405/1‑148 (MQ=255)
                                                                  tgctgctgTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCt                                                                            <  1:931405/148‑1 (MQ=255)
                                                                   gctgctgTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTg                                                                           >  2:3351758/1‑148 (MQ=255)
                                                                      gctgTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCgt                                                                        >  1:843287/1‑148 (MQ=255)
                                                                        tgTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCgtgt                                                                      <  2:2977633/148‑1 (MQ=255)
                                                                        tgTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCgtgt                                                                      <  2:1981899/148‑1 (MQ=255)
                                                                           cTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTgg                                                                   >  2:3418105/1‑148 (MQ=255)
                                                                               aGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGAc                                                               <  2:1539202/148‑1 (MQ=255)
                                                                               aGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGAc                                                               <  2:3403879/148‑1 (MQ=255)
                                                                                gCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGAcc                                                              <  1:259050/148‑1 (MQ=255)
                                                                                 cATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCt                                                             >  2:1044093/1‑148 (MQ=255)
                                                                                    gTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCtttt                                                          >  2:3108175/1‑148 (MQ=255)
                                                                                    gTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCtttt                                                          <  1:2799783/148‑1 (MQ=255)
                                                                                    gTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCtttt                                                          <  1:1530022/148‑1 (MQ=255)
                                                                                          tGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACaa                                                    <  1:483917/148‑1 (MQ=255)
                                                                                              ccTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGaa                                                >  1:2858967/1‑148 (MQ=255)
                                                                                              ccTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGaa                                                >  2:1969243/1‑148 (MQ=255)
                                                                                                ttGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAAc                                              <  1:1123204/148‑1 (MQ=255)
                                                                                                   aTTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCtt                                           <  2:760330/148‑1 (MQ=255)
                                                                                                   aTTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCtt                                           <  2:1861208/148‑1 (MQ=255)
                                                                                                    ttATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTc                                          >  2:2795838/1‑148 (MQ=255)
                                                                                                     tATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCg                                         <  1:95365/148‑1 (MQ=255)
                                                                                                      aTCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGc                                        <  1:1634474/148‑1 (MQ=255)
                                                                                                      aTCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGc                                        >  1:1733332/1‑148 (MQ=255)
                                                                                                       tCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCt                                       >  2:1867102/1‑148 (MQ=255)
                                                                                                          cGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGaa                                    >  1:3281840/1‑148 (MQ=255)
                                                                                                          cGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGaa                                    >  1:3534439/1‑148 (MQ=255)
                                                                                                            ttGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGaaaa                                  <  1:504763/148‑1 (MQ=255)
                                                                                                             tGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGaaaaa                                 <  1:1786275/148‑1 (MQ=255)
                                                                                                             tGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGaaaaa                                 <  1:1953678/148‑1 (MQ=255)
                                                                                                              gcgcTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAAc                                <  1:3244358/148‑1 (MQ=255)
                                                                                                                 cTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGa                             >  2:2883054/1‑148 (MQ=255)
                                                                                                                  tGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGaa                            <  2:2506109/148‑1 (MQ=255)
                                                                                                                    gATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAaag                          <  2:3043151/148‑1 (MQ=255)
                                                                                                                    gATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTAACGGGCCAGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTACGCTGAAAAACTGAaag                          >  1:2159406/1‑148 (MQ=255)
                                                                                                                       aaCGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAaagaag                       >  2:1459676/1‑148 (MQ=255)
                                                                                                                       aaCGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAaagaag                       <  2:3544375/148‑1 (MQ=255)
                                                                                                                        aCGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGt                      <  1:2734241/148‑1 (MQ=255)
                                                                                                                          gcgcgTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGa                    <  2:1298696/148‑1 (MQ=255)
                                                                                                                           cgcgTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGAt                   <  2:2457927/148‑1 (MQ=255)
                                                                                                                              gTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGa                <  2:620595/148‑1 (MQ=255)
                                                                                                                               tGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGaa               <  1:1377109/148‑1 (MQ=255)
                                                                                                                                gagaAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAat              <  1:135130/148‑1 (MQ=255)
                                                                                                                                gagaAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAat              >  2:1147324/1‑148 (MQ=255)
                                                                                                                                 agaAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAata             <  2:411270/148‑1 (MQ=255)
                                                                                                                                   aaaGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATc           >  2:1834589/1‑148 (MQ=255)
                                                                                                                                   aaaGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATc           >  2:1834661/1‑148 (MQ=255)
                                                                                                                                    ccgcgcGTGGCGCGCCAGCGATCGCCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGGTGCCGACATTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATca          <  2:1941229/146‑1 (MQ=255)
                                                                                                                                     aGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATcac         >  2:3243813/1‑148 (MQ=255)
                                                                                                                                      gcgcGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATcaca        <  2:361597/148‑1 (MQ=255)
                                                                                                                                       cgcgTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACaa       <  2:432593/148‑1 (MQ=255)
                                                                                                                                       cgcgTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACaa       >  2:2579586/1‑148 (MQ=255)
                                                                                                                                          gTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACtt    <  2:2195640/148‑1 (MQ=255)
                                                                                                                                            ggCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTcc  <  1:687481/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CTCCGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTCC  >  CP009273/4394730‑4395017

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: