Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,902,982 C→A 23.0% S259S (TCG→TCT pstC ← phosphate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,902,9820CA23.0% 818.2 / 209.1 321S259S (TCG→TCTpstCphosphate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (122/125);  new base A (42/32);  total (164/157)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 322 total aligned reads displayed.

ATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTG  >  CP009273/3902833‑3903129
                                                                                                                                                     |                                                                                                                                                   
aTCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCACTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGc                                                                                                                                                     >  2:2081376/1‑150 (MQ=255)
      aaCTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGAGAGgcg                                                                                                                                               >  1:2512432/1‑150 (MQ=255)
       aCTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGgcgc                                                                                                                                              >  2:3794928/1‑150 (MQ=255)
          tCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGAGAGGCGCTGt                                                                                                                                           <  1:3867806/150‑1 (MQ=255)
           cGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTc                                                                                                                                          <  1:475712/150‑1 (MQ=255)
             aTGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGa                                                                                                                                        >  1:3785980/1‑150 (MQ=255)
             aTGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGAGAGGCGCTGTCGa                                                                                                                                        >  2:1371311/1‑150 (MQ=255)
              tGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAg                                                                                                                                       >  2:514273/1‑150 (MQ=255)
              tGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAg                                                                                                                                       <  1:4415171/150‑1 (MQ=255)
                 ggcgAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTg                                                                                                                                    >  2:3991003/1‑150 (MQ=255)
                   cgAGGACGATGAAGGTAATCACACCCAGGATCAGGCCCAGGTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGt                                                                                                                                  <  2:174051/150‑1 (MQ=255)
                     aGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAg                                                                                                                                <  1:1820886/150‑1 (MQ=255)
                     aGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAg                                                                                                                                <  2:691690/150‑1 (MQ=255)
                        aCGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGgtg                                                                                                                             >  1:131093/1‑150 (MQ=255)
                        aCGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGgtg                                                                                                                             >  2:2429393/1‑150 (MQ=255)
                          gATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGtt                                                                                                                           >  2:1463254/1‑150 (MQ=255)
                           aTGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTa                                                                                                                          <  1:3832865/150‑1 (MQ=255)
                            tGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGAGAGGCGCTGTCGAGCTGGTAGGTGTTAc                                                                                                                         >  2:289913/1‑150 (MQ=255)
                              aaGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCg                                                                                                                       >  2:305946/1‑150 (MQ=255)
                                ggTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCgat                                                                                                                     <  1:4267980/150‑1 (MQ=255)
                                 gTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCgatg                                                                                                                    <  2:2066278/150‑1 (MQ=255)
                                       acaAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAg                                                                                                              <  2:2367899/150‑1 (MQ=255)
                                       acaAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAg                                                                                                              <  1:1118333/150‑1 (MQ=255)
                                            cAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGAGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCAc                                                                                                         >  1:2329512/1‑150 (MQ=255)
                                              ggATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGAGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCg                                                                                                       >  2:4497156/1‑150 (MQ=255)
                                               gATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGc                                                                                                      >  1:2111832/1‑150 (MQ=255)
                                                aTCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGcc                                                                                                     >  1:368579/1‑150 (MQ=255)
                                                aTCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGcc                                                                                                     >  1:369437/1‑150 (MQ=255)
                                                   aGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATg                                                                                                  >  2:4059584/1‑150 (MQ=255)
                                                    ggCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGAGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATgg                                                                                                 <  1:3933728/150‑1 (MQ=255)
                                                       ccAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGttt                                                                                              <  1:3311020/150‑1 (MQ=255)
                                                       ccAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGttt                                                                                              <  2:2688505/150‑1 (MQ=255)
                                                        cAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTc                                                                                             <  2:3405696/150‑1 (MQ=255)
                                                        cAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTc                                                                                             <  1:797663/150‑1 (MQ=255)
                                                          gTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCAc                                                                                           <  2:1983206/150‑1 (MQ=255)
                                                           ttCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGAGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCAcc                                                                                          <  2:1975818/150‑1 (MQ=255)
                                                             ccATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGa                                                                                        >  1:1065182/1‑150 (MQ=255)
                                                              cATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAg                                                                                       >  1:2768679/1‑150 (MQ=255)
                                                              cATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGAGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAg                                                                                       <  1:1438126/150‑1 (MQ=255)
                                                                  aGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAgcgcg                                                                                   >  1:2649177/1‑150 (MQ=255)
                                                                       ggCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGcc                                                                              >  2:56966/1‑150 (MQ=255)
                                                                         cAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGAGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGccca                                                                            <  1:4493776/150‑1 (MQ=255)
                                                                          aaCGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGcccag                                                                           >  1:4737196/1‑150 (MQ=255)
                                                                             gTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGcccagccc                                                                        >  1:2231686/1‑150 (MQ=255)
                                                                               gCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGAGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGcccagcccca                                                                      >  1:3487572/1‑150 (MQ=255)
                                                                                 aGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGcccagccccagc                                                                    >  2:1492275/1‑150 (MQ=255)
                                                                                   aCCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCat                                                                  >  2:4430692/1‑150 (MQ=255)
                                                                                   aCCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGAGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCat                                                                  <  1:493237/150‑1 (MQ=255)
                                                                                    ccGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGAGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCatg                                                                 <  1:1142794/150‑1 (MQ=255)
                                                                                     cGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCatga                                                                <  1:3119357/150‑1 (MQ=255)
                                                                                     cGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCatga                                                                <  1:687898/150‑1 (MQ=255)
                                                                                      ggATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCatgat                                                               <  1:2057783/150‑1 (MQ=255)
                                                                                         ttccgcttccgcAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATgcc                                                            <  2:1663688/150‑1 (MQ=255)
                                                                                         ttccgcttccgcAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATgcc                                                            >  2:1241073/1‑150 (MQ=255)
                                                                                         ttccgcttccgcAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGAGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATgcc                                                            <  2:2803231/150‑1 (MQ=255)
                                                                                          tccgcttccgcAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATgccg                                                           <  2:1042918/150‑1 (MQ=255)
                                                                                          tccgcttccgcAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATgccg                                                           <  2:3898043/150‑1 (MQ=255)
                                                                                          tccgcttccgcAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATgccg                                                           <  1:1607823/150‑1 (MQ=255)
                                                                                          tccgcttccgcAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACAGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATgccg                                                           >  1:2989715/1‑150 (MQ=255)
                                                                                           ccgcttccgcAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATgccgc                                                          >  2:1422780/1‑150 (MQ=255)
                                                                                            cgcttccgcAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATgccgcc                                                         <  2:1177904/150‑1 (MQ=255)
                                                                                            cgcttccgcAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATgccgcc                                                         <  1:3560608/150‑1 (MQ=255)
                                                                                                tccgcAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATa                                                     >  2:2653677/1‑150 (MQ=255)
                                                                                                tccgcAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGAGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATa                                                     <  1:289913/150‑1 (MQ=255)
                                                                                                 ccgcAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCAGCCCAGCCCCAGCATGATGCCGCCGATaa                                                    >  2:1096807/1‑150 (MQ=255)
                                                                                                 ccgcAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGAGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATaa                                                    >  1:1269520/1‑150 (MQ=255)
                                                                                                     aaaTTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACAcc                                                <  1:1745585/150‑1 (MQ=255)
                                                                                                        ttcgttcgCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCAtt                                             <  2:3163542/150‑1 (MQ=255)
                                                                                                         tcgttcgCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCAttt                                            <  1:3021880/150‑1 (MQ=255)
                                                                                                          cgttcgCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCAtttt                                           >  1:467490/1‑150 (MQ=255)
                                                                                                          cgttcgCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGAGAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCAtttt                                           >  2:777148/1‑150 (MQ=255)
                                                                                                             tcgCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTgg                                        >  1:3457975/1‑150 (MQ=255)
                                                                                                             tcgCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTgg                                        >  1:83413/1‑150 (MQ=255)
                                                                                                               gCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTg                                      >  1:4086270/1‑150 (MQ=255)
                                                                                                                  aGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAAc                                   >  1:3380917/1‑150 (MQ=255)
                                                                                                                    cgcAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACgg                                 <  1:153196/150‑1 (MQ=255)
                                                                                                                    cgcAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACgg                                 <  1:2712974/150‑1 (MQ=255)
                                                                                                                       agagGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGgaag                              >  1:4162306/1‑150 (MQ=255)
                                                                                                                         agGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGgaagaa                            >  1:605074/1‑150 (MQ=255)
                                                                                                                          ggTGATACTGTTGCCCGGCATATACAGAGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAAc                           >  1:3987849/1‑150 (MQ=255)
                                                                                                                           gTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACg                          <  2:2033174/150‑1 (MQ=255)
                                                                                                                             gATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGAt                        <  1:1482334/150‑1 (MQ=255)
                                                                                                                             gATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGAt                        <  1:3653582/150‑1 (MQ=255)
                                                                                                                             gATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGAt                        <  2:3690307/150‑1 (MQ=255)
                                                                                                                             gATACTGTTGCCCGGCATATACAGAGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGAt                        <  2:1635211/150‑1 (MQ=255)
                                                                                                                               tACTGTTGCCCGGCATATACAGAGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATAc                      <  2:4487723/150‑1 (MQ=255)
                                                                                                                                  tGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGcc                   <  2:4768652/150‑1 (MQ=255)
                                                                                                                                    ttGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAg                 <  2:2234481/150‑1 (MQ=255)
                                                                                                                                       cccGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATa              <  2:1490235/150‑1 (MQ=255)
                                                                                                                                         cGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAAc            <  1:90287/150‑1 (MQ=255)
                                                                                                                                         cGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAAc            >  2:4357345/1‑150 (MQ=255)
                                                                                                                                          ggCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACt           >  2:3950795/1‑150 (MQ=255)
                                                                                                                                          ggCATATACAGAGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACt           <  2:2329512/150‑1 (MQ=255)
                                                                                                                                           gCATATACAGAGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACtt          <  2:3631695/150‑1 (MQ=255)
                                                                                                                                           gCATATACAGAGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACtt          >  2:2921641/1‑150 (MQ=255)
                                                                                                                                            cATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTc         >  2:3580222/1‑150 (MQ=255)
                                                                                                                                            cATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTc         <  2:1699836/150‑1 (MQ=255)
                                                                                                                                             atatACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTcc        >  2:2581792/1‑150 (MQ=255)
                                                                                                                                              tataCAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTccc       >  2:4091933/1‑150 (MQ=255)
                                                                                                                                                taCAGAGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAg     <  2:1034360/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTG  >  CP009273/3902833‑3903129

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: