Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,869,269 C→A 7.1% 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 7.1% 1130.0 / 45.7 308E140D (GAG→GATcysDsulfate adenylyltransferase, subunit 2
Reads supporting (aligned to +/- strand):  ref base C (147/137);  new base A (8/16);  total (155/153)
Fisher's exact test for biased strand distribution p-value = 9.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 321 total aligned reads displayed.

TCCAGTTAGAGAGCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTT  >  CP009273/2869121‑2869417
                                                                                                                                                    |                                                                                                                                                    
 ccAGTTAGAGAGCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTtct                                                                                                                                                     <  1:3694156/149‑1 (MQ=255)
 ccAGTTAGAGAGCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTtct                                                                                                                                                     <  2:4294427/149‑1 (MQ=255)
   aGTTAGAGAGCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTctct                                                                                                                                                   >  2:4405827/1‑149 (MQ=255)
    gTTAGAGAGCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCtt                                                                                                                                                  >  2:4717069/1‑149 (MQ=255)
      tAGAGAGCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCg                                                                                                                                                <  1:2730127/149‑1 (MQ=255)
       agagagCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGt                                                                                                                                               <  2:2082319/149‑1 (MQ=255)
         agagCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCa                                                                                                                                             >  1:4469090/1‑149 (MQ=255)
          gagCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCAc                                                                                                                                            >  1:4637251/1‑149 (MQ=255)
               ggAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGgcgc                                                                                                                                       >  1:4304798/1‑149 (MQ=255)
               ggAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGGCACGGAAAGAGTAAATTAGCTCGTTATCGCGGGATTTCTCTTAGTCACCgcgc                                                                                                                                       >  1:4799053/1‑149 (MQ=255)
                 aaGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGgcgcgc                                                                                                                                     <  2:2657744/149‑1 (MQ=255)
                  aGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGca                                                                                                                                    >  1:4224779/1‑149 (MQ=255)
                        ggATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCg                                                                                                                              >  2:2301809/1‑149 (MQ=255)
                         gATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGa                                                                                                                             >  2:2216153/1‑149 (MQ=255)
                          aTGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCTCAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGaa                                                                                                                            <  2:427376/149‑1 (MQ=255)
                          aTGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGaa                                                                                                                            <  1:1621421/149‑1 (MQ=255)
                          aTGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGaa                                                                                                                            <  2:2468492/149‑1 (MQ=255)
                           tGCTTTCGCATTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAg                                                                                                                           <  2:4099528/149‑1 (MQ=255)
                            gCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAgg                                                                                                                          <  1:2301809/149‑1 (MQ=255)
                             cTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggc                                                                                                                         <  2:4385923/149‑1 (MQ=255)
                             cTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggc                                                                                                                         >  2:303085/1‑149 (MQ=255)
                                  gCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCAt                                                                                                                    >  2:4052956/1‑149 (MQ=255)
                                   ccTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATc                                                                                                                   <  1:5198047/149‑1 (MQ=255)
                                   ccTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAAATCGCTCTTTAGCGCGGGATTTCTCATCGTCACCGCGCGCACCACCGAAGGCGGCATc                                                                                                                   >  1:4974562/1‑149 (MQ=255)
                                    cTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCa                                                                                                                  >  1:2097701/1‑149 (MQ=255)
                                     ttttatAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTTATTTTTCGGATCAAATCGATGGAAACGGACACGGAAAGAAAAAAATCGCTCATTAGCGCGGGAATTATCTTCGTCACGGCGCGCACCACCGAAGGCCGcaacaa                                                                                                                 >  1:3466097/6‑149 (MQ=255)
                                     tttGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCaa                                                                                                                 >  2:652827/1‑149 (MQ=255)
                                       tGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTATCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCaaaa                                                                                                               >  2:1708204/1‑149 (MQ=255)
                                        gTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTATCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAAc                                                                                                              >  1:1909029/1‑149 (MQ=255)
                                            aTTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTa                                                                                                          <  1:4403062/149‑1 (MQ=255)
                                              ttGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTAtt                                                                                                        <  2:15263/149‑1 (MQ=255)
                                                 cccGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTATCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGt                                                                                                     <  1:4906598/149‑1 (MQ=255)
                                                  ccGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGtt                                                                                                    >  1:550902/1‑149 (MQ=255)
                                                    gTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTAGGATACCAGCGATGGTAGCGGTCACGGAAAGAGTAAATTCGCTCTTTATCGCGGGGTTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCCTATTTGTTCa                                                                                                  >  1:2178703/1‑149 (MQ=255)
                                                     ttGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCCCCACCGAAGGCGGCATCAAAACCGTATTTGTTCAg                                                                                                 >  1:4604181/1‑149 (MQ=255)
                                                         aGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGcc                                                                                             <  2:1027563/149‑1 (MQ=255)
                                                          gTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCt                                                                                            >  2:2845203/1‑149 (MQ=255)
                                                              tgCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGttt                                                                                        <  2:828330/149‑1 (MQ=255)
                                                               gCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTc                                                                                       >  2:3061039/1‑149 (MQ=255)
                                                                 cacaGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAg                                                                                     <  1:3139631/149‑1 (MQ=255)
                                                                 cacaGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAg                                                                                     <  2:3321964/149‑1 (MQ=255)
                                                                  acaGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCGGTTTCAgg                                                                                    >  2:4150600/1‑149 (MQ=255)
                                                                     gCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCt                                                                                 >  2:1999588/1‑149 (MQ=255)
                                                                     gCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTATCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCt                                                                                 <  1:2824885/149‑1 (MQ=255)
                                                                       tCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTc                                                                               <  1:1647003/149‑1 (MQ=255)
                                                                         cGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAAATCGCTCATTAGCGCGGGATTTCTATTCGTCACGGCGCGCACCACCCAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCGGTTTCAGGCCGTCAg                                                                             >  1:1135228/1‑149 (MQ=255)
                                                                          gggCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGt                                                                            >  2:4509703/1‑149 (MQ=255)
                                                                           ggCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGtt                                                                           <  2:1092731/149‑1 (MQ=255)
                                                                            gcgcTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGttt                                                                          <  2:4771009/149‑1 (MQ=255)
                                                                            gcgcTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGttt                                                                          <  1:3954715/149‑1 (MQ=255)
                                                                              gcTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTc                                                                        <  2:4674606/149‑1 (MQ=255)
                                                                                 gATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCata                                                                     <  2:2821772/149‑1 (MQ=255)
                                                                                 gATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCata                                                                     <  1:2065726/149‑1 (MQ=255)
                                                                                     tttCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAatat                                                                 <  1:3191538/149‑1 (MQ=255)
                                                                                       tCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCg                                                               <  1:829721/149‑1 (MQ=255)
                                                                                         ggATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGt                                                             >  1:3271626/1‑149 (MQ=255)
                                                                                          gATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTa                                                            >  1:630295/1‑149 (MQ=255)
                                                                                             cccAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTATCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGt                                                         <  2:3776518/149‑1 (MQ=255)
                                                                                              ccAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGtt                                                        <  2:3511873/149‑1 (MQ=255)
                                                                                              ccAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTATCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGtt                                                        >  2:3978684/1‑149 (MQ=255)
                                                                                               cAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGttt                                                       <  2:4885140/149‑1 (MQ=255)
                                                                                                aGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTc                                                      >  2:1608251/1‑149 (MQ=255)
                                                                                                aGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTc                                                      <  2:737321/149‑1 (MQ=255)
                                                                                                 gCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTcg                                                     >  1:2665360/1‑149 (MQ=255)
                                                                                                        aaGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGcc                                              <  1:1017857/149‑1 (MQ=255)
                                                                                                         aGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCg                                             >  2:880088/1‑149 (MQ=255)
                                                                                                          gCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGt                                            >  1:1746151/1‑149 (MQ=255)
                                                                                                          gCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGt                                            >  1:5121218/1‑149 (MQ=255)
                                                                                                            ggTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGc                                          >  1:1160544/1‑149 (MQ=255)
                                                                                                            ggTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGc                                          >  2:1829528/1‑149 (MQ=255)
                                                                                                            ggTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGc                                          >  1:1246579/1‑149 (MQ=255)
                                                                                                            ggTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTATCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGc                                          <  1:3133916/149‑1 (MQ=255)
                                                                                                             gTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCa                                         >  2:3089588/1‑149 (MQ=255)
                                                                                                              tCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCAc                                        >  2:1256710/1‑149 (MQ=255)
                                                                                                               cACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACg                                       <  2:325281/149‑1 (MQ=255)
                                                                                                                 cGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGaa                                     <  2:122377/149‑1 (MQ=255)
                                                                                                                 cGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTACCGTGCACGaa                                     >  1:2042092/1‑149 (MQ=255)
                                                                                                                  ggAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAAt                                    <  2:550902/149‑1 (MQ=255)
                                                                                                                  ggAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAAt                                    >  2:726591/1‑149 (MQ=255)
                                                                                                                  ggAAAGAGTAAAATCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAAt                                    >  2:288085/1‑149 (MQ=255)
                                                                                                                   gAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTATCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATg                                   >  1:30647/1‑149 (MQ=255)
                                                                                                                    aaaGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATgg                                  <  2:4965166/149‑1 (MQ=255)
                                                                                                                     aaGAGTAAATTCGCTCTTTAGCGCGGGATTTATCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGa                                 <  1:4748496/149‑1 (MQ=255)
                                                                                                                       gagTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGAtt                               <  1:2224562/149‑1 (MQ=255)
                                                                                                                         gTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTaa                             <  1:666706/149‑1 (MQ=255)
                                                                                                                          tAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAAt                            >  1:3040312/1‑149 (MQ=255)
                                                                                                                                cGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCAt                      <  1:3652896/149‑1 (MQ=255)
                                                                                                                                  ctctTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCg                    >  2:5134707/1‑149 (MQ=255)
                                                                                                                                   tctTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGc                   <  2:2770115/149‑1 (MQ=255)
                                                                                                                                    ctTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGcc                  <  1:291423/149‑1 (MQ=255)
                                                                                                                                     tttAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCa                 <  1:4079118/149‑1 (MQ=255)
                                                                                                                                     tttAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCa                 >  2:14627/1‑149 (MQ=255)
                                                                                                                                      ttAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCAc                >  1:1125328/1‑149 (MQ=255)
                                                                                                                                         gcgcGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGcc             <  1:3858560/149‑1 (MQ=255)
                                                                                                                                         gcgcGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTATTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGcc             <  2:2607037/149‑1 (MQ=255)
                                                                                                                                            cgGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTc          <  1:2714233/149‑1 (MQ=255)
                                                                                                                                               gATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCAACAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTGCAGTTTTCAAAATATCGGTATGTTTCGCGCTGCCGTGCACGAAAGGATTAATCCCCATCGCCACGCCTTCCgg       >  2:2035724/1‑149 (MQ=255)
                                                                                                                                                aTTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCggg      >  1:3925196/1‑149 (MQ=255)
                                                                                                                                                aTTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCggg      <  2:3447800/149‑1 (MQ=255)
                                                                                                                                                  ttCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGtt    <  2:4786658/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCCAGTTAGAGAGCGGGAAGACGCGGATGCTTTCGCCTTTGTTAATTTGCCCGTTGTAGTTGTGCCACAGCTCCGGGCGCTGATTTTTCGGATCCCAGCGATGGAAGCGGTCACGGAAAGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGGCATCAAAACCGTATTTGTTCAGCGCCTGTTTCAGGCCTTCAGTTTTCATAATATCGGTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTT  >  CP009273/2869121‑2869417

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