Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,029,664 C→T 6.2% L319L (CTG→TTG)  hyaB → hydrogenase 1, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,029,6640CT 6.2% 1230.5 / 44.4 340L319L (CTG→TTG) hyaBhydrogenase 1, large subunit
Reads supporting (aligned to +/- strand):  ref base C (165/154);  new base T (8/13);  total (173/167)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 355 total aligned reads displayed.

TTCATTAACAACGTGATGATCCCCGACGCCTTAGCCATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGG  >  CP009273/1029517‑1029812
                                                                                                                                                   |                                                                                                                                                    
  cATTAACAACGTGATGATCCCCGACGCCTTAGCCATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTctgc                                                                                                                                                   >  1:1293896/1‑149 (MQ=255)
  cATTAACAACGTGATGATCCCCGACGCCTTAGCCATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCACATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTctgc                                                                                                                                                   <  1:11345/149‑1 (MQ=255)
     tAACAACGTGATGATCCCCGACGCCTTAGCCATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGa                                                                                                                                                <  2:1552669/149‑1 (MQ=255)
      aacaacGTGATGATCCCCGACGCCTTAGCCATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTTTGCTGAt                                                                                                                                               >  2:3636979/1‑149 (MQ=255)
      aacaacGTGATGATCCCCGACGCCTTAGCCATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGAt                                                                                                                                               >  2:1107970/1‑149 (MQ=255)
        caacGTGATGATCCCCGACGCCTTAGCCATCGGTCAGTTCAACAACCCGTGGAGCTAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTTTGCTGATGc                                                                                                                                             <  1:3928783/149‑1 (MQ=255)
          acGTGATGATCCCCGACGCCTTAGCCATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCt                                                                                                                                           <  1:5194977/149‑1 (MQ=255)
           cGTGATGATCCCCGACGCCTTAGCCATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTg                                                                                                                                          <  2:892629/149‑1 (MQ=255)
            gTGATGATCCCCGACGCCTTAGCCATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTgg                                                                                                                                         <  1:2476700/149‑1 (MQ=255)
                    ccccGACGCCTTAGCCATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCgg                                                                                                                                 >  2:2437698/1‑149 (MQ=255)
                     cccGACGCCTTAGCCATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGt                                                                                                                                <  2:5092213/149‑1 (MQ=255)
                      ccGACGCCTTAGCCATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTg                                                                                                                               >  1:3139395/1‑149 (MQ=255)
                       cGACGCCTTAGCCATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGa                                                                                                                              <  1:3399895/149‑1 (MQ=255)
                           gCCTTAGCCATCGGTCAGTTCAACAACCCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTaa                                                                                                                          <  2:780782/149‑1 (MQ=255)
                           gCCTTAGCCATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTaa                                                                                                                          >  1:994769/1‑149 (MQ=255)
                           gCCTTAGCCATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTaa                                                                                                                          >  1:1444000/1‑149 (MQ=255)
                              ttAGCCATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACgg                                                                                                                       <  1:4916833/149‑1 (MQ=255)
                                 gCCATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGa                                                                                                                    >  1:3582061/1‑149 (MQ=255)
                                   cATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACt                                                                                                                  <  2:4613908/149‑1 (MQ=255)
                                   cATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACt                                                                                                                  >  2:263818/1‑149 (MQ=255)
                                       ggTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTcaa                                                                                                              >  2:931539/1‑149 (MQ=255)
                                        gTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTTTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTcaac                                                                                                             <  2:1317644/149‑1 (MQ=255)
                                        gTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTcaac                                                                                                             <  2:201844/149‑1 (MQ=255)
                                        gTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTcaac                                                                                                             >  2:1123181/1‑149 (MQ=255)
                                          cAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGTGCGGTGATTAACGGCGACTTcaacaa                                                                                                           <  1:3837241/149‑1 (MQ=255)
                                          cAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTcaacaa                                                                                                           <  1:4751865/149‑1 (MQ=255)
                                            gTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAAtg                                                                                                         <  2:3460287/149‑1 (MQ=255)
                                            gTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAAtg                                                                                                         <  1:2010361/149‑1 (MQ=255)
                                             ttCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAAtgt                                                                                                        <  1:1053515/149‑1 (MQ=255)
                                              tCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAAtgtg                                                                                                       >  1:3488873/1‑149 (MQ=255)
                                                aacaaACCGTGGAGCGTTATCGGCACTGGTCTTTCTGATAAATGCGGTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGTCTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGtgct                                                                                                     <  2:2453342/149‑1 (MQ=255)
                                                aacaaACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGtgct                                                                                                     <  2:748530/149‑1 (MQ=255)
                                                  caaACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGtgctgc                                                                                                   >  1:2560582/1‑149 (MQ=255)
                                                  caaACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGtgctgc                                                                                                   >  1:4581754/1‑149 (MQ=255)
                                                   aaaCCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGcc                                                                                                  <  2:4269249/149‑1 (MQ=255)
                                                    aaCCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCa                                                                                                 <  1:2330780/149‑1 (MQ=255)
                                                     aCCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAg                                                                                                >  1:947398/1‑149 (MQ=255)
                                                     aCCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAg                                                                                                >  2:3855450/1‑149 (MQ=255)
                                                      cagTGGAGCGAAATCGTCACTTGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCTACGACTTTGGCGAGAAAAGTTTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGt                                                                                               <  2:3732969/147‑1 (MQ=255)
                                                       cGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTg                                                                                              <  2:4152705/149‑1 (MQ=255)
                                                           gAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGAtt                                                                                          <  1:3792660/149‑1 (MQ=255)
                                                           gAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGAtt                                                                                          <  1:1617300/149‑1 (MQ=255)
                                                            aGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGAttt                                                                                         >  1:3673226/1‑149 (MQ=255)
                                                               gAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTTTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGt                                                                                      >  1:1372088/1‑149 (MQ=255)
                                                                aaaTCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGtt                                                                                     <  1:303290/149‑1 (MQ=255)
                                                                aaaTCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGtt                                                                                     >  1:4088971/1‑149 (MQ=255)
                                                                 aaTCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTg                                                                                    >  2:3526853/1‑149 (MQ=255)
                                                                  aTCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGa                                                                                   <  1:964045/149‑1 (MQ=255)
                                                                    cGGCACTGGTCTTTCTGATACATGCGTTCTCAGCTACGGCGCAGTCCCGGATTTTGCCAACGACTTGGGCGAGACAAGTTTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATc                                                                                 <  2:591972/149‑1 (MQ=255)
                                                                     ggCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCAGTCCCGGTTATTGCCATCGTCTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATcc                                                                                <  2:2733311/149‑1 (MQ=255)
                                                                       cACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGTCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCgc                                                                              >  2:3205373/1‑149 (MQ=255)
                                                                        aCTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCgca                                                                             <  2:619628/149‑1 (MQ=255)
                                                                         cTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCgcag                                                                            >  1:4256919/1‑149 (MQ=255)
                                                                           ggTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCgcagca                                                                          <  1:649616/149‑1 (MQ=255)
                                                                              cTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTGGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGt                                                                       <  2:3792365/149‑1 (MQ=255)
                                                                               tttCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTg                                                                      <  2:856962/149‑1 (MQ=255)
                                                                                  cTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAg                                                                   <  1:4036925/149‑1 (MQ=255)
                                                                                   tGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAgg                                                                  >  1:1593466/1‑149 (MQ=255)
                                                                                        aaTGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGttt                                                             <  1:1541703/149‑1 (MQ=255)
                                                                                          tGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGt                                                           <  1:3261/149‑1 (MQ=255)
                                                                                          tGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGt                                                           >  1:184583/1‑149 (MQ=255)
                                                                                          tGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGt                                                           >  1:2788440/1‑149 (MQ=255)
                                                                                          tGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGt                                                           <  2:4520806/149‑1 (MQ=255)
                                                                                          tGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGt                                                           >  2:2939348/1‑149 (MQ=255)
                                                                                           gCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTc                                                          >  2:522838/1‑149 (MQ=255)
                                                                                              ttCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGAc                                                       >  2:2304044/1‑149 (MQ=255)
                                                                                                ctcAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCa                                                     <  2:1563077/149‑1 (MQ=255)
                                                                                                     cTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCt                                                <  1:4331072/149‑1 (MQ=255)
                                                                                                       aCGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTgg                                              >  2:2761211/1‑149 (MQ=255)
                                                                                                       aCGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTgg                                              <  2:1400596/149‑1 (MQ=255)
                                                                                                         ggCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTa                                            >  2:3921922/1‑149 (MQ=255)
                                                                                                             cATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGa                                        <  2:4030221/149‑1 (MQ=255)
                                                                                                               ttCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGata                                      <  2:3691969/149‑1 (MQ=255)
                                                                                                                tCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGatat                                     <  1:4830042/149‑1 (MQ=255)
                                                                                                                 cccGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATc                                    >  1:5023275/1‑149 (MQ=255)
                                                                                                                  ccGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATcc                                   <  1:4879523/149‑1 (MQ=255)
                                                                                                                   cGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATccc                                  <  1:4452514/149‑1 (MQ=255)
                                                                                                                      atatTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAAc                               <  1:2931367/149‑1 (MQ=255)
                                                                                                                       tatTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACg                              >  1:2200758/1‑149 (MQ=255)
                                                                                                                         ttGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGAt                            <  1:931539/149‑1 (MQ=255)
                                                                                                                         ttGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGAt                            >  2:5044794/1‑149 (MQ=255)
                                                                                                                         ttGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGAt                            <  2:4221170/149‑1 (MQ=255)
                                                                                                                         ttGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGAt                            >  2:2382875/1‑149 (MQ=255)
                                                                                                                          tGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATc                           <  2:3068239/149‑1 (MQ=255)
                                                                                                                           gCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCa                          >  1:4406388/1‑149 (MQ=255)
                                                                                                                             cAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAgg                        >  2:4663629/1‑149 (MQ=255)
                                                                                                                             cAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAgg                        <  1:2568620/149‑1 (MQ=255)
                                                                                                                             cAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAgg                        <  1:1775393/149‑1 (MQ=255)
                                                                                                                             cAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAgg                        >  1:3213899/1‑149 (MQ=255)
                                                                                                                             cAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGGGGATTTGGTTGATCCGCAGCAGGTGCTGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAgg                        >  2:3379049/1‑149 (MQ=255)
                                                                                                                              aaCGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGt                       <  1:2225086/149‑1 (MQ=255)
                                                                                                                               aCGACTTTGGCGAGAAAAGTTTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTc                      <  2:3177477/149‑1 (MQ=255)
                                                                                                                                  aCTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCggg                   >  2:2954487/1‑149 (MQ=255)
                                                                                                                                      tGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTc               >  2:1633934/1‑149 (MQ=255)
                                                                                                                                       ggCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCa              >  1:4973733/1‑149 (MQ=255)
                                                                                                                                            gaAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGt         <  2:4775297/149‑1 (MQ=255)
                                                                                                                                                aGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGa     <  2:4903274/149‑1 (MQ=255)
                                                                                                                                                aGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGa     <  1:2784520/149‑1 (MQ=255)
                                                                                                                                                 gTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGAt    >  2:1916822/1‑149 (MQ=255)
                                                                                                                                                   ctgctgATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATgg  <  1:5117455/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTCATTAACAACGTGATGATCCCCGACGCCTTAGCCATCGGTCAGTTCAACAAACCGTGGAGCGAAATCGGCACTGGTCTTTCTGATAAATGCGTTCTCAGCTACGGCGCATTCCCGGATATTGCCAACGACTTTGGCGAGAAAAGTCTGCTGATGCCTGGCGGCGCGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGG  >  CP009273/1029517‑1029812

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: