Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,534,961 C→T 23.8% G26G (GGG→GGA cysW ← sulfate/thiosulfate ABC transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,534,9610CT23.8% 654.8 / 212.9 312G26G (GGG→GGAcysWsulfate/thiosulfate ABC transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (101/136);  new base T (33/41);  total (135/177)
Fisher's exact test for biased strand distribution p-value = 7.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 320 total aligned reads displayed.

CAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTGAGAACAGCAGCAGCAGAGATGCC  >  CP009273/2534813‑2535110
                                                                                                                                                    |                                                                                                                                                     
cAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGACACCAGCATcc                                                                                                                                                      >  2:659565/1‑150 (MQ=255)
 aaTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATccc                                                                                                                                                     <  1:2638217/150‑1 (MQ=255)
 aaTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATccc                                                                                                                                                     <  1:4697735/150‑1 (MQ=255)
     aGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATg                                                                                                                                                 <  2:2059691/150‑1 (MQ=255)
      gcgcGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGc                                                                                                                                                <  2:2657475/150‑1 (MQ=255)
          gATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAAt                                                                                                                                            <  1:3314270/150‑1 (MQ=255)
          gATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAAt                                                                                                                                            <  2:3584563/150‑1 (MQ=255)
          gATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAAt                                                                                                                                            <  2:1265135/150‑1 (MQ=255)
           atcatcACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGCCAATc                                                                                                                                           <  1:3632217/150‑1 (MQ=255)
            tcatcaCCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGCCAATCa                                                                                                                                          >  1:1396076/1‑150 (MQ=255)
            tcatcaCCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCa                                                                                                                                          >  2:3319474/1‑150 (MQ=255)
             catcaCCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGCCAATCAg                                                                                                                                         <  2:1529230/150‑1 (MQ=255)
             catcaCCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAg                                                                                                                                         <  2:3904154/150‑1 (MQ=255)
             catcaCCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAg                                                                                                                                         >  2:1985375/1‑150 (MQ=255)
              atcaCCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGa                                                                                                                                        >  2:1695370/1‑150 (MQ=255)
                 aCCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAAc                                                                                                                                     <  1:511488/150‑1 (MQ=255)
                    gTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCAt                                                                                                                                  >  1:3667391/1‑150 (MQ=255)
                     tCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCAtt                                                                                                                                 <  1:2577279/150‑1 (MQ=255)
                      cAGTCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGCCAATCAGAAACCAttt                                                                                                                                <  2:3639654/150‑1 (MQ=255)
                      cAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCAttt                                                                                                                                >  1:4107648/1‑150 (MQ=255)
                       aGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGCCAATCAGAAACCATTTg                                                                                                                               <  1:3456361/150‑1 (MQ=255)
                       aGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTg                                                                                                                               <  2:1782059/150‑1 (MQ=255)
                          cAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGccc                                                                                                                            >  1:2985772/1‑150 (MQ=255)
                           aGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGcccc                                                                                                                           >  1:4766181/1‑150 (MQ=255)
                               ggCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGtt                                                                                                                       <  2:1452453/150‑1 (MQ=255)
                                  gTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAAt                                                                                                                    <  2:3743619/150‑1 (MQ=255)
                                   tGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATc                                                                                                                   <  2:1036021/150‑1 (MQ=255)
                                      agcaTGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCggg                                                                                                                >  1:4440800/1‑150 (MQ=255)
                                       gcaTGTCCGGATCGGCCAGATTCTGTAAAACAGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGgc                                                                                                               >  1:2670721/1‑150 (MQ=255)
                                        caTGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGgcg                                                                                                              >  2:4381534/1‑150 (MQ=255)
                                        caTGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGAAAAAAAAATAATAAAGATTACAATCAGAAACCATTTGCCCCAGTTAATAGGgcg                                                                                                              >  2:4333218/1‑150 (MQ=255)
                                              cGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTc                                                                                                        >  2:3533644/1‑150 (MQ=255)
                                                 aTCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGACTGCCTGCACGAAGATGTAACTCATCGGCACCAGCAGGATGAACGCCGTAACCAGCATTCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATa                                                                                                     <  1:3796323/150‑1 (MQ=255)
                                                 aTCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATa                                                                                                     <  1:3576337/150‑1 (MQ=255)
                                                 aTCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATa                                                                                                     <  2:59174/150‑1 (MQ=255)
                                                   cGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAAc                                                                                                   <  2:2334517/150‑1 (MQ=255)
                                                   cGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAAc                                                                                                   >  1:3854653/1‑150 (MQ=255)
                                                     gCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGc                                                                                                 >  1:680287/1‑150 (MQ=255)
                                                      ccAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCt                                                                                                <  1:966947/150‑1 (MQ=255)
                                                       cAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCtt                                                                                               <  2:519221/150‑1 (MQ=255)
                                                          aTTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCaa                                                                                            >  1:1957349/1‑150 (MQ=255)
                                                          aTTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCaa                                                                                            <  2:3561422/150‑1 (MQ=255)
                                                            tCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAAtt                                                                                          <  1:488385/150‑1 (MQ=255)
                                                            tCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAAtt                                                                                          <  1:4523539/150‑1 (MQ=255)
                                                                tAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGt                                                                                      >  2:1540104/1‑150 (MQ=255)
                                                                tAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGt                                                                                      <  2:3266959/150‑1 (MQ=255)
                                                                tAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGt                                                                                      >  2:4841496/1‑150 (MQ=255)
                                                                 aaaaCCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTa                                                                                     <  2:3161906/150‑1 (MQ=255)
                                                                  aaaCCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTaa                                                                                    >  1:3984854/1‑150 (MQ=255)
                                                                     ccGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACtt                                                                                 <  2:4690592/150‑1 (MQ=255)
                                                                      cGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTc                                                                                >  1:107924/1‑150 (MQ=255)
                                                                       ggCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTcc                                                                               <  1:1731094/150‑1 (MQ=255)
                                                                         cATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGc                                                                             >  1:802960/1‑150 (MQ=255)
                                                                          aTCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGcc                                                                            <  2:3039292/150‑1 (MQ=255)
                                                                             aGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCAtt                                                                         >  1:1286649/1‑150 (MQ=255)
                                                                                 ccTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATg                                                                     <  1:4333218/150‑1 (MQ=255)
                                                                                 ccTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATg                                                                     <  2:3947980/150‑1 (MQ=255)
                                                                                  cTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGa                                                                    >  1:2259721/1‑150 (MQ=255)
                                                                                   ttGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGAc                                                                   <  1:2764593/150‑1 (MQ=255)
                                                                                   ttGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGAc                                                                   <  1:2134354/150‑1 (MQ=255)
                                                                                     gCTGAATGCCTGCACGAAGTTGTAAATCATCGGCACCAGCTGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCTGTTAATCGGGCGCGCGTCATATCGCTTCAATTGGGTTACTTCCGCCTTTAATGACCt                                                                 <  2:1096994/150‑1 (MQ=255)
                                                                                     gCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCt                                                                 <  1:1337850/150‑1 (MQ=255)
                                                                                      cTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTa                                                                <  1:4598117/150‑1 (MQ=255)
                                                                                        gAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTAcc                                                              >  2:2823616/1‑150 (MQ=255)
                                                                                        gAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTAcc                                                              <  1:3071188/150‑1 (MQ=255)
                                                                                         aaTGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACca                                                             >  1:2922923/1‑150 (MQ=255)
                                                                                         aaTGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACca                                                             >  2:4702305/1‑150 (MQ=255)
                                                                                                 cACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGa                                                     <  2:2957854/150‑1 (MQ=255)
                                                                                                    gAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCa                                                  <  1:1823017/150‑1 (MQ=255)
                                                                                                     aaGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCaa                                                 >  2:3625063/1‑150 (MQ=255)
                                                                                                       gATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAg                                               <  1:2324850/150‑1 (MQ=255)
                                                                                                        aTGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGc                                              <  1:1889678/150‑1 (MQ=255)
                                                                                                         tGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCg                                             >  1:4286242/1‑150 (MQ=255)
                                                                                                           tAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGAc                                           <  1:3250979/150‑1 (MQ=255)
                                                                                                                catcGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGc                                      <  2:2282146/150‑1 (MQ=255)
                                                                                                                    ggCACCAGCAGGATGAACGCCGAAACCAGCATTCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCagag                                  <  2:3918407/150‑1 (MQ=255)
                                                                                                                     gCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAgt                                 >  1:4838364/1‑150 (MQ=255)
                                                                                                                       aCCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAgtgt                               <  2:215679/150‑1 (MQ=255)
                                                                                                                         cagcagGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTa                             >  2:4490516/1‑150 (MQ=255)
                                                                                                                          agcagGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTaa                            >  1:4389379/1‑150 (MQ=255)
                                                                                                                           gcagGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAAt                           <  2:1990587/150‑1 (MQ=255)
                                                                                                                             agGATGAACGCCGAAACCAGCATTCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTg                         >  2:4201774/1‑150 (MQ=255)
                                                                                                                              ggATGAACGCCGAAACCAGCATCCCTATGCCAATCAGAAACCATTTGCCCCATTTAATCGGGCGCGCGCAATAACGCTTCAATTGGGTAACTACCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTga                        <  1:2458865/150‑1 (MQ=255)
                                                                                                                              ggATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTga                        >  1:1270422/1‑150 (MQ=255)
                                                                                                                                aTGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTgaga                      >  1:314099/1‑150 (MQ=255)
                                                                                                                                aTGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTgaga                      >  2:1850673/1‑150 (MQ=255)
                                                                                                                                 tGAACGCCGAAACCAGCATCCCGATGCCCATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTCACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAAGTGAGaa                     >  1:314998/1‑150 (MQ=255)
                                                                                                                                   aaCGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTGAGACca                   >  2:2024389/1‑150 (MQ=255)
                                                                                                                                     cGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTGAGAAcagc                 <  2:1286649/150‑1 (MQ=255)
                                                                                                                                        cGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTGAGAAcagcagc              >  2:3835563/1‑150 (MQ=255)
                                                                                                                                        cGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTGAGAAcagcagc              >  1:545113/1‑150 (MQ=255)
                                                                                                                                         gAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTGAGAAcagcagca             <  1:1581511/150‑1 (MQ=255)
                                                                                                                                          aaaCCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTGAGAAcagcagcag            <  1:3492002/150‑1 (MQ=255)
                                                                                                                                           aaCCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTGAGAAcagcagcagc           >  1:283737/1‑150 (MQ=255)
                                                                                                                                           aaCCAGCATCCCGATGCCAATCAGAAACCATTGGCCCCAGTTAATCGGGCGCGCGTCATACCGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTGAGAAcagcagcagc           <  2:2880569/150‑1 (MQ=255)
                                                                                                                                              cAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTGAGAACAGCAGCAGCaga        <  1:1742241/150‑1 (MQ=255)
                                                                                                                                                 cATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTGAGAACAGCAGCAGCAGAGAt     >  1:1460925/1‑150 (MQ=255)
                                                                                                                                                 cATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTGAGAACAGCAGCAGCAGAGAt     <  2:2493445/150‑1 (MQ=255)
                                                                                                                                                  aTCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTGAGAACAGCAGCAGCAGAGATg    <  2:4752631/150‑1 (MQ=255)
                                                                                                                                                    cccGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTGAGAACAGCAGCAGCAGAGATGcc  <  1:1540133/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGAAGATGTAAATCATCGGCACCAGCAGGATGAACGCCGAAACCAGCATCCCGATGCCAATCAGAAACCATTTGCCCCAGTTAATCGGGCGCGCGTCATAACGCTTCAATTGGGTAACTTCCGCCATTAATGACCTACCACACGCCGACCAAAGCGACTTTGCAGAGTGTTAATTGAGAACAGCAGCAGCAGAGATGCC  >  CP009273/2534813‑2535110

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: