Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 829,009 G→T 94.2% D162Y (GAT→TAT)  dinG → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273829,0090GT94.2% 538.4 / 20.9 194D162Y (GAT→TAT) dinGATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (2/10);  new base T (91/91);  total (93/101)
Fisher's exact test for biased strand distribution p-value = 3.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 194 total aligned reads displayed.

TGGGCGTGGGCGCTACGTTTGTCCGCGTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCGATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCT  >  CP009273/828864‑829155
                                                                                                                                                 |                                                                                                                                                  
      tgggcgCTACGTTTGTCCGCGTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGAt                                                                                                                                            >  2:1440329/1‑148 (MQ=255)
       gggcgCTACGTTTGTCCGCGTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATa                                                                                                                                           <  1:3325550/148‑1 (MQ=255)
          cgcTACGTTTGTCCGCGTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCTAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACtt                                                                                                                                        <  2:1100267/148‑1 (MQ=255)
          cgcTACGTTTGTCCGCGTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACtt                                                                                                                                        <  1:510118/148‑1 (MQ=255)
          cgcTACGTTTGTCCGCGTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACtt                                                                                                                                        <  1:396471/148‑1 (MQ=255)
            cTACGTTTGTCCGCGTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTtat                                                                                                                                      >  1:3059130/1‑148 (MQ=255)
             tACGTTTGTCCGCGTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTtata                                                                                                                                     >  2:3003857/1‑148 (MQ=255)
               cGTTTGTCCGCGTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATaaa                                                                                                                                   <  2:1999520/148‑1 (MQ=255)
                 tttGTCCGCGTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATg                                                                                                                                 <  2:2216593/148‑1 (MQ=255)
                 tttGTCCGCGTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATg                                                                                                                                 >  2:1764006/1‑148 (MQ=255)
                 tttGTCCGCGTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATg                                                                                                                                 >  1:2088728/1‑148 (MQ=255)
                 tttGTCCGCGTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACGTATAAATg                                                                                                                                 >  1:1488376/1‑148 (MQ=255)
                  ttGTCCGCGTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATgg                                                                                                                                >  2:488326/1‑148 (MQ=255)
                   tGTCCGCGTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATggg                                                                                                                               <  1:1149479/148‑1 (MQ=255)
                   tGTCCGCGTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATggg                                                                                                                               >  1:2177479/1‑148 (MQ=255)
                    gTCCGCGTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGa                                                                                                                              >  2:2736506/1‑148 (MQ=255)
                      ccGCGTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATg                                                                                                                            <  2:1348032/148‑1 (MQ=255)
                        gcgTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGt                                                                                                                          <  2:701480/148‑1 (MQ=255)
                           tAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGTCGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCGATCTCGATACTTATAAATGGGATGGTCTg                                                                                                                       >  1:2074068/1‑148 (MQ=255)
                           tAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTg                                                                                                                       <  2:1523515/148‑1 (MQ=255)
                             aTCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTgcg                                                                                                                     <  2:2585777/148‑1 (MQ=255)
                               cTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGcgcg                                                                                                                   <  1:2554739/148‑1 (MQ=255)
                               cTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGcgcg                                                                                                                   >  1:2235096/1‑148 (MQ=255)
                                 gACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTGTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGAt                                                                                                                 >  1:40072/1‑148 (MQ=255)
                                   cGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCGATCTCGATACTTATAAATGGGCTGGTCTGCGCGATCa                                                                                                               >  1:2211276/1‑148 (MQ=255)
                                     gcgcTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCGATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATa                                                                                                             <  1:1045498/148‑1 (MQ=255)
                                        cTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGCTGGTCTGCGCGATCATACTg                                                                                                          >  2:1050522/1‑148 (MQ=255)
                                        cTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTg                                                                                                          <  2:2186344/148‑1 (MQ=255)
                                        cTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTg                                                                                                          <  1:2159841/148‑1 (MQ=255)
                                            ccAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGatat                                                                                                      <  2:1597136/148‑1 (MQ=255)
                                              aGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCg                                                                                                    <  2:392516/148‑1 (MQ=255)
                                                 aCTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCta                                                                                                 <  2:1209146/148‑1 (MQ=255)
                                                 aCTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCta                                                                                                 <  2:1720404/148‑1 (MQ=255)
                                                  cTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCCAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGCTACTTATAAATGGGATGGTCTGCGCGCTCATACTGATATAGCgat                                                                                                >  1:728226/1‑148 (MQ=255)
                                                   tGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCtata                                                                                               >  2:1830606/1‑148 (MQ=255)
                                                    gAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAg                                                                                              <  1:2883138/148‑1 (MQ=255)
                                                    gAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAg                                                                                              >  1:2385128/1‑148 (MQ=255)
                                                      aCCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGAt                                                                                            >  2:1708004/1‑148 (MQ=255)
                                                       cccACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATg                                                                                           >  2:2787129/1‑148 (MQ=255)
                                                        ccACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATTGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGa                                                                                          >  1:2415209/1‑148 (MQ=255)
                                                        ccACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAAAGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGa                                                                                          >  1:626197/1‑148 (MQ=255)
                                                           cGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGa                                                                                       <  1:2990105/148‑1 (MQ=255)
                                                              aaCAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGAtct                                                                                    <  1:3276323/148‑1 (MQ=255)
                                                              aaCAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGAtct                                                                                    <  2:469839/148‑1 (MQ=255)
                                                                cAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATctct                                                                                  <  2:3150076/148‑1 (MQ=255)
                                                                cAGGATCTGCTGGCGTGTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATAATTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATctct                                                                                  >  1:497498/1‑148 (MQ=255)
                                                                   gATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTCTAGATGACGATCTCTGGc                                                                               >  2:3464465/1‑148 (MQ=255)
                                                                            gCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTaa                                                                      <  1:3540553/148‑1 (MQ=255)
                                                                             cGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAg                                                                     <  1:1014035/148‑1 (MQ=255)
                                                                             cGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAg                                                                     >  1:1609552/1‑148 (MQ=255)
                                                                             cGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAg                                                                     <  1:1014002/148‑1 (MQ=255)
                                                                               tttCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTa                                                                   <  1:1866485/148‑1 (MQ=255)
                                                                                ttCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGATCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTAc                                                                  >  2:2519032/1‑148 (MQ=255)
                                                                                ttCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGATCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTAc                                                                  >  2:2519108/1‑148 (MQ=255)
                                                                                   ttGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGa                                                               >  2:774084/1‑148 (MQ=255)
                                                                                    tGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGAc                                                              <  2:3316418/148‑1 (MQ=255)
                                                                                          cgaACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGcc                                                        <  2:1428719/148‑1 (MQ=255)
                                                                                            aaCTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAg                                                      >  1:2336591/1‑148 (MQ=255)
                                                                                            aaCTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAg                                                      >  1:2336355/1‑148 (MQ=255)
                                                                                              cTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAATCCAGCt                                                    >  1:2099809/1‑148 (MQ=255)
                                                                                               tGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCACCTg                                                   >  1:2989861/1‑148 (MQ=255)
                                                                                                gACGCCTAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGc                                                  <  2:346981/148‑1 (MQ=255)
                                                                                                gACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGc                                                  <  2:3276974/148‑1 (MQ=255)
                                                                                                 aCGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGcc                                                 >  1:2561834/1‑148 (MQ=255)
                                                                                                 aCGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGcc                                                 >  1:1946602/1‑148 (MQ=255)
                                                                                                 aCGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCGATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGcc                                                 <  1:288070/148‑1 (MQ=255)
                                                                                                  cGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCt                                                <  2:3552523/148‑1 (MQ=255)
                                                                                                           aTCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGCTACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTCAGTACCGACAAAGCCAGCTGccgcaaccgcaa                                       >  2:2275652/1‑148 (MQ=255)
                                                                                                           aTCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCaa                                       <  1:561484/148‑1 (MQ=255)
                                                                                                             cAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAAtt                                     <  2:2673578/148‑1 (MQ=255)
                                                                                                             cAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTGTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAAtt                                     >  1:751976/1‑148 (MQ=255)
                                                                                                               ggAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGt                                   >  2:694732/1‑148 (MQ=255)
                                                                                                               ggAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCGATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGt                                   <  2:636061/148‑1 (MQ=255)
                                                                                                                 aaGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTta                                 <  2:3128821/148‑1 (MQ=255)
                                                                                                                 aaGAGCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTta                                 >  1:1587171/1‑148 (MQ=255)
                                                                                                                   gagCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTtact                               <  2:2270656/148‑1 (MQ=255)
                                                                                                                   gagCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTtact                               <  2:480233/148‑1 (MQ=255)
                                                                                                                    agCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTtacta                              >  1:2601075/1‑148 (MQ=255)
                                                                                                                    agCAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTtacta                              >  1:2606685/1‑148 (MQ=255)
                                                                                                                      cAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTactact                            <  2:2606685/148‑1 (MQ=255)
                                                                                                                      cAAAAACGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTacaact                            >  1:3132478/1‑148 (MQ=255)
                                                                                                                       aaaaaCGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTActacta                           <  2:164168/148‑1 (MQ=255)
                                                                                                                       aaaaaCGTTGTGCGAAGCTGAAGGGCGATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTActacta                           <  1:388438/148‑1 (MQ=255)
                                                                                                                       aaaaaCGTTGTGCGAAGCTGAAGGGCGATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTActacta                           <  1:1835373/148‑1 (MQ=255)
                                                                                                                        aaaaCGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTAt                          <  2:99673/148‑1 (MQ=255)
                                                                                                                            cGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTg                      <  1:1705958/148‑1 (MQ=255)
                                                                                                                            cGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTg                      <  1:1348452/148‑1 (MQ=255)
                                                                                                                            cGTTGTGCGAAGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTg                      <  1:1291442/148‑1 (MQ=255)
                                                                                                                                     aaGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGt             <  1:2957794/148‑1 (MQ=255)
                                                                                                                                     aaGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGt             <  1:1764006/148‑1 (MQ=255)
                                                                                                                                     aaGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTGTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAAATGTTACTACTATCGTGAATGCCCGt             >  1:34930/1‑148 (MQ=255)
                                                                                                                                      aGCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGtt            <  2:1299676/148‑1 (MQ=255)
                                                                                                                                       gCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGttt           <  1:2649877/148‑1 (MQ=255)
                                                                                                                                       gCTGAAGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGttt           <  1:864126/148‑1 (MQ=255)
                                                                                                                                           aaGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTg       >  1:229497/1‑148 (MQ=255)
                                                                                                                                           aaGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTg       >  2:3114548/1‑148 (MQ=255)
                                                                                                                                           aaGGGCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTg       >  2:1914277/1‑148 (MQ=255)
                                                                                                                                             gggCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTc     >  2:1812130/1‑148 (MQ=255)
                                                                                                                                              ggCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCg    >  2:2131497/1‑148 (MQ=255)
                                                                                                                                              ggCTATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCg    >  2:1847697/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TGGGCGTGGGCGCTACGTTTGTCCGCGTAATCTGACGGCGCTCGCCAGTACTGAACCCACGCAACAGGATCTGCTGGCGTTTCTTGACGACGAACTGACGCCGAACAATCAGGAAGAGCAAAAACGTTGTGCGAAGCTGAAGGGCGATCTCGATACTTATAAATGGGATGGTCTGCGCGATCATACTGATATCGCTATAGATGACGATCTCTGGCGTCGTTTAAGTACCGACAAAGCCAGCTGCCTCAACCGCAATTGTTACTACTATCGTGAATGCCCGTTTTTTGTCGCT  >  CP009273/828864‑829155

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: