Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,869,269 C→A 6.5% E140D (GAG→GAT cysD ← sulfate adenylyltransferase, subunit 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,869,2690CA 6.5% 749.7 / 24.5 216E140D (GAG→GATcysDsulfate adenylyltransferase, subunit 2
Reads supporting (aligned to +/- strand):  ref base C (85/117);  new base A (3/11);  total (88/128)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

Only 100 of 231 total aligned reads displayed.

GTCCAGTTAGAGAGCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTT  >  CP009273/2869120‑2869418
                                                                                                                                                     |                                                                                                                                                     
gTCCAGTTAGAGAGCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTtc                                                                                                                                                       <  2:1025498/150‑1 (MQ=255)
    aGTTAGAGAGCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCtt                                                                                                                                                   >  2:196059/1‑150 (MQ=255)
    aGTTAGAGAGCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCtt                                                                                                                                                   <  2:5174653/150‑1 (MQ=255)
     gTTAGAGAGCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTc                                                                                                                                                  <  2:3896653/150‑1 (MQ=255)
      ttAGAGAGCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCg                                                                                                                                                 >  1:1692922/1‑150 (MQ=255)
        agagagCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTc                                                                                                                                               <  2:2827333/150‑1 (MQ=255)
        agagagCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTc                                                                                                                                               >  1:2561878/1‑150 (MQ=255)
          agagCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCAc                                                                                                                                             <  2:2799933/150‑1 (MQ=255)
              cGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGgcg                                                                                                                                         <  2:2980854/150‑1 (MQ=255)
              cGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGgcg                                                                                                                                         <  2:3399392/150‑1 (MQ=255)
              cGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGgcg                                                                                                                                         <  2:3572130/150‑1 (MQ=255)
              cGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCGCTTTAGCGCGGGATTTCTCTTCGTCAAGgcg                                                                                                                                         >  1:4879307/1‑150 (MQ=255)
                ggAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGgcgcg                                                                                                                                       <  2:3033340/150‑1 (MQ=255)
                   aGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGcac                                                                                                                                    >  2:3006783/1‑150 (MQ=255)
                   aGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGcac                                                                                                                                    >  1:3468684/1‑150 (MQ=255)
                    gACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGcacc                                                                                                                                   >  1:206074/1‑150 (MQ=255)
                     aCGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGcacca                                                                                                                                  >  2:3448229/1‑150 (MQ=255)
                       gcgGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCaccacc                                                                                                                                >  1:271826/1‑150 (MQ=255)
                         ggATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGa                                                                                                                              >  2:934600/1‑150 (MQ=255)
                         ggATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGa                                                                                                                              >  2:3592951/1‑150 (MQ=255)
                          gATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGaa                                                                                                                             <  1:984377/150‑1 (MQ=255)
                             gCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTATCTTCGTCACGGCGCGCACCACCGAAggc                                                                                                                          <  1:3116703/150‑1 (MQ=255)
                             gCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTATCTTCGTCACGGCGCGCACCACCGAAggc                                                                                                                          <  1:1629868/150‑1 (MQ=255)
                               tttCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg                                                                                                                        <  1:546281/150‑1 (MQ=255)
                                    ccTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCa                                                                                                                   <  2:3181400/150‑1 (MQ=255)
                                     cTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTATCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCaa                                                                                                                  >  1:4591697/1‑150 (MQ=255)
                                      tttGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCaaa                                                                                                                 <  1:2187309/150‑1 (MQ=255)
                                       ttGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCaaaa                                                                                                                >  1:3936145/1‑150 (MQ=255)
                                         gTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAAcc                                                                                                              <  2:4053391/150‑1 (MQ=255)
                                         gTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAAcc                                                                                                              >  1:4093635/1‑150 (MQ=255)
                                             aTTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTAt                                                                                                          <  1:4541790/150‑1 (MQ=255)
                                               ttGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTAttt                                                                                                        >  2:4606855/1‑150 (MQ=255)
                                               ttGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTAttt                                                                                                        <  2:462785/150‑1 (MQ=255)
                                                tGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGCTTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCCTCAAAACCGTATTTg                                                                                                       >  1:1520641/1‑150 (MQ=255)
                                                  cccGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGtt                                                                                                     <  2:4734487/150‑1 (MQ=255)
                                                   ccGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTc                                                                                                    <  2:537204/150‑1 (MQ=255)
                                                    cGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCa                                                                                                   >  2:2349344/1‑150 (MQ=255)
                                                          aGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCt                                                                                             >  2:4906774/1‑150 (MQ=255)
                                                          aGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTATCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCt                                                                                             <  2:1478613/150‑1 (MQ=255)
                                                           gTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTg                                                                                            <  2:4426176/150‑1 (MQ=255)
                                                            ttGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGt                                                                                           <  2:3487801/150‑1 (MQ=255)
                                                             tgtgCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGtt                                                                                          <  2:220581/150‑1 (MQ=255)
                                                               tgCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTc                                                                                        <  1:2965438/150‑1 (MQ=255)
                                                                gCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCa                                                                                       <  2:2139358/150‑1 (MQ=255)
                                                                  cacaGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAgg                                                                                     <  2:2633371/150‑1 (MQ=255)
                                                                  cacaGCTCCGGGCGCTGATTTTTCGGATCCCAGCCATGGAAGCGGTCACGGAAAGAGTAAATTAGCTCTTTATCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAgg                                                                                     >  1:504779/1‑150 (MQ=255)
                                                                   acaGCTCCGGGCGCTGATTTTTCGGATCCCAGAGATGGAAGCGGTCACGGAAATATTAAATTCGCTCTTTATCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGc                                                                                    >  1:816916/1‑150 (MQ=255)
                                                                    caGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGcc                                                                                   <  1:3025826/150‑1 (MQ=255)
                                                                     aGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCGGTTTCAGGCCt                                                                                  >  2:459189/1‑150 (MQ=255)
                                                                         ccGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAg                                                                              <  2:1334864/150‑1 (MQ=255)
                                                                          cGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCGGTTTCAGGCCTTCAGt                                                                             >  1:194983/1‑150 (MQ=255)
                                                                          cGGGCGCTGATTTTTCGGAGCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGt                                                                             <  2:2535277/150‑1 (MQ=255)
                                                                               gcTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCa                                                                        >  2:228476/1‑150 (MQ=255)
                                                                               gcTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCa                                                                        >  2:3543151/1‑150 (MQ=255)
                                                                                cTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCat                                                                       <  1:3311399/150‑1 (MQ=255)
                                                                                 tGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCata                                                                      <  2:843290/150‑1 (MQ=255)
                                                                                  gATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCataa                                                                     >  2:868072/1‑150 (MQ=255)
                                                                                     ttttCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAatat                                                                  <  2:4631403/150‑1 (MQ=255)
                                                                                        tCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCgg                                                               <  2:2548895/150‑1 (MQ=255)
                                                                                        tCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCgg                                                               <  2:1452946/150‑1 (MQ=255)
                                                                                         cGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGt                                                              <  2:1697861/150‑1 (MQ=255)
                                                                                         cGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGt                                                              <  2:1131834/150‑1 (MQ=255)
                                                                                           gATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTAt                                                            <  2:4494729/150‑1 (MQ=255)
                                                                                            aTCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATg                                                           <  1:1269697/150‑1 (MQ=255)
                                                                                             tCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCGGTTTTCATAATATCGGTATGt                                                          >  1:4653721/1‑150 (MQ=255)
                                                                                              cccAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGtt                                                         <  1:1644656/150‑1 (MQ=255)
                                                                                              cccAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTATCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGtt                                                         <  1:1859972/150‑1 (MQ=255)
                                                                                               ccAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGttt                                                        >  2:3142055/1‑150 (MQ=255)
                                                                                                cAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTc                                                       <  2:1598189/150‑1 (MQ=255)
                                                                                                   cGATGGAAGCGGTCACGGAAAGAGTAAAGTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTcgcg                                                    <  2:3298088/150‑1 (MQ=255)
                                                                                                     aTGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCt                                                  <  1:639601/150‑1 (MQ=255)
                                                                                                     aTGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCt                                                  >  2:381565/1‑150 (MQ=255)
                                                                                                      tGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTg                                                 <  2:1624208/150‑1 (MQ=255)
                                                                                                         aaGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCg                                              <  1:1105818/150‑1 (MQ=255)
                                                                                                         aaGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCg                                              <  2:3382010/150‑1 (MQ=255)
                                                                                                         aaGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCg                                              <  1:1068784/150‑1 (MQ=255)
                                                                                                          aGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGt                                             >  2:2425787/1‑150 (MQ=255)
                                                                                                             ggTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCa                                          >  1:4286643/1‑150 (MQ=255)
                                                                                                                  cGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAAt                                     <  2:3977569/150‑1 (MQ=255)
                                                                                                                   ggAAAGAGTAAATTCGCTCTGTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATg                                    <  2:3697213/150‑1 (MQ=255)
                                                                                                                    gAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATgg                                   >  1:5146565/1‑150 (MQ=255)
                                                                                                                      aaGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGAt                                 <  1:4549952/150‑1 (MQ=255)
                                                                                                                       agagTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGAtt                                <  1:3428766/150‑1 (MQ=255)
                                                                                                                       agagTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGAtt                                <  2:3423557/150‑1 (MQ=255)
                                                                                                                       agagTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGAtt                                >  2:2349711/1‑150 (MQ=255)
                                                                                                                       agagTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGAtt                                <  2:4823965/150‑1 (MQ=255)
                                                                                                                          gTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAAt                             >  2:1912213/1‑150 (MQ=255)
                                                                                                                          gTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAAt                             <  2:2563958/150‑1 (MQ=255)
                                                                                                                          gTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCGTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAAt                             >  2:4504918/1‑150 (MQ=255)
                                                                                                                           tAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATc                            <  2:1518470/150‑1 (MQ=255)
                                                                                                                                 cGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATc                      <  1:2101062/150‑1 (MQ=255)
                                                                                                                                 cGCTCTTTAGCGCGGGATTTATCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATc                      <  2:4127009/150‑1 (MQ=255)
                                                                                                                                 cGCTCTTTAGCGCGGGATTTATCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATc                      <  2:3074206/150‑1 (MQ=255)
                                                                                                                                    tctTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGcc                   >  2:649356/1‑150 (MQ=255)
                                                                                                                                          gcgcGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCt             <  2:2930934/150‑1 (MQ=255)
                                                                                                                                           cgcgGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCtt            <  1:1643512/150‑1 (MQ=255)
                                                                                                                                             cgGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTcc          <  2:271826/150‑1 (MQ=255)
                                                                                                                                             cgGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTcc          <  1:3656158/150‑1 (MQ=255)
                                                                                                                                                    tctcTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGtttt   <  2:1234320/150‑1 (MQ=255)
                                                                                                                                                    tctcTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGtttt   <  1:3223638/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCCAGTTAGAGAGCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTT  >  CP009273/2869120‑2869418

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: