Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,617,498 G→T 27.1% R204S (CGC→AGC)  dgcZ ← diguanylate cyclase, zinc‑sensing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,617,4980GT27.1% 521.1 / 160.0 243R204S (CGC→AGC) dgcZdiguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base G (101/76);  new base T (26/40);  total (127/116)
Fisher's exact test for biased strand distribution p-value = 2.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 243 total aligned reads displayed.

CCTGCTGTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGATTCAGAGGCTC  >  CP009273/1617353‑1617645
                                                                                                                                                 |                                                                                                                                                   
ccTGCTGTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATa                                                                                                                                                  >  2:2751190/1‑149 (MQ=255)
ccTGCTGTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATa                                                                                                                                                  >  1:3926209/1‑149 (MQ=255)
  tgctgTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATaaa                                                                                                                                                <  2:3374895/149‑1 (MQ=255)
   gctgTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAAc                                                                                                                                               <  1:3107286/149‑1 (MQ=255)
    ctgTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAAcc                                                                                                                                              <  1:1908043/149‑1 (MQ=255)
      gTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGt                                                                                                                                            >  2:3842555/1‑149 (MQ=255)
            gTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTa                                                                                                                                      >  1:2282009/1‑149 (MQ=255)
                  ttGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCAcg                                                                                                                                >  1:1074979/1‑149 (MQ=255)
                      tatGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTc                                                                                                                            >  2:1601666/1‑149 (MQ=255)
                       atGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTcc                                                                                                                           >  2:1598706/1‑149 (MQ=255)
                        tGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCa                                                                                                                          <  2:3404661/149‑1 (MQ=255)
                        tGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCa                                                                                                                          <  2:3583311/149‑1 (MQ=255)
                         gCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCaa                                                                                                                         <  1:265031/149‑1 (MQ=255)
                         gCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCaa                                                                                                                         >  2:2322571/1‑149 (MQ=255)
                          cccTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAAc                                                                                                                        <  1:1281165/149‑1 (MQ=255)
                          cccTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAAc                                                                                                                        <  2:365056/149‑1 (MQ=255)
                            cTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTg                                                                                                                      <  1:3826574/149‑1 (MQ=255)
                            cTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTg                                                                                                                      >  2:2543404/1‑149 (MQ=255)
                             ttCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTgg                                                                                                                     >  1:2676872/1‑149 (MQ=255)
                              tCAGAATGTGTGATGGCATGGTTATCGACTAACTGACAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGc                                                                                                                    >  1:3584841/1‑149 (MQ=255)
                                aGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCta                                                                                                                  >  1:1567144/1‑149 (MQ=255)
                                  aaTGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaag                                                                                                                >  1:2989400/1‑149 (MQ=255)
                                        tgATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGtt                                                                                                          >  1:2787903/1‑149 (MQ=255)
                                        tgATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGtt                                                                                                          >  1:1611080/1‑149 (MQ=255)
                                          aTGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGc                                                                                                        >  2:3332481/1‑149 (MQ=255)
                                              cATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAgg                                                                                                    >  2:2920703/1‑149 (MQ=255)
                                                tGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGt                                                                                                  >  2:2508835/1‑149 (MQ=255)
                                                tGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGt                                                                                                  >  1:2024882/1‑149 (MQ=255)
                                                tGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGt                                                                                                  <  2:448530/149‑1 (MQ=255)
                                                 ggTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTg                                                                                                 <  1:2033960/149‑1 (MQ=255)
                                                  gTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGc                                                                                                <  1:2295047/149‑1 (MQ=255)
                                                        gACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAAta                                                                                          >  1:3183231/1‑149 (MQ=255)
                                                         aCTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAAtac                                                                                         >  1:2494060/1‑149 (MQ=255)
                                                          cTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAAtact                                                                                        >  1:2070695/1‑149 (MQ=255)
                                                            aaCTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAAtactac                                                                                      <  1:1985504/149‑1 (MQ=255)
                                                            aaCTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAAtactac                                                                                      >  1:331206/1‑149 (MQ=255)
                                                              cTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACAt                                                                                    >  1:865295/1‑149 (MQ=255)
                                                               tGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATc                                                                                   >  1:461676/1‑149 (MQ=255)
                                                                 gCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGc                                                                                 >  2:2904447/1‑149 (MQ=255)
                                                                 gCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGc                                                                                 >  2:2741126/1‑149 (MQ=255)
                                                                     aTTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGAt                                                                             >  2:2640492/1‑149 (MQ=255)
                                                                     aTTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGAt                                                                             >  2:914175/1‑149 (MQ=255)
                                                                       tCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTa                                                                           <  1:921114/149‑1 (MQ=255)
                                                                         tGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTaaa                                                                         <  1:141893/149‑1 (MQ=255)
                                                                         tGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTaaa                                                                         >  1:1741309/1‑149 (MQ=255)
                                                                          gACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAAt                                                                        >  1:1203401/1‑149 (MQ=255)
                                                                          gACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAAt                                                                        <  2:818652/149‑1 (MQ=255)
                                                                               cTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCg                                                                   <  1:733818/149‑1 (MQ=255)
                                                                                 gCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTa                                                                 <  1:4086415/149‑1 (MQ=255)
                                                                                   aCGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAgg                                                               <  1:226135/149‑1 (MQ=255)
                                                                                         tGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCAt                                                         >  1:3943535/1‑149 (MQ=255)
                                                                                           cttcttCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTa                                                       <  1:3499470/149‑1 (MQ=255)
                                                                                            ttcttcATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTaa                                                      <  1:3456850/149‑1 (MQ=255)
                                                                                             tcttcATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAAc                                                     >  2:2192248/1‑149 (MQ=255)
                                                                                             tcttcATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAAc                                                     >  2:1924957/1‑149 (MQ=255)
                                                                                              cttcATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAAcc                                                    >  2:3226599/1‑149 (MQ=255)
                                                                                              cttcATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAAcc                                                    <  2:1980327/149‑1 (MQ=255)
                                                                                               ttcATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCa                                                   <  1:556044/149‑1 (MQ=255)
                                                                                                tcatcaTTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCaa                                                  <  1:3486600/149‑1 (MQ=255)
                                                                                                tcatcaTTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCaa                                                  <  1:1238948/149‑1 (MQ=255)
                                                                                                 catcatTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAAt                                                 <  1:1985053/149‑1 (MQ=255)
                                                                                                  atcatTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAAtt                                                >  2:451238/1‑149 (MQ=255)
                                                                                                  atcatTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAAtt                                                <  2:703909/149‑1 (MQ=255)
                                                                                                   tcatTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAAttt                                               <  1:2005333/149‑1 (MQ=255)
                                                                                                     atTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTaa                                             >  2:1482916/1‑149 (MQ=255)
                                                                                                     atTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTaa                                             >  2:4024366/1‑149 (MQ=255)
                                                                                                      ttAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTaaa                                            >  1:1374859/1‑149 (MQ=255)
                                                                                                       tAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAAt                                           <  2:1203401/149‑1 (MQ=255)
                                                                                                       tAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAAt                                           <  2:2919958/149‑1 (MQ=255)
                                                                                                        aGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATc                                          <  1:3253537/149‑1 (MQ=255)
                                                                                                            gCTGTGACAAGAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTc                                      <  2:2115462/149‑1 (MQ=255)
                                                                                                             cTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCa                                     >  2:2645254/1‑149 (MQ=255)
                                                                                                               ttGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAat                                   >  2:1716148/1‑149 (MQ=255)
                                                                                                                tGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAata                                  <  1:2598531/149‑1 (MQ=255)
                                                                                                                   cAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATcc                               <  1:1028393/149‑1 (MQ=255)
                                                                                                                     ataatGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCcaa                             <  1:1657441/149‑1 (MQ=255)
                                                                                                                     ataatGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCcaa                             >  1:440137/1‑149 (MQ=255)
                                                                                                                      taatGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCcaac                            >  2:1347681/1‑149 (MQ=255)
                                                                                                                          gATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAacaaca                        >  2:1306335/1‑149 (MQ=255)
                                                                                                                           aTAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACAt                       >  2:1041437/1‑149 (MQ=255)
                                                                                                                           aTAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACAt                       <  2:3425882/149‑1 (MQ=255)
                                                                                                                            tAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACAtt                      >  2:1662711/1‑149 (MQ=255)
                                                                                                                            tAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACAtt                      >  1:3539292/1‑149 (MQ=255)
                                                                                                                                ttcttcGCCCCCGTAGAGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAAACAATTTAAATCGGTCAATATCCAACAACATtaaat                  <  1:3440581/149‑1 (MQ=255)
                                                                                                                                  cttcGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATtaaataa                >  2:3924025/1‑149 (MQ=255)
                                                                                                                                  cttcGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATtaaataa                <  2:1612555/149‑1 (MQ=255)
                                                                                                                                   ttcGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATtaaataaa               <  2:3400168/149‑1 (MQ=255)
                                                                                                                                   ttcGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATtaaataaa               <  2:1511116/149‑1 (MQ=255)
                                                                                                                                   ttcGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATtaaataaa               <  2:1699456/149‑1 (MQ=255)
                                                                                                                                    tcGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAg              >  1:300135/1‑149 (MQ=255)
                                                                                                                                      gCCCCAGTAGAGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGTGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGAt            <  2:2676872/149‑1 (MQ=255)
                                                                                                                                       cccccGTAGCGATAAACCGTTTCTTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGAtt           <  1:2332443/149‑1 (MQ=255)
                                                                                                                                       cccccGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGAtt           >  2:3109587/1‑149 (MQ=255)
                                                                                                                                       cccccGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGAtt           >  2:1557606/1‑149 (MQ=255)
                                                                                                                                        ccccGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGATTc          >  2:3828436/1‑149 (MQ=255)
                                                                                                                                           cGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTTAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCTATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGATTCaga       <  1:1493846/149‑1 (MQ=255)
                                                                                                                                            gTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGATTCagag      <  2:393716/149‑1 (MQ=255)
                                                                                                                                             tAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGATTCAGAgg     >  2:578408/1‑149 (MQ=255)
                                                                                                                                               gCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGATTCAGAGGCt   >  1:3250067/1‑149 (MQ=255)
                                                                                                                                                cGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGATTCAGAGGCTc  <  1:1685383/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CCTGCTGTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGATTCAGAGGCTC  >  CP009273/1617353‑1617645

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: