Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,534,363 C→A 7.4% G226C (GGC→TGC)  cysW ← sulfate/thiosulfate ABC transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,534,3630CA 7.4% 596.3 / 17.0 161G226C (GGC→TGC) cysWsulfate/thiosulfate ABC transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (87/62);  new base A (9/3);  total (96/65)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

CACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGTC  >  CP009273/2534218‑2534510
                                                                                                                                                 |                                                                                                                                                   
cACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCgg                                                                                                                                                   <  1:2463258/148‑1 (MQ=255)
 aCTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGa                                                                                                                                                  >  1:1765523/1‑148 (MQ=255)
  cTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGaa                                                                                                                                                 >  2:678659/1‑148 (MQ=255)
       aCATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAaccac                                                                                                                                            >  2:3063168/1‑148 (MQ=255)
          tACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGa                                                                                                                                         <  2:3004226/148‑1 (MQ=255)
              ttttAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACAcc                                                                                                                                     >  2:923716/1‑148 (MQ=255)
                 tAAATACCGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACcgcg                                                                                                                                  <  2:2595769/148‑1 (MQ=255)
                     aaCAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCaa                                                                                                                              >  2:1854409/1‑148 (MQ=255)
                       cAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAAc                                                                                                                            >  1:2887518/1‑148 (MQ=255)
                          ggTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCCGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCg                                                                                                                         >  2:113679/1‑148 (MQ=255)
                          ggTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCg                                                                                                                         >  2:1139645/1‑148 (MQ=255)
                             taaaaccGCCATCCGCGTTAACCGCGCCGCCGGGGTAACGCCGCCGACGTTTTTGTTGTCCTGCACCAGAAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCa                                                                                                                      <  1:3061156/142‑1 (MQ=255)
                             gATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCa                                                                                                                      >  2:2505390/1‑148 (MQ=255)
                             gATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCa                                                                                                                      >  2:2091862/1‑148 (MQ=255)
                             gATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCa                                                                                                                      >  2:2091889/1‑148 (MQ=255)
                             gATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCa                                                                                                                      >  2:2505386/1‑148 (MQ=255)
                             gATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCa                                                                                                                      >  2:426969/1‑148 (MQ=255)
                                aaTCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAAtt                                                                                                                   >  1:98829/1‑148 (MQ=255)
                                  tCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTgc                                                                                                                 >  1:1243655/1‑148 (MQ=255)
                                  tCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTgc                                                                                                                 >  2:3207201/1‑148 (MQ=255)
                                  tCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTgc                                                                                                                 >  1:1319754/1‑148 (MQ=255)
                                   cGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTgcg                                                                                                                >  2:157475/1‑148 (MQ=255)
                                    gCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGAGGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTgcgc                                                                                                               <  2:563384/148‑1 (MQ=255)
                                       aTCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCggg                                                                                                            >  1:1310717/1‑148 (MQ=255)
                                        tCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGc                                                                                                           >  1:2942216/1‑148 (MQ=255)
                                        tCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGc                                                                                                           >  1:1925104/1‑148 (MQ=255)
                                        tCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGc                                                                                                           >  2:671034/1‑148 (MQ=255)
                                           gCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGtt                                                                                                        >  2:1580222/1‑148 (MQ=255)
                                            cGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTg                                                                                                       >  1:1918776/1‑148 (MQ=255)
                                                aaCCGCGCCGCCGCGGTACAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCa                                                                                                   <  2:1932423/148‑1 (MQ=255)
                                                 aCAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCaa                                                                                                  >  2:1957323/1‑148 (MQ=255)
                                                 aCAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCaa                                                                                                  >  2:3170904/1‑148 (MQ=255)
                                                  cAGCTCCGCAGCGGTAAAGGAGACGACGGTGTTGTAGTACTGCTCCAGCAATTCAATCTGTAACGGAAGCGACAGGGTTTCGCCGCGAATCGAGCAGTAAACCACCGACACCTCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAac                                                                                                 >  1:10619/1‑148 (MQ=255)
                                                  cAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAac                                                                                                 >  2:73671/1‑148 (MQ=255)
                                                  cAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAac                                                                                                 <  2:3248901/148‑1 (MQ=255)
                                                   aGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAaca                                                                                                >  1:670933/1‑148 (MQ=255)
                                                   aGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCGGTACCGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAaca                                                                                                <  1:2368655/148‑1 (MQ=255)
                                                     cgccgcAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAAcacc                                                                                              <  2:984009/148‑1 (MQ=255)
                                                      gccgcAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAAcacca                                                                                             <  1:1771700/148‑1 (MQ=255)
                                                          cAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGcc                                                                                         >  1:760774/1‑148 (MQ=255)
                                                              ggTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATaa                                                                                     <  2:2904728/148‑1 (MQ=255)
                                                                 aaaGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAagc                                                                                  >  2:1106890/1‑148 (MQ=255)
                                                                   aGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAagcag                                                                                >  2:2120823/1‑148 (MQ=255)
                                                                    ggAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAagcagc                                                                               <  1:3133597/148‑1 (MQ=255)
                                                                     gAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAgcg                                                                              >  1:2917310/1‑148 (MQ=255)
                                                                     gAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAgcg                                                                              >  1:622685/1‑148 (MQ=255)
                                                                       gCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTAAACACCACGCCATAAAGCAGCGcc                                                                            <  2:1918776/148‑1 (MQ=255)
                                                                        ccGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGccc                                                                           >  2:2786750/1‑148 (MQ=255)
                                                                        ccGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGccc                                                                           >  2:2787463/1‑148 (MQ=255)
                                                                               tgttgtAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGAt                                                                    >  2:2493514/1‑148 (MQ=255)
                                                                                gttgtAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATg                                                                   <  1:896974/148‑1 (MQ=255)
                                                                                 ttgtAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGt                                                                  <  2:3232065/148‑1 (MQ=255)
                                                                                 ttgtAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGt                                                                  <  1:1382365/148‑1 (MQ=255)
                                                                                  tgtAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGtt                                                                 >  1:1620342/1‑148 (MQ=255)
                                                                                  tgtAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGtt                                                                 >  1:3257534/1‑148 (MQ=255)
                                                                                  tgtAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGtt                                                                 >  2:2562628/1‑148 (MQ=255)
                                                                                  tgtAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGtt                                                                 <  1:2642865/148‑1 (MQ=255)
                                                                                   ttAGTCCTGCTCCATCAATTCAATGTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGACGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTc                                                                <  1:2556449/147‑1 (MQ=255)
                                                                                   gtAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTc                                                                >  1:931495/1‑148 (MQ=255)
                                                                                    tAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCg                                                               >  2:2892231/1‑148 (MQ=255)
                                                                                     aGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCgg                                                              <  1:2481256/148‑1 (MQ=255)
                                                                                      gTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGt                                                             <  1:2447692/148‑1 (MQ=255)
                                                                                      gTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGt                                                             <  2:1243695/148‑1 (MQ=255)
                                                                                      gTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGt                                                             <  2:1243655/148‑1 (MQ=255)
                                                                                       tCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTa                                                            >  1:2162238/1‑148 (MQ=255)
                                                                                        ccTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTaa                                                           >  2:1647906/1‑148 (MQ=255)
                                                                                        ccTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCAGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTaa                                                           >  2:3258099/1‑148 (MQ=255)
                                                                                                aGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCAGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGAcgc                                                   >  1:1094017/1‑148 (MQ=255)
                                                                                                aGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCAGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGAcgc                                                   >  1:1094026/1‑148 (MQ=255)
                                                                                                    aTTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACg                                               >  1:412515/1‑148 (MQ=255)
                                                                                                       cAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGaa                                            <  1:1033122/148‑1 (MQ=255)
                                                                                                          tCTGTAACGGCAGCGACAGGGTTTCGCCGCGACTCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACAt                                         <  2:2850940/148‑1 (MQ=255)
                                                                                                          tCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACAt                                         >  1:167862/1‑148 (MQ=255)
                                                                                                          tCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCAGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACAt                                         <  2:2874731/148‑1 (MQ=255)
                                                                                                          gaggTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACAt                                         <  1:2844024/145‑1 (MQ=255)
                                                                                                           cTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATc                                        <  2:2301580/148‑1 (MQ=255)
                                                                                                           cTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATc                                        >  2:715309/1‑148 (MQ=255)
                                                                                                           cTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATc                                        >  2:2650379/1‑148 (MQ=255)
                                                                                                            tGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCt                                       >  1:2435486/1‑148 (MQ=255)
                                                                                                             gTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTg                                      <  1:2296323/148‑1 (MQ=255)
                                                                                                               aaCGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGcc                                    >  1:2813527/1‑148 (MQ=255)
                                                                                                               aaCGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGcc                                    >  1:2813535/1‑148 (MQ=255)
                                                                                                                    cAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCg                               >  2:2829749/1‑148 (MQ=255)
                                                                                                                    cAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCg                               >  2:784836/1‑148 (MQ=255)
                                                                                                                    cAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCg                               >  1:1280622/1‑148 (MQ=255)
                                                                                                                     aGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCgg                              <  2:3257534/148‑1 (MQ=255)
                                                                                                                        gACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGAcg                           <  1:3207201/148‑1 (MQ=255)
                                                                                                                         aCAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGAcgc                          >  2:557686/1‑148 (MQ=255)
                                                                                                                          cAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGAcgcg                         >  2:2711328/1‑148 (MQ=255)
                                                                                                                           aGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACgcgc                        <  1:3046518/148‑1 (MQ=255)
                                                                                                                                     cgcgAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAAt              <  1:923716/148‑1 (MQ=255)
                                                                                                                                       cgAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATcg            <  1:1531482/148‑1 (MQ=255)
                                                                                                                                       cgAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATcg            >  1:1272427/1‑148 (MQ=255)
                                                                                                                                         aaTCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCACGCAAAATcgcc          >  2:1914065/1‑148 (MQ=255)
                                                                                                                                          aTCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATcgccg         >  1:1848552/1‑148 (MQ=255)
                                                                                                                                            cGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCt       >  1:3193487/1‑148 (MQ=255)
                                                                                                                                            cGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCt       <  1:986012/148‑1 (MQ=255)
                                                                                                                                             gAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCtt      <  2:2404034/148‑1 (MQ=255)
                                                                                                                                                ccGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGt   >  1:402503/1‑148 (MQ=255)
                                                                                                                                                ccGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGt   <  2:1310717/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CACTGCAACATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGCAAAATCGCCGCTTCGTC  >  CP009273/2534218‑2534510

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: