Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,536,466 A→T 5.5% V141E (GTG→GAG)  cysP ← thiosulfate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,536,4660AT 5.5% 754.5 / 24.6 225V141E (GTG→GAG) cysPthiosulfate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (108/104);  new base T (6/7);  total (114/111)
Fisher's exact test for biased strand distribution p-value = 7.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 229 total aligned reads displayed.

CGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGTGCAGGATTTGTACGTC  >  CP009273/2536317‑2536614
                                                                                                                                                     |                                                                                                                                                    
cGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCa                                                                                                                                                      >  1:2879574/1‑150 (MQ=255)
cGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCa                                                                                                                                                      >  1:2359450/1‑150 (MQ=255)
 gTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCAc                                                                                                                                                     >  2:3331300/1‑150 (MQ=255)
   ttttCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCa                                                                                                                                                   <  1:4442720/150‑1 (MQ=255)
   ttttCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCa                                                                                                                                                   <  2:3936355/150‑1 (MQ=255)
   ttttCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCa                                                                                                                                                   >  2:1791947/1‑150 (MQ=255)
    tttCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAg                                                                                                                                                  <  2:2387013/150‑1 (MQ=255)
    tttCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGAGAGGTATAACGCGCGTTACCCGACGTTTACGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAg                                                                                                                                                  >  2:789065/1‑150 (MQ=255)
    gttCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAg                                                                                                                                                  <  2:3417136/149‑1 (MQ=255)
       cAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTc                                                                                                                                               >  1:4240718/1‑150 (MQ=255)
         ggAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGt                                                                                                                                             >  2:2431246/1‑150 (MQ=255)
         ggAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGt                                                                                                                                             >  2:2242763/1‑150 (MQ=255)
           aaCTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTc                                                                                                                                           >  2:4176563/1‑150 (MQ=255)
           aaCTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTc                                                                                                                                           >  2:663159/1‑150 (MQ=255)
            aCTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTcc                                                                                                                                          <  1:283352/150‑1 (MQ=255)
            aCTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTcc                                                                                                                                          >  2:1025999/1‑150 (MQ=255)
             cTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCa                                                                                                                                         >  1:1623356/1‑150 (MQ=255)
             cTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCa                                                                                                                                         <  1:3452216/150‑1 (MQ=255)
              tGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCaa                                                                                                                                        >  1:865879/1‑150 (MQ=255)
              tGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCaa                                                                                                                                        >  1:4557803/1‑150 (MQ=255)
               gggTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAAt                                                                                                                                       <  2:1918967/150‑1 (MQ=255)
                ggTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATc                                                                                                                                      <  2:1251295/150‑1 (MQ=255)
                  tCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGt                                                                                                                                    >  2:208372/1‑150 (MQ=255)
                    aTAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGAGAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCTCCAGGTCGTTCCAATCGTgg                                                                                                                                  >  2:1168949/1‑150 (MQ=255)
                      aaaCTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGat                                                                                                                                >  2:2975039/1‑150 (MQ=255)
                      aaaCTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGat                                                                                                                                >  1:3704286/1‑150 (MQ=255)
                       aaCTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGata                                                                                                                               >  1:4106857/1‑150 (MQ=255)
                           gTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATAttct                                                                                                                           >  2:2447993/1‑150 (MQ=255)
                            ttCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATAttctt                                                                                                                          <  1:2614799/150‑1 (MQ=255)
                              cGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCg                                                                                                                        <  2:3075083/150‑1 (MQ=255)
                              cGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCg                                                                                                                        <  2:2380538/150‑1 (MQ=255)
                                gTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCggg                                                                                                                      <  2:4412627/150‑1 (MQ=255)
                                  tttGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCTCCAGGTCGTTCCAATCGTGGATATTCTTCGGGtt                                                                                                                    <  2:3628651/150‑1 (MQ=255)
                                  tttGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGtt                                                                                                                    >  1:4490163/1‑150 (MQ=255)
                                    tGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTAc                                                                                                                  >  1:303808/1‑150 (MQ=255)
                                       cTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATCTTCTTCGGGTTACCCt                                                                                                               >  2:1329323/1‑150 (MQ=255)
                                       cTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCt                                                                                                               <  1:3612330/150‑1 (MQ=255)
                                        tttGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTCTAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCtt                                                                                                              >  2:36949/1‑150 (MQ=255)
                                        tttGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCtt                                                                                                              >  1:4109823/1‑150 (MQ=255)
                                               ccaccGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCAcc                                                                                                       >  2:541094/1‑150 (MQ=255)
                                               ccaccGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCAcc                                                                                                       >  2:2521727/1‑150 (MQ=255)
                                               ccaccGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCAcc                                                                                                       >  2:1459916/1‑150 (MQ=255)
                                               ccaccGTCAGCTTTATCCGCTGCGCCCCATGCCGCCACATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCAcc                                                                                                       >  1:2177114/1‑150 (MQ=255)
                                                 accGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAg                                                                                                     >  1:1783562/1‑150 (MQ=255)
                                                   cGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGa                                                                                                   >  2:1612443/1‑150 (MQ=255)
                                                   cGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGa                                                                                                   >  1:1067432/1‑150 (MQ=255)
                                                         cTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCa                                                                                             >  2:1006158/1‑150 (MQ=255)
                                                           ttATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCagg                                                                                           >  1:1256728/1‑148 (MQ=255)
                                                            tATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCAtgg                                                                                          >  1:3299638/1‑150 (MQ=255)
                                                              tCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCAtggtg                                                                                        <  2:1834709/150‑1 (MQ=255)
                                                               ccGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCAtggtgg                                                                                       <  2:1962347/150‑1 (MQ=255)
                                                                 gCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCTCCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAg                                                                                     <  2:10577/150‑1 (MQ=255)
                                                                  cTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGt                                                                                    >  2:1410991/1‑150 (MQ=255)
                                                                     cgcCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGa                                                                                 <  1:1841098/150‑1 (MQ=255)
                                                                      gcCCCATGCCGCCAGATAGGTATAACGCGCGTGACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGaa                                                                                <  1:2860208/150‑1 (MQ=255)
                                                                        cccATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACg                                                                              <  2:2619552/150‑1 (MQ=255)
                                                                        cccATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACg                                                                              >  2:277422/1‑150 (MQ=255)
                                                                          cATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGc                                                                            >  1:4360442/1‑150 (MQ=255)
                                                                          cATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGc                                                                            <  2:398430/150‑1 (MQ=255)
                                                                            tGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGa                                                                          >  1:3580344/1‑150 (MQ=255)
                                                                               cgccAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCt                                                                       >  1:3814923/1‑150 (MQ=255)
                                                                                 ccAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCtat                                                                     <  2:2662105/150‑1 (MQ=255)
                                                                                   aGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCtatta                                                                   <  2:1793861/150‑1 (MQ=255)
                                                                                   aGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCtatta                                                                   >  2:3122133/1‑150 (MQ=255)
                                                                                      tAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTc                                                                >  2:20610/1‑150 (MQ=255)
                                                                                      tAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTc                                                                <  1:4583211/150‑1 (MQ=255)
                                                                                       aGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCg                                                               <  1:1979597/150‑1 (MQ=255)
                                                                                           ataACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAg                                                           >  1:4325959/1‑150 (MQ=255)
                                                                                             aaCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGgc                                                         >  2:8696/1‑150 (MQ=255)
                                                                                             aaCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGgc                                                         >  1:2741887/1‑150 (MQ=255)
                                                                                              aCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGgcg                                                        <  2:2166185/150‑1 (MQ=255)
                                                                                                   cgTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACt                                                   <  1:3287108/150‑1 (MQ=255)
                                                                                                   cgTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACt                                                   <  1:3962553/150‑1 (MQ=255)
                                                                                                    gTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGGTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTg                                                  >  1:3100711/1‑150 (MQ=255)
                                                                                                       aCCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCa                                               >  1:1072517/1‑150 (MQ=255)
                                                                                                          cGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTc                                            >  2:20410/1‑150 (MQ=255)
                                                                                                           gACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCg                                           <  1:3884611/150‑1 (MQ=255)
                                                                                                            aCGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCgg                                          <  2:318943/150‑1 (MQ=255)
                                                                                                               ttttCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCg                                       >  2:419854/1‑150 (MQ=255)
                                                                                                               ttttCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCg                                       <  2:1681080/150‑1 (MQ=255)
                                                                                                                  tcgggttcgggAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGa                                    <  2:630157/150‑1 (MQ=255)
                                                                                                                     ggttcgggAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCa                                 >  2:2347224/1‑150 (MQ=255)
                                                                                                                     ggttcgggAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCa                                 >  2:4540006/1‑150 (MQ=255)
                                                                                                                          gggAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTg                            >  2:2688400/1‑150 (MQ=255)
                                                                                                                          gggAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTg                            >  2:4343885/1‑150 (MQ=255)
                                                                                                                            gAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGcc                          <  1:751194/150‑1 (MQ=255)
                                                                                                                            gAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGcc                          <  1:3145283/150‑1 (MQ=255)
                                                                                                                            gAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGcc                          <  1:3998777/150‑1 (MQ=255)
                                                                                                                                  cAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATc                    <  2:3874858/150‑1 (MQ=255)
                                                                                                                                  cAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATc                    <  1:3323294/150‑1 (MQ=255)
                                                                                                                                  cAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCATGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATc                    <  1:669915/150‑1 (MQ=255)
                                                                                                                                       tCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGTGCa               <  1:3766416/150‑1 (MQ=255)
                                                                                                                                       tCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGTGCa               <  1:1400993/150‑1 (MQ=255)
                                                                                                                                        cACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGTGCAg              <  1:1083486/150‑1 (MQ=255)
                                                                                                                                              ggAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGTGCAGGATTTg        <  2:3742849/150‑1 (MQ=255)
                                                                                                                                              ggAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGTGCAGGATTTg        >  2:1955964/1‑150 (MQ=255)
                                                                                                                                                aGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGTGCAGGATTTGTa      >  2:3541312/1‑150 (MQ=255)
                                                                                                                                                aGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCAGGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGTGCAGGATTTGTa      >  1:1443980/1‑150 (MQ=255)
                                                                                                                                                  cgcACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGTGCAGGATTTGTACg    <  2:1216746/150‑1 (MQ=255)
                                                                                                                                                    cACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGTGCAGGATTTGTACGTc  <  1:1861986/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGTGCAGGATTTGTACGTC  >  CP009273/2536317‑2536614

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: