Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,617,498 G→T 11.7% 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,4980GT11.7% 770.8 / 50.4 230R204S (CGC→AGC) dgcZdiguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base G (99/104);  new base T (14/13);  total (113/117)
Fisher's exact test for biased strand distribution p-value = 8.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 233 total aligned reads displayed.

CACCTGCTGTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGATTCAGAGGC  >  CP009273/1617351‑1617643
                                                                                                                                                   |                                                                                                                                                 
cacCTGCTGTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCg                                                                                                                                                   <  1:663919/148‑1 (MQ=255)
cacCTGCTGTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCg                                                                                                                                                   <  2:1688438/148‑1 (MQ=255)
 acCTGCTGTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGa                                                                                                                                                  <  2:2808506/148‑1 (MQ=255)
  ccTGCTGTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGAt                                                                                                                                                 >  2:2707218/1‑148 (MQ=255)
   ctgctgTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATa                                                                                                                                                >  2:1128207/1‑148 (MQ=255)
    tgctgTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATaa                                                                                                                                               >  1:3094791/1‑148 (MQ=255)
     gctgTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTATCTGGCAATTTATGACACGTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATaaa                                                                                                                                              <  1:2878404/148‑1 (MQ=255)
     gctgTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATaaa                                                                                                                                              <  1:24423/148‑1 (MQ=255)
        gTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCg                                                                                                                                           >  2:2355875/1‑148 (MQ=255)
          cACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGtt                                                                                                                                         >  2:150535/1‑148 (MQ=255)
           aCGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGttt                                                                                                                                        <  2:780252/148‑1 (MQ=255)
              gTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGt                                                                                                                                     >  1:3464270/1‑148 (MQ=255)
              gTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGt                                                                                                                                     >  1:2560541/1‑148 (MQ=255)
              gTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGt                                                                                                                                     <  1:2308442/148‑1 (MQ=255)
               tAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTa                                                                                                                                    >  2:45914/1‑148 (MQ=255)
                aaTGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTaa                                                                                                                                   <  1:758708/148‑1 (MQ=255)
                  tGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATc                                                                                                                                 >  1:1470933/1‑148 (MQ=255)
                  tGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATc                                                                                                                                 <  1:2316622/148‑1 (MQ=255)
                  tGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATc                                                                                                                                 <  2:3212302/148‑1 (MQ=255)
                   gTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCa                                                                                                                                <  1:949047/148‑1 (MQ=255)
                   gTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCa                                                                                                                                >  1:442792/1‑148 (MQ=255)
                         atGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTc                                                                                                                          >  1:524416/1‑148 (MQ=255)
                           gCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCa                                                                                                                        >  1:946234/1‑148 (MQ=255)
                               ttCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTg                                                                                                                    >  2:86854/1‑148 (MQ=255)
                                       gtgtgATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaagta                                                                                                            <  1:3344436/148‑1 (MQ=255)
                                        tgtgATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaagtaa                                                                                                           >  1:3582216/1‑148 (MQ=255)
                                        tgtgATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaagtaa                                                                                                           <  2:346973/148‑1 (MQ=255)
                                         gtgATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaagtaag                                                                                                          <  1:83432/148‑1 (MQ=255)
                                            aTGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTg                                                                                                       <  2:1131184/148‑1 (MQ=255)
                                            aTGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTg                                                                                                       >  2:3618948/1‑148 (MQ=255)
                                            aTGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTg                                                                                                       >  2:2524806/1‑148 (MQ=255)
                                             tGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCCTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGc                                                                                                      <  1:892112/148‑1 (MQ=255)
                                               gCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCa                                                                                                    <  2:3134857/148‑1 (MQ=255)
                                                   ggTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGt                                                                                                >  1:1865698/1‑148 (MQ=255)
                                                      tATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCg                                                                                             <  1:2638170/148‑1 (MQ=255)
                                                      tATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCg                                                                                             >  1:1356418/1‑148 (MQ=255)
                                                       aTCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGt                                                                                            >  2:1639517/1‑148 (MQ=255)
                                                        tCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTa                                                                                           >  2:2368535/1‑148 (MQ=255)
                                                        tCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTa                                                                                           <  1:3680379/148‑1 (MQ=255)
                                                         cGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTaa                                                                                          >  1:2113427/1‑148 (MQ=255)
                                                             tAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAAtact                                                                                      >  2:2790358/1‑148 (MQ=255)
                                                              aaCTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAAtacta                                                                                     >  1:1778067/1‑148 (MQ=255)
                                                               aCTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTAGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAAtactac                                                                                    <  1:2840782/148‑1 (MQ=255)
                                                                 tGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACAt                                                                                  <  2:1073421/148‑1 (MQ=255)
                                                                  ggCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATc                                                                                 >  2:708920/1‑148 (MQ=255)
                                                                   gCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCg                                                                                <  2:3652108/148‑1 (MQ=255)
                                                                      aaTTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCg                                                                             <  2:734229/148‑1 (MQ=255)
                                                                      aaTTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCg                                                                             <  2:1709875/148‑1 (MQ=255)
                                                                       aTTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGa                                                                            <  1:510548/148‑1 (MQ=255)
                                                                       aTTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGa                                                                            <  2:2160187/148‑1 (MQ=255)
                                                                       aTTCTGACACCTGCACGACATGCTTCTTCATAATTAGCCGCTATGACAATAATGATAAATTCATCGCCCCCGTAGCGATAAACAGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGa                                                                            >  1:2688682/1‑148 (MQ=255)
                                                                           tGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTaa                                                                        <  2:293906/148‑1 (MQ=255)
                                                                            gACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTaaa                                                                       >  2:804307/1‑148 (MQ=255)
                                                                               acCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGc                                                                    <  1:150535/148‑1 (MQ=255)
                                                                               acCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGc                                                                    >  2:1329605/1‑148 (MQ=255)
                                                                                   gctcGACATGATTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGt                                                                <  1:2096256/145‑1 (MQ=255)
                                                                                   gCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGt                                                                <  1:78546/148‑1 (MQ=255)
                                                                                        aCATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTa                                                           <  1:1475528/148‑1 (MQ=255)
                                                                                        aCATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTa                                                           <  2:204766/148‑1 (MQ=255)
                                                                                         cATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTAt                                                          >  1:3014950/1‑148 (MQ=255)
                                                                                         cATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTAt                                                          <  2:3069239/148‑1 (MQ=255)
                                                                                         cATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTAt                                                          >  1:2199308/1‑148 (MQ=255)
                                                                                         cATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTAt                                                          >  1:2199337/1‑148 (MQ=255)
                                                                                           tGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCa                                                        <  1:3542431/148‑1 (MQ=255)
                                                                                            gCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCAt                                                       >  2:406043/1‑148 (MQ=255)
                                                                                            gCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCAt                                                       >  2:2329401/1‑148 (MQ=255)
                                                                                              ttcttcATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTa                                                     <  2:3794682/148‑1 (MQ=255)
                                                                                               tcttcATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTaa                                                    >  2:2843581/1‑148 (MQ=255)
                                                                                               tcttcATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTaa                                                    <  2:1283071/148‑1 (MQ=255)
                                                                                                cttcATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAAc                                                   >  1:1662091/1‑148 (MQ=255)
                                                                                                     tcatTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAAtt                                              >  1:1845015/1‑148 (MQ=255)
                                                                                                        ttAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTaa                                           <  1:1648988/148‑1 (MQ=255)
                                                                                                        ttAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTaa                                           <  2:128959/148‑1 (MQ=255)
                                                                                                        ttAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTaa                                           <  1:2926719/148‑1 (MQ=255)
                                                                                                         tAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTaaa                                          >  2:1296913/1‑148 (MQ=255)
                                                                                                          aGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAAt                                         <  2:2560541/148‑1 (MQ=255)
                                                                                                          aGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAAt                                         <  2:1565400/148‑1 (MQ=255)
                                                                                                            ccGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCg                                       <  2:2482823/148‑1 (MQ=255)
                                                                                                            ccGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCg                                       >  1:497759/1‑148 (MQ=255)
                                                                                                             cGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCgg                                      <  1:1985731/148‑1 (MQ=255)
                                                                                                              gCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGt                                     <  1:641664/148‑1 (MQ=255)
                                                                                                                tttGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCa                                   <  1:2355875/148‑1 (MQ=255)
                                                                                                                    aCAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAatat                               >  2:2000000/1‑148 (MQ=255)
                                                                                                                    aCAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAatat                               >  1:3637642/1‑148 (MQ=255)
                                                                                                                       ataatGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCca                            >  2:1956480/1‑148 (MQ=255)
                                                                                                                        taatGATAAATTCTTCGCCCCCGTAGCGAAAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCAGTAGGTATCATTAACCAATTTAAATCGGTCAATAACcaa                           >  2:539269/1‑148 (MQ=255)
                                                                                                                               aaaTTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACAtt                    >  2:760646/1‑148 (MQ=255)
                                                                                                                                aaTTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATta                   <  2:442792/148‑1 (MQ=255)
                                                                                                                                 aTTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATtaa                  >  1:10231/1‑148 (MQ=255)
                                                                                                                                   tcttcGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATtaaat                >  1:187774/1‑148 (MQ=255)
                                                                                                                                     ttcGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATtaaataa              <  1:3189048/148‑1 (MQ=255)
                                                                                                                                      tcGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATtaaataaa             >  2:1804179/1‑148 (MQ=255)
                                                                                                                                        gCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGa           <  1:708262/148‑1 (MQ=255)
                                                                                                                                         cccccGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGAt          <  1:1237429/148‑1 (MQ=255)
                                                                                                                                         cccccGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGAt          >  1:1091137/1‑148 (MQ=255)
                                                                                                                                           cccGTAGCTATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGATTc        <  2:3582216/148‑1 (MQ=255)
                                                                                                                                           cccGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGATTc        <  1:3246632/148‑1 (MQ=255)
                                                                                                                                             cGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGATTCag      >  2:2171270/1‑148 (MQ=255)
                                                                                                                                                aGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGATTCAGAgg   >  2:689708/1‑148 (MQ=255)
                                                                                                                                                 gCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGATTCAGAGGc  >  1:3520567/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CACCTGCTGTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCGTAATACTACATCGCCGATTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGATTCAGAGGC  >  CP009273/1617351‑1617643

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: