Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,647,653 C→A 7.4% D214Y (GAT→TAT)  sseB ← rhodanase‑like enzyme, sulfur transfer from thiosulfate

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,647,6530CA 7.4% 1128.0 / 43.7 297D214Y (GAT→TAT) sseBrhodanase‑like enzyme, sulfur transfer from thiosulfate
Reads supporting (aligned to +/- strand):  ref base C (158/116);  new base A (11/12);  total (169/128)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 298 total aligned reads displayed.

GGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATAAAAGCACGCTTCACCGGCTTAAT  >  CP009273/2647505‑2647800
                                                                                                                                                    |                                                                                                                                                   
    gcATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCa                                                                                                                                                 >  1:4922471/1‑149 (MQ=255)
    gcATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATAGCCa                                                                                                                                                 >  2:2713790/1‑149 (MQ=255)
       tatTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGt                                                                                                                                              >  2:571837/1‑149 (MQ=255)
         ttttAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATAGCCAGGTaa                                                                                                                                            <  1:2642278/149‑1 (MQ=255)
           ttAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATg                                                                                                                                          >  2:4829020/1‑149 (MQ=255)
                ttATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCg                                                                                                                                     >  1:1164986/1‑149 (MQ=255)
                ttATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATAGCCAGGTAATGTATCg                                                                                                                                     >  2:202806/1‑149 (MQ=255)
                 tATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCgg                                                                                                                                    >  1:4217134/1‑149 (MQ=255)
                  aTCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGt                                                                                                                                   >  1:453605/1‑149 (MQ=255)
                  aTCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGt                                                                                                                                   <  1:3793096/149‑1 (MQ=255)
                  aTCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGt                                                                                                                                   >  1:581105/1‑149 (MQ=255)
                  aTCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGt                                                                                                                                   >  1:1556608/1‑149 (MQ=255)
                   tCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTc                                                                                                                                  >  1:4633229/1‑149 (MQ=255)
                      gATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCt                                                                                                                               <  1:1487216/149‑1 (MQ=255)
                       aTTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATAGCCAGGTAATGTATCGGTCGCTa                                                                                                                              >  1:3200215/1‑149 (MQ=255)
                        ttCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTac                                                                                                                             >  2:766072/1‑149 (MQ=255)
                        ttCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTac                                                                                                                             >  2:2837031/1‑149 (MQ=255)
                             tttAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTc                                                                                                                        <  2:3799846/149‑1 (MQ=255)
                              ttAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGTTTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTcc                                                                                                                       <  1:1377302/149‑1 (MQ=255)
                              ttAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTcc                                                                                                                       <  2:3428959/149‑1 (MQ=255)
                              ttAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATAGCCAGGTAATGTATCGGTCGCTACACTTcc                                                                                                                       <  1:2068831/149‑1 (MQ=255)
                                 aaGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGt                                                                                                                    >  2:1606333/1‑149 (MQ=255)
                                  aGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTc                                                                                                                   >  2:111327/1‑149 (MQ=255)
                                   gTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCg                                                                                                                  >  2:2317352/1‑149 (MQ=255)
                                     cACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGcc                                                                                                                <  1:4834634/149‑1 (MQ=255)
                                      aCGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATAGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCt                                                                                                               >  1:1000287/1‑149 (MQ=255)
                                        gCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGa                                                                                                             <  2:3469506/149‑1 (MQ=255)
                                          aaaaaaCCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTGAGTGCGGTGTCCCCATTTTTCACCTGACACATATCGATTGGTTCATCGCCAGGTAAGGGCGCGGTCGCTCCCCGGCCCCTCGCCCGaca                                                                                                           >  1:17988/1‑147 (MQ=255)
                                              aaCCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCATTGTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAAtt                                                                                                       <  2:3658712/149‑1 (MQ=255)
                                                ccaccCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATttc                                                                                                     >  1:1458399/1‑149 (MQ=255)
                                                ccaccCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATttc                                                                                                     <  1:3540145/149‑1 (MQ=255)
                                                  accCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATttctt                                                                                                   <  1:3712240/149‑1 (MQ=255)
                                                  accCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATAGGTCGCGACACATCCCGTCGCCAGAATAAattatt                                                                                                   >  1:1767757/1‑149 (MQ=255)
                                                   ccccAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATAGATTGGTTCATCGCCAGGTAATGTATCGGGCGCTACACATCCCGTCGCCTGAATAATttattc                                                                                                  >  1:2198607/1‑149 (MQ=255)
                                                      cAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGat                                                                                               <  1:4667603/149‑1 (MQ=255)
                                                       aGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGata                                                                                              >  2:2748690/1‑149 (MQ=255)
                                                       aGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTACCGTCGCCTGAATAATTTCTTCGata                                                                                              >  2:5554/1‑149 (MQ=255)
                                                        gCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGatat                                                                                             >  2:4978878/1‑149 (MQ=255)
                                                         cGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATc                                                                                            <  2:2184357/149‑1 (MQ=255)
                                                              ttCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCAt                                                                                       <  2:883489/149‑1 (MQ=255)
                                                                 aTAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCgg                                                                                    >  2:2229583/1‑149 (MQ=255)
                                                                  tAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGc                                                                                   >  1:4487568/1‑149 (MQ=255)
                                                                  tAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATAGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGc                                                                                   <  1:3676400/149‑1 (MQ=255)
                                                                   aGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCt                                                                                  <  1:3419206/149‑1 (MQ=255)
                                                                        ggCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAAt                                                                             <  2:2705183/149‑1 (MQ=255)
                                                                           gcAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGcc                                                                          <  2:1212710/149‑1 (MQ=255)
                                                                              ctatGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAAt                                                                       >  2:4101775/2‑149 (MQ=255)
                                                                                atGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATaa                                                                     <  1:3224951/149‑1 (MQ=255)
                                                                                   ttCGGTAATAAAGTGGCTCATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTa                                                                  >  1:2333012/1‑149 (MQ=255)
                                                                                   ttCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTa                                                                  >  1:381198/1‑149 (MQ=255)
                                                                                   ttCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTa                                                                  >  2:4780633/1‑149 (MQ=255)
                                                                                     cGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATAGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTaaa                                                                <  1:2481543/149‑1 (MQ=255)
                                                                                       gTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAAtt                                                              <  2:4123311/149‑1 (MQ=255)
                                                                                        taataaAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTa                                                             >  2:5092255/1‑149 (MQ=255)
                                                                                         aataaAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAg                                                            <  1:241648/149‑1 (MQ=255)
                                                                                         aataaAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAg                                                            >  1:406777/1‑149 (MQ=255)
                                                                                         aataaAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAg                                                            <  2:1750275/149‑1 (MQ=255)
                                                                                          ataaAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAgg                                                           >  2:1180397/1‑149 (MQ=255)
                                                                                          ataaAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATGGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAgg                                                           >  2:3855596/1‑149 (MQ=255)
                                                                                           taaAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGc                                                          <  2:2890287/149‑1 (MQ=255)
                                                                                           taaAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATAACCATCGGCTTCAATGCCAATAAGTAAATTAGGc                                                          >  1:2294966/1‑149 (MQ=255)
                                                                                             aaGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCtg                                                        <  2:1381539/149‑1 (MQ=255)
                                                                                             aaGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCtg                                                        <  1:2996161/149‑1 (MQ=255)
                                                                                                tGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGc                                                     >  2:1885751/1‑149 (MQ=255)
                                                                                                 ggCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGcc                                                    >  2:3951094/1‑149 (MQ=255)
                                                                                                 ggCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGcc                                                    <  2:2091869/149‑1 (MQ=255)
                                                                                                   cTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCctc                                                  >  1:1807527/1‑149 (MQ=255)
                                                                                                     aaTTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTATGCTGTGCCTCtt                                                >  1:3088122/1‑149 (MQ=255)
                                                                                                          cTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTtc                                           <  2:4486420/149‑1 (MQ=255)
                                                                                                          cTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTtc                                           <  2:244710/149‑1 (MQ=255)
                                                                                                           tttttCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTtct                                          >  2:3078651/1‑149 (MQ=255)
                                                                                                             tttCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTctct                                        <  2:2602489/149‑1 (MQ=255)
                                                                                                             tttCCCCTTGTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTCTCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTctct                                        <  2:1783100/149‑1 (MQ=255)
                                                                                                              ttCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCtt                                       <  1:5057426/149‑1 (MQ=255)
                                                                                                                 cccTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTaa                                    >  2:3904872/1‑149 (MQ=255)
                                                                                                                  ccTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAAt                                   >  1:5078842/1‑149 (MQ=255)
                                                                                                                      ttttCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGaa                               <  1:1933898/149‑1 (MQ=255)
                                                                                                                        ttCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACa                             >  2:3167283/1‑149 (MQ=255)
                                                                                                                         tCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACaa                            <  2:4119803/149‑1 (MQ=255)
                                                                                                                           aCCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAAt                          >  2:3111736/1‑149 (MQ=255)
                                                                                                                           aCCTGACAAATATCGATTGGTTCATAGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATCAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACACAt                          >  2:4722573/1‑149 (MQ=255)
                                                                                                                            cgagACAAATATCTATTGGTTCATCGACAGGTAATGTATCGGTCGCTACACTTCCAGTCGAATGAATAATTTCTTCTATATCACCATCGGCTTCAATTCCAATAAGTAAATTATGATGTGAATCTTCGTTATCTTTAATTGAACAAATa                         <  2:2392494/146‑1 (MQ=255)
                                                                                                                            ccTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATa                         <  1:1080331/149‑1 (MQ=255)
                                                                                                                              tGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATaaa                       >  2:4799294/1‑149 (MQ=255)
                                                                                                                                aCATATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATAAAAg                     <  1:2837031/149‑1 (MQ=255)
                                                                                                                                aCAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATAAAAg                     <  1:766072/149‑1 (MQ=255)
                                                                                                                                 cAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATAAAAGc                    <  1:5092255/149‑1 (MQ=255)
                                                                                                                                 cAAATATCGATTGGTTCATAGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATAAAAGc                    >  1:1871166/1‑149 (MQ=255)
                                                                                                                                       tCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATAAAAGCACGCtt              >  1:807221/1‑149 (MQ=255)
                                                                                                                                       tCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATAAAAGCACGCtt              >  1:2148327/1‑149 (MQ=255)
                                                                                                                                       tCGATTGGTTCATAGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATAAAAGCACGCtt              >  2:1642420/1‑149 (MQ=255)
                                                                                                                                          aTTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATAAAAGCACGCTTCAc           <  2:4217134/149‑1 (MQ=255)
                                                                                                                                          aTTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATAAAAGCACGCTTCAc           <  1:3962708/149‑1 (MQ=255)
                                                                                                                                                tCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATAAAAGCACGCTTCACCGGCtt     <  1:4151695/149‑1 (MQ=255)
                                                                                                                                                tCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCACCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATAAAAGCACGCTTCACCGGCtt     <  1:261938/149‑1 (MQ=255)
                                                                                                                                                tCATAGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATAAAAGCACGCTTCACCGGCtt     <  1:2439764/149‑1 (MQ=255)
                                                                                                                                                 cATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATAAAAGCACGCTTCACCGGCTTa    >  2:1085348/1‑149 (MQ=255)
                                                                                                                                                  aTCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATAAAAGCACGCTTCACCGGCTTaa   <  2:3223426/149‑1 (MQ=255)
                                                                                                                                                  aTCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATAAAAGCACGCTTCACCGGCTTaa   >  2:913199/1‑149 (MQ=255)
                                                                                                                                                   tCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATAAAAGCACGCTTCACCGGCTTAAt  <  1:2317330/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAATTAGGCTGTGCCTCTTCGTTCTCTTTAATTGAACAAATAAAAGCACGCTTCACCGGCTTAAT  >  CP009273/2647505‑2647800

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: