Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,305,116 C→A 92.5% intergenic (‑103/+9) apbE ← / ← ompC putative thiamine‑synthetic flavin transferase lipoprotein/outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,305,1160CA92.5% 731.3 / 48.1 240intergenic (‑103/+9)apbE/ompCputative thiamine‑synthetic flavin transferase lipoprotein/outer membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base C (12/6);  new base A (113/109);  total (125/115)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

Only 100 of 242 total aligned reads displayed.

GAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAG  >  CP009273/2304969‑2305261
                                                                                                                                                   |                                                                                                                                                 
gagCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGATTATTGCCATCCGTGATTATCCTCATGCGAACGGTCGCAAGAGAACACCAAAAAACTGAGTTTGTACGCTGAAAAAAATCAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATaaa                                                                                                                                                 <  1:736412/150‑1 (MQ=255)
   cgcgGCAGCCAGCAGTGCCACGCGGGTACAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGATTTTGTACGCTTAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCg                                                                                                                                              <  2:3795990/150‑1 (MQ=255)
      ggCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCgaga                                                                                                                                           >  2:4007403/1‑150 (MQ=255)
      ggCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCgaga                                                                                                                                           >  2:1054243/1‑150 (MQ=255)
      ggCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCgaga                                                                                                                                           >  2:1672109/1‑150 (MQ=255)
       gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTAAACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGAt                                                                                                                                          <  2:2325085/150‑1 (MQ=255)
         agccagcTGTGCCACGCGGGTAAAGCTTTTTTCCATACGTTATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGATTTTGTACGCTGAAAACAATCAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTa                                                                                                                                        <  1:4256657/150‑1 (MQ=255)
         agccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCACGAGTACACAAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTa                                                                                                                                        <  2:3890465/150‑1 (MQ=255)
         agccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTa                                                                                                                                        <  2:3890970/150‑1 (MQ=255)
            cagcagTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGaa                                                                                                                                     >  2:1411606/1‑150 (MQ=255)
            cagcagTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGaa                                                                                                                                     >  1:2985405/1‑150 (MQ=255)
             agcagTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAAc                                                                                                                                    >  2:4291563/1‑150 (MQ=255)
              gcagTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGTGTTTGTACGCTTAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACt                                                                                                                                   <  1:3577128/150‑1 (MQ=255)
               cagTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTg                                                                                                                                  <  2:3498410/150‑1 (MQ=255)
               cagGGCCACGCGGGTAACGCTTATTTCCATCCGGGATTAACCTCATGCGACCGGTCGCACGAGTACCCCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAACTCGAGATTAGAACTg                                                                                                                                  <  2:1172641/150‑1 (MQ=255)
                 gTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTggc                                                                                                                                >  1:2372256/1‑149 (MQ=255)
                  tGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTggca                                                                                                                               >  1:1044448/1‑148 (MQ=255)
                    ccACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCaaa                                                                                                                             >  2:3119106/1‑150 (MQ=255)
                     cACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAAc                                                                                                                            >  2:1177736/1‑150 (MQ=255)
                        gcgGGTAAAGCTTATTTCCATAGGTGATTATCCTCATGCGAAGGGTCGCAAGAGTACACCAAAAACCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAg                                                                                                                         <  1:3883197/150‑1 (MQ=255)
                             tAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCa                                                                                                                    >  2:3308344/1‑150 (MQ=255)
                               aaGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCaga                                                                                                                  >  1:3758570/1‑150 (MQ=255)
                                aGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCagag                                                                                                                 >  2:3759318/1‑150 (MQ=255)
                                  cTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCt                                                                                                               >  1:4691806/1‑150 (MQ=255)
                                   ttATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTa                                                                                                              >  2:1302220/1‑150 (MQ=255)
                                    tATTTCCATACGTGATTATCCTCATGCGAACGGTCTCAAGAGTACACAAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCGGTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTAc                                                                                                             >  1:2283932/1‑150 (MQ=255)
                                    tATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTAGGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTAc                                                                                                             <  2:2425494/150‑1 (MQ=255)
                                       ttCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGAt                                                                                                          <  2:1601284/150‑1 (MQ=255)
                                       tacaaTACTTGATTATCCTCCTGCGACCGGTCGCAAGAGTACACCAAAAAACTGTTTTTGTACGCTGAAAACATTGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGAt                                                                                                          <  2:872045/146‑1 (MQ=255)
                                          cATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGtt                                                                                                       <  2:1962215/150‑1 (MQ=255)
                                          cATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGtt                                                                                                       >  1:2962674/1‑150 (MQ=255)
                                           aTACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTa                                                                                                      >  1:3514743/1‑150 (MQ=255)
                                           aTACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTa                                                                                                      >  1:1807848/1‑150 (MQ=255)
                                           aTACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTa                                                                                                      >  1:4212032/1‑150 (MQ=255)
                                           aTACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTa                                                                                                      >  1:897685/1‑150 (MQ=255)
                                           aTACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTa                                                                                                      >  1:4351872/1‑150 (MQ=255)
                                              cGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCa                                                                                                   >  2:2965514/1‑150 (MQ=255)
                                               gTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAg                                                                                                  >  2:2266866/1‑150 (MQ=255)
                                                  aTTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAgtgt                                                                                               <  2:788945/150‑1 (MQ=255)
                                                      tCCTCATGCGACCGGTCGCAAGAGAACACCAAAAAACGGAGTTTGTCCGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGAt                                                                                           <  1:2164529/150‑1 (MQ=255)
                                                      tCCTCATGCGAACGGTCGCAAGATTACACAAAAAAACTGAGTTTTTACGCTGAAAACAATCAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGAt                                                                                           <  1:1409590/150‑1 (MQ=255)
                                                      tCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGAt                                                                                           <  2:370040/150‑1 (MQ=255)
                                                       ccTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGGTGATg                                                                                          >  1:3629320/1‑150 (MQ=255)
                                                        cTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGc                                                                                         >  1:55717/1‑150 (MQ=255)
                                                        cTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGc                                                                                         >  1:4102385/1‑150 (MQ=255)
                                                        cTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGc                                                                                         >  2:2885079/1‑150 (MQ=255)
                                                           aTGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAg                                                                                      >  1:4354137/1‑150 (MQ=255)
                                                           aTGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAg                                                                                      >  2:3437333/1‑150 (MQ=255)
                                                            tGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTGCGATAGAAATCGAGATTAGACCTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGc                                                                                     >  2:1875057/1‑150 (MQ=255)
                                                              cGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGt                                                                                   >  1:4589632/1‑150 (MQ=255)
                                                               gAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTc                                                                                  >  2:3565415/1‑150 (MQ=255)
                                                                aaCGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCa                                                                                 >  2:4743388/1‑150 (MQ=255)
                                                                  cGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTTAAAAAAATCAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACg                                                                               <  1:1438565/150‑1 (MQ=255)
                                                                    gTCGCAAGAGTACCCCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAg                                                                             <  1:331950/150‑1 (MQ=255)
                                                                      cGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATACATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTg                                                                           <  2:4691806/150‑1 (MQ=255)
                                                                      cGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTg                                                                           <  1:2766615/150‑1 (MQ=255)
                                                                       gCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGa                                                                          <  2:4381326/150‑1 (MQ=255)
                                                                           gagTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTg                                                                      <  2:2297965/150‑1 (MQ=255)
                                                                            agTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTgg                                                                     >  1:1931445/1‑150 (MQ=255)
                                                                            agTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTgg                                                                     >  2:4498447/1‑150 (MQ=255)
                                                                            agTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTgg                                                                     >  2:2239059/1‑150 (MQ=255)
                                                                             gTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGt                                                                    <  2:1735879/150‑1 (MQ=255)
                                                                                 acCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtc                                                                >  1:3553866/1‑150 (MQ=255)
                                                                                 acCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtc                                                                >  1:1244221/1‑150 (MQ=255)
                                                                                   cAAAAAACAGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCATGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtcga                                                              >  2:1641559/1‑149 (MQ=255)
                                                                                     aaaaaCTTAGTTTGTACGATGAAAAAAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTcc                                                            <  2:73729/150‑1 (MQ=255)
                                                                                     aaaaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTcc                                                            <  2:2466538/150‑1 (MQ=255)
                                                                                     aaaaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTcc                                                            >  2:2069680/1‑150 (MQ=255)
                                                                                      aaaaCTTAGTTTGAACGCTTAAAAAAATCAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCca                                                           <  1:30424/150‑1 (MQ=255)
                                                                                       acaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcag                                                          <  1:3952531/148‑1 (MQ=255)
                                                                                         aCTGATTTTGTCCGCTCAAAACAATCAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcagca                                                        <  1:1210739/150‑1 (MQ=255)
                                                                                         aCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcagca                                                        >  2:470647/1‑150 (MQ=255)
                                                                                          cTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcagcag                                                       >  1:2608140/1‑150 (MQ=255)
                                                                                           tGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAgg                                                      >  1:3313641/1‑150 (MQ=255)
                                                                                             aGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGtt                                                    <  2:3015963/150‑1 (MQ=255)
                                                                                                  gTACGCTCAAAACAATTAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGAttt                                               <  1:2041538/150‑1 (MQ=255)
                                                                                                     cGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGt                                            >  1:1282377/1‑150 (MQ=255)
                                                                                                        tGAAAACAATTAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGt                                         <  1:4175318/150‑1 (MQ=255)
                                                                                                         gAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCCGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTc                                        >  1:926476/1‑150 (MQ=255)
                                                                                                           aaaCAATGAAAAATGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTGGTCGTCCAGCAGGTTGATTTTGTAGTCaa                                      >  2:1807640/1‑150 (MQ=255)
                                                                                                               aaTGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGGTGATGTTGTAGTCAACGtc                                  >  2:593486/1‑149 (MQ=255)
                                                                                                                aTGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAg                                 <  1:142045/150‑1 (MQ=255)
                                                                                                                 tGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAgg                                >  1:3366192/1‑150 (MQ=255)
                                                                                                                  gAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGt                               >  2:3812838/1‑150 (MQ=255)
                                                                                                                  gAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGt                               <  1:1116813/150‑1 (MQ=255)
                                                                                                                   aaaaaaGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTg                              >  2:1962370/1‑150 (MQ=255)
                                                                                                                   aaaaaaGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTg                              >  1:3863792/1‑150 (MQ=255)
                                                                                                                     aaaaGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGa                            >  1:4298493/1‑150 (MQ=255)
                                                                                                                       aaGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACa                          >  2:3920507/1‑150 (MQ=255)
                                                                                                                       aaGGGCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACa                          >  2:1509153/1‑150 (MQ=255)
                                                                                                                           gCCCGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGANATGtt                      <  1:187065/150‑1 (MQ=255)
                                                                                                                             ccGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGtttt                    >  2:3039967/1‑150 (MQ=255)
                                                                                                                              cGCAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGttttt                   <  1:4286410/150‑1 (MQ=255)
                                                                                                                                cAGGCCCTTTGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttgt                 >  1:4476446/1‑150 (MQ=255)
                                                                                                                                        ttGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtag         <  1:2222110/150‑1 (MQ=255)
                                                                                                                                         tGTTCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtagt        <  2:2770677/150‑1 (MQ=255)
                                                                                                                                            tCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtagg     <  1:2885079/150‑1 (MQ=255)
                                                                                                                                            tCGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtagg     <  1:2200714/150‑1 (MQ=255)
                                                                                                                                             cGATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtaggt    <  1:3119106/150‑1 (MQ=255)
                                                                                                                                              gATATAAATCGAGATTAGAACTGGCAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtaggta   >  2:2502863/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAG  >  CP009273/2304969‑2305261

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: