Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,178,919 G→T 6.7% intergenic (‑139/+84) thiM ← / ← rcnR hydoxyethylthiazole kinase/transcriptional repressor of rcnA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,178,9190GT 6.7% 484.0 / 4.1 166intergenic (‑139/+84)thiM/rcnRhydoxyethylthiazole kinase/transcriptional repressor of rcnA
Reads supporting (aligned to +/- strand):  ref base G (44/104);  new base T (7/11);  total (51/115)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 176 total aligned reads displayed.

CGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCTTCTTCACGTTTTAGCTC  >  CP009273/2178771‑2179065
                                                                                                                                                     |                                                                                                                                                   
cGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCgatg                                                                                                                                                     <  2:895574/149‑1 (MQ=255)
 gACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCgatga                                                                                                                                                    <  1:1087409/149‑1 (MQ=255)
  aCTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCgatgat                                                                                                                                                   <  2:2116743/149‑1 (MQ=255)
  aCTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCgatgat                                                                                                                                                   <  2:2117160/149‑1 (MQ=255)
  aCTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCgatgat                                                                                                                                                   <  1:201726/149‑1 (MQ=255)
  aCTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTGCgattat                                                                                                                                                   >  1:1804765/1‑149 (MQ=255)
    ttGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGGTTGGATGATaa                                                                                                                                                 >  1:1349406/1‑149 (MQ=255)
      gCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTCCGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAAt                                                                                                                                               <  2:57614/149‑1 (MQ=255)
       cATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTTCGATGATAAATc                                                                                                                                              >  2:4367465/1‑149 (MQ=255)
         tAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTTGGATGATaaaaaac                                                                                                                                            >  1:4238155/1‑145 (MQ=255)
          aGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCa                                                                                                                                           <  2:814807/149‑1 (MQ=255)
          aGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCa                                                                                                                                           <  2:814379/149‑1 (MQ=255)
           gTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTGCGATGATAAATAGCag                                                                                                                                          >  2:1221992/1‑149 (MQ=255)
           gTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGGTGCGGGTTTAAATCGCag                                                                                                                                          >  1:2767801/1‑149 (MQ=255)
            tttGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCaga                                                                                                                                         <  2:4253560/149‑1 (MQ=255)
            tttGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCaga                                                                                                                                         <  2:4253554/149‑1 (MQ=255)
            tttGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACGTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCaga                                                                                                                                         <  2:3098030/149‑1 (MQ=255)
               gCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAgagga                                                                                                                                     <  1:3218453/149‑1 (MQ=255)
                 tCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTGCGCTGATAAATAGCAGaggagg                                                                                                                                   >  2:262079/1‑149 (MQ=255)
                  ccTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAggagga                                                                                                                                  <  1:1097128/149‑1 (MQ=255)
                   cTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGAt                                                                                                                                 <  1:1863761/149‑1 (MQ=255)
                   cTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGAt                                                                                                                                 <  2:1722046/149‑1 (MQ=255)
                     gCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATgg                                                                                                                               <  2:2939631/149‑1 (MQ=255)
                      ccATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGt                                                                                                                              <  2:124382/149‑1 (MQ=255)
                         tAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATAGCAGAGGAGGATGGTAAt                                                                                                                           >  2:3237079/1‑149 (MQ=255)
                         tAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGGGGGTGCGCTGATAAAAACCCGAGGAGGGTGGGAAt                                                                                                                           >  2:793360/1‑149 (MQ=255)
                           aCGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGGTGCGATGATAAAACGCAGAGGGGGTTGGTAATGt                                                                                                                         >  2:2619615/1‑149 (MQ=255)
                              tGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCa                                                                                                                      <  1:3756635/149‑1 (MQ=255)
                               gaagaaGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTTGGTTCGGTGATAAAAAGCAGATGAGGGTGGTAAtttgcac                                                                                                                     >  1:3683306/1‑143 (MQ=255)
                               gaagaaGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAg                                                                                                                     <  1:3169366/149‑1 (MQ=255)
                                aagaagCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAgc                                                                                                                    <  2:3395631/149‑1 (MQ=255)
                                 agaagCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTGCGATGTTAAAAAGCAGAGGAGGCTGGTAATTTACAgcg                                                                                                                   >  2:2929566/1‑149 (MQ=255)
                                        aTGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTGGGGTGCGATGATAAAAAGCAGAGGAGGAGGGTAATGTCCAgcgcccccgg                                                                                                            >  2:2530587/1‑143 (MQ=255)
                                        aTGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTTCGATGATAAATCGCAGAGGAGGATGGTAATTTCCAGCGCACGCGt                                                                                                            >  2:3571511/1‑149 (MQ=255)
                                           aCCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGGGGTTTGGTTTATAAAAAGCAGATGAGGATGGTAATTTCCAGCCCAccccgtgt                                                                                                         >  1:2074558/1‑142 (MQ=255)
                                           aCCTGGTGGTCCGTGACGTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGt                                                                                                         <  1:2295483/149‑1 (MQ=255)
                                              tggtggTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTaaa                                                                                                      <  2:4165709/149‑1 (MQ=255)
                                               ggtggtCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAAc                                                                                                     <  1:2522781/149‑1 (MQ=255)
                                               ggtggtCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCAGAGTCGTTTGGTTGCGATGTTAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAAc                                                                                                     <  1:75087/149‑1 (MQ=255)
                                                 tggtCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGa                                                                                                   <  2:2059981/149‑1 (MQ=255)
                                                 tggtCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGGTGGGCTTTTAAAACTCAGAGGAGGATGGTTAATTGCCCCGCACCCGTTGTTAAAGa                                                                                                   >  1:1923241/1‑149 (MQ=255)
                                                    tCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGa                                                                                                <  1:1979394/149‑1 (MQ=255)
                                                    tCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGa                                                                                                <  1:1979362/149‑1 (MQ=255)
                                                      cGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGGTGCGCTTATTAAAATCAGAGGAGGTTGGTTATGTCCCGCGCACGCCGTGTAAAAGAGGaat                                                                                              >  1:1316064/1‑148 (MQ=255)
                                                        tGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGaaaac                                                                                            <  1:3962141/149‑2 (MQ=255)
                                                            ttCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGc                                                                                        <  2:4099471/149‑1 (MQ=255)
                                                            ttCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTTCGATGTTAAATAGCAGAGGAGGATGGTAATGTCCCCCGCACGCGTTGTAAACGAGGGAAAACGGc                                                                                        >  2:1675182/1‑149 (MQ=255)
                                                             tCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTGCGATGTTAAAACGCAGAGGAGGCTTGTAAATTACCGCGCACGCGTTGTAAACGAGGAAAAACGGCa                                                                                       >  2:2419797/1‑149 (MQ=255)
                                                                 tACGCTGGCAGTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTACTGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGttt                                                                                   <  2:2107124/149‑1 (MQ=255)
                                                                      tGGCATTATCCAGATCTGGTGATACGGGTATTTCTCCGCCTTCACGCAGATGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGTGGTTGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACaa                                                                              <  1:632477/149‑1 (MQ=255)
                                                                        gCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATc                                                                            <  1:4301250/149‑1 (MQ=255)
                                                                          aTTATCCAGATCCGGTGATACGGGTATTTCTCCGCCTTCACGCAGACGGGCACCCCGTGTCTTTTGGTTGCGTTTATAAATCGCAGTGGTGGATGGTAATGTCCAGCGCACGCGTTGTAAACGCGTAAAAACGGCAGTTTTACAATcgc                                                                          <  1:3717992/149‑1 (MQ=255)
                                                                             aTCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGtt                                                                       <  2:2595886/149‑1 (MQ=255)
                                                                               ccAGATCTGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGACCAGCGCACGCGTTGTAAACGTGCAAACACGGCAGTTTTACAATCGCGTtaa                                                                     <  1:429751/149‑1 (MQ=255)
                                                                                cAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTtaat                                                                    <  2:3051977/149‑1 (MQ=255)
                                                                                 aGATCAGGTGATACGGTTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTtaata                                                                   <  1:1695502/149‑1 (MQ=255)
                                                                                 aGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTtaata                                                                   <  1:3913237/149‑1 (MQ=255)
                                                                                   aTCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATaaa                                                                 <  2:1970341/149‑1 (MQ=255)
                                                                                     cAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAatt                                                               <  1:258319/149‑1 (MQ=255)
                                                                                       ggTGATACGGGTATTTCTCAGCCTTCACGCAGACGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGTTGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAattat                                                             <  1:1329506/149‑1 (MQ=255)
                                                                                          gATACGGGTATTTCTCAGCCTTCACGCAGAAGGGAACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGTTAATGTCCAGCGCACGCTTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTg                                                          <  2:2294129/149‑1 (MQ=255)
                                                                                            tACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGat                                                        <  2:933328/149‑1 (MQ=255)
                                                                                              cGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGatat                                                      <  2:1959423/149‑1 (MQ=255)
                                                                                                gttATTTCTCAGCCTTCACGCAGAAGGGCACCCCGTGTCGTTTGGTTGCGATTATAAATCGCAGTGGTGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGatatat                                                    <  1:1517288/147‑1 (MQ=255)
                                                                                                ggTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGatatat                                                    <  2:1799568/149‑1 (MQ=255)
                                                                                                  tATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTTGGATGTTAAAAAGCAGAGGTGGGTTGGTAATTCACGCCCCCGCGGTGTTAAACAGGGAAAACAGCAGATTTAAAAAAGCGTTAAAAAATTAATTTttataatt                                                  >  1:1788044/1‑146 (MQ=255)
                                                                                                    tttCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGaac                                                <  1:3949171/149‑2 (MQ=255)
                                                                                                    tttCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTTGGATGATAAATAGCAGAGGAGGATGGTAAATTCAACCGCCCGCGTTGTAAACGAGGAAAAAAGGCAGTTTTTAAATCGCGGTAATAAAATATTTTATATATTaaa                                                >  2:1241262/1‑148 (MQ=255)
                                                                                                      tctcAGCCTTCACGCATAAGGGCACCCCTTGTGGTTTGGTTGCTATTATAAATCGCTGTGGTGGTTGTTAATGTCCTGCGCACGCGTTTTAAACGAGTAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATcc                                              <  1:2974683/149‑1 (MQ=255)
                                                                                                      tctcAGCCTTCACGCAGATGGGCACCCCGAGTGTTTTGGTTGCTATTATAAATCGCTGAGGTGTATGGTAATGTCCAGCGCACGGGTTGTAAACTAGTAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATcc                                              <  1:4113653/149‑1 (MQ=255)
                                                                                                      tctcAGCCTTCACGCAGAAGGGCACCCCGCGTCTTTTGGTTGCGATGATAAATCGCTGTGGTGGATGGTAATGTCCAGCGCACGCGTTGTAAACGCGTAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATcc                                              <  1:2639113/149‑1 (MQ=255)
                                                                                                      tctcAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCTATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGTGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATcc                                              <  2:1349406/149‑1 (MQ=255)
                                                                                                      tctcAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATcc                                              <  2:3742029/149‑1 (MQ=255)
                                                                                                         cAGCCTTCACTCAGAAGGGCACCCCGAGTCTTTTGGTTGCTATTATAAATCGAAGTGGAGGATGGTATTGTCCAGCGCACGCTTTTTAAACTAGCAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGc                                           <  1:3358571/149‑1 (MQ=255)
                                                                                                          aGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCa                                          <  1:999945/149‑1 (MQ=255)
                                                                                                              ttCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGTGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCtt                                      <  1:3571511/149‑1 (MQ=255)
                                                                                                                 aCGCAGAAGGGCACCCCGAGTCGTGTGGTTTCGATTATAAAAAGCAGAGGAGGATGGTAATGTCCACCGCACCCGGTTTAAAAGAGGAAAAAACGCAGTTTTAAAAAAGCGGTAAAAAAAAATTTGATATATTAATCCACCACCATCac                                   >  2:1581626/1‑148 (MQ=255)
                                                                                                                 aCGCAGAAGGGCACCCCGAGTCGTGTGGGTGCGATGATAAAAAGCAGAGGGGGGTGGTTATGTACACCGCACCCGGTGTAAAAAAGGAAAAAACGCAGTTTTTCAATAGCCTTAATAAATTATTTTATATTTGAATACACCACCATCaa                                   >  2:3801860/1‑148 (MQ=255)
                                                                                                                  cGCAGACGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGTTAATGTCCAGCGCACGCGTTTTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGa                                  <  1:3387549/149‑1 (MQ=255)
                                                                                                                   gCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGaa                                 <  2:1324791/149‑1 (MQ=255)
                                                                                                                       aaGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCTGAGGTGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGAc                             <  1:938730/149‑1 (MQ=255)
                                                                                                                         gggCACCCCGAGTCGTTTGGTTGCGATTATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACAt                           <  1:326046/149‑1 (MQ=255)
                                                                                                                            cACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGTATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCa                        <  1:501764/149‑1 (MQ=255)
                                                                                                                                    gtgttttGGTTGCGATTATAAATCGCTGGGGAGGTTGGTATTGTCCAGCGCACGCGTTTTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGAtcttct                <  1:2227649/145‑1 (MQ=255)
                                                                                                                                     tCGTTTGGTTGCGATGATAAATCGCAGAGGAGGTTGGTATTGTCCAGCGCACGCTTTGTAAACGAGTAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATcttctt               <  1:2247089/149‑1 (MQ=255)
                                                                                                                                      cGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCttcttc              <  2:4296649/149‑1 (MQ=255)
                                                                                                                                      cGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAAGGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCttcttc              <  2:1136471/149‑1 (MQ=255)
                                                                                                                                       gTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCTTCTTCa             >  2:2519809/1‑149 (MQ=255)
                                                                                                                                        tttGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCTTCTTCAc            >  1:1553030/1‑149 (MQ=255)
                                                                                                                                           tGGTTGCGATTATAAATGGCTGAGGAGTATGGTAATGTCCAGCGCACGAGTTGTAATCGAGGAAAAACGGCAGTTTTACAATCGCGTTAATTAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCTTCTTCACGt          <  2:2392308/149‑1 (MQ=255)
                                                                                                                                            ggTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCTTCTTCACGtt         >  2:4112023/1‑149 (MQ=255)
                                                                                                                                             gTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCTTCTTCACGttt        <  1:4367465/149‑1 (MQ=255)
                                                                                                                                              ttGCGATTATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCTTCTTCACGtttt       <  2:1250887/149‑1 (MQ=255)
                                                                                                                                              ttGCGATGATAAATCGCAGAGGAGGTTGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCTTCTTCACGtttt       <  1:562206/149‑1 (MQ=255)
                                                                                                                                              ttGCGATGATAAATCGCAGAGGAGGATGTTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCTTCTTCACGtttt       <  1:562204/149‑1 (MQ=255)
                                                                                                                                              ttGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCTTCTTCACGtttt       <  2:3821127/149‑1 (MQ=255)
                                                                                                                                                gCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCTTCTTCACGTTTTAg     >  2:255430/1‑149 (MQ=255)
                                                                                                                                                gCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCTTCTTCACGTTTTAg     >  1:3940974/1‑149 (MQ=255)
                                                                                                                                                 cGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCTTCTTCACGTTTTAGc    <  1:2841301/149‑1 (MQ=255)
                                                                                                                                                  gatgatAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCTTCTTCACGTTTTAGCt   >  1:3420710/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGTAATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCTTCTTCACGTTTTAGCTC  >  CP009273/2178771‑2179065

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: