Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,788,981 A→T 6.3% S143R (AGT→AGA btuC ← vitamin B12 transporter subunit: membrane component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,788,9810AT 6.3% 1051.8 / 44.7 314S143R (AGT→AGAbtuCvitamin B12 transporter subunit: membrane component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base A (157/137);  new base T (10/10);  total (167/147)
Fisher's exact test for biased strand distribution p-value = 8.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 320 total aligned reads displayed.

AGCCAGCTTTGCCGCCAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGCCTGCGCCGTTAGA  >  CP009273/1788832‑1789130
                                                                                                                                                     |                                                                                                                                                     
agccagcTTTGCCGCCAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGa                                                                                                                                                       <  1:1713491/150‑1 (MQ=255)
  ccagcTTTGCCGCCAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACt                                                                                                                                                     >  1:4977698/1‑150 (MQ=255)
  ccagcTTTGCCGCCAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACt                                                                                                                                                     <  2:4576460/150‑1 (MQ=255)
   cagcTTTGCCGCCAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTg                                                                                                                                                    <  2:1167259/150‑1 (MQ=255)
       tttGCCGCCAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCg                                                                                                                                                <  2:4058244/150‑1 (MQ=255)
         tGCCGCCAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGaa                                                                                                                                              <  2:3047816/150‑1 (MQ=255)
          gccgccAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGaaa                                                                                                                                             <  2:3321722/150‑1 (MQ=255)
             gccAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGAt                                                                                                                                          >  2:2477653/1‑150 (MQ=255)
              ccAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATg                                                                                                                                         >  1:1127678/1‑150 (MQ=255)
                aGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTCGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAAGAACCGACTGGTCGCAAGATgac                                                                                                                                       >  2:15450/1‑150 (MQ=255)
                aGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATgac                                                                                                                                       >  1:3530029/1‑150 (MQ=255)
                  tCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATgacga                                                                                                                                     <  2:440210/150‑1 (MQ=255)
                   cTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATgacgac                                                                                                                                    <  1:1081870/150‑1 (MQ=255)
                    tACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGTCTGGTCGAAAGATGacgacg                                                                                                                                   <  1:682202/150‑1 (MQ=255)
                     aCGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACgg                                                                                                                                  <  2:1915292/150‑1 (MQ=255)
                       gccgccAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGc                                                                                                                                >  1:1818305/1‑150 (MQ=255)
                         cgccAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGa                                                                                                                              >  1:1337703/1‑150 (MQ=255)
                         cgccAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGa                                                                                                                              <  1:471284/150‑1 (MQ=255)
                           ccAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGaaa                                                                                                                            >  1:3776727/1‑150 (MQ=255)
                           ccAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGaaa                                                                                                                            >  2:2630226/1‑150 (MQ=255)
                              aaaCCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGt                                                                                                                         <  1:2899929/150‑1 (MQ=255)
                              aaaCCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGt                                                                                                                         <  1:2900023/150‑1 (MQ=255)
                                 ccGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAag                                                                                                                      >  2:14221/1‑150 (MQ=255)
                                    cccATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGt                                                                                                                   >  2:3480205/1‑150 (MQ=255)
                                      catcatCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTAt                                                                                                                 <  2:4200038/150‑1 (MQ=255)
                                      catcatCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTAt                                                                                                                 >  1:2675472/1‑150 (MQ=255)
                                       atcatcCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTAtt                                                                                                                >  2:1300343/1‑150 (MQ=255)
                                         catcCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGTCTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTaa                                                                                                              >  1:4195914/1‑150 (MQ=255)
                                            ccAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGt                                                                                                           >  2:2273348/1‑150 (MQ=255)
                                            ccAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGt                                                                                                           <  1:2558728/150‑1 (MQ=255)
                                            ccAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGt                                                                                                           <  1:1164579/150‑1 (MQ=255)
                                             cAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCCGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTCAAGtg                                                                                                          >  2:4719674/1‑150 (MQ=255)
                                              aGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGtga                                                                                                         >  1:2198050/1‑150 (MQ=255)
                                               gTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGtgat                                                                                                        >  2:923468/1‑150 (MQ=255)
                                                tACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGtgatg                                                                                                       <  2:4826491/150‑1 (MQ=255)
                                                 aCATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGtgatga                                                                                                      <  2:4469241/150‑1 (MQ=255)
                                                  cATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGTCTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTgatgat                                                                                                     <  1:3222714/150‑1 (MQ=255)
                                                  cATCAGCTGACGCAAATCAACGGAGGTGGAAAAGTAGAGAGCCCACGTCCTTAGTGCGCTCCAGATACTCCCTAATGCAACGCCAGCCCGCACTAACCGACTGGTCGAAAGATGACGACGGGCGAAACGGCCGAGTCGGAAAGTgatgat                                                                                                     >  2:1640105/1‑150 (MQ=255)
                                                      aGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAgc                                                                                                 <  1:3050521/150‑1 (MQ=255)
                                                      aGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAgc                                                                                                 <  1:1932288/150‑1 (MQ=255)
                                                        cTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAgcgc                                                                                               <  2:5034492/150‑1 (MQ=255)
                                                             gCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAg                                                                                          >  2:3372640/1‑150 (MQ=255)
                                                              cAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGAAGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGc                                                                                         >  2:86136/1‑150 (MQ=255)
                                                               aaaTCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCa                                                                                        >  1:2519187/1‑150 (MQ=255)
                                                               aaaTCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCa                                                                                        >  1:5048883/1‑150 (MQ=255)
                                                                aaTCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCaa                                                                                       >  1:303700/1‑150 (MQ=255)
                                                                 aTCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAAt                                                                                      <  1:3107642/150‑1 (MQ=255)
                                                                   cAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCg                                                                                    <  1:843230/150‑1 (MQ=255)
                                                                   cAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCg                                                                                    >  2:719835/1‑150 (MQ=255)
                                                                    aaCTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGc                                                                                   >  2:2058395/1‑150 (MQ=255)
                                                                        gAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGcacac                                                                               >  2:799284/1‑150 (MQ=255)
                                                                          ggTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAg                                                                             <  1:1595096/150‑1 (MQ=255)
                                                                            tGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGcc                                                                           >  1:3177893/1‑150 (MQ=255)
                                                                             ggAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGccc                                                                          >  2:4090686/1‑150 (MQ=255)
                                                                              gAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCt                                                                         >  2:3806042/1‑150 (MQ=255)
                                                                               aaaaGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTa                                                                        >  1:2452433/1‑150 (MQ=255)
                                                                                aaaGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCtcg                                                                       >  2:3945820/1‑148 (MQ=255)
                                                                                aaaGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAg                                                                       <  2:2387159/150‑1 (MQ=255)
                                                                                  aGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAgcg                                                                     >  1:2298650/1‑150 (MQ=255)
                                                                                  aGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAgcg                                                                     >  1:4810203/1‑150 (MQ=255)
                                                                                   gTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAgcgc                                                                    >  1:2222476/1‑150 (MQ=255)
                                                                                    tAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGcc                                                                   >  2:2692980/1‑150 (MQ=255)
                                                                                       aTAGCCCACGTCACTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAg                                                                >  2:2145637/1‑150 (MQ=255)
                                                                                         aGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGtt                                                              <  1:1418924/150‑1 (MQ=255)
                                                                                         aGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGtt                                                              <  1:4076802/150‑1 (MQ=255)
                                                                                         aGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGtt                                                              >  1:4248288/1‑150 (MQ=255)
                                                                                           cccACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTACGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTgg                                                            >  2:4121489/1‑150 (MQ=255)
                                                                                            ccACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTggg                                                           <  1:3718698/150‑1 (MQ=255)
                                                                                              aCGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGa                                                         >  1:1015884/1‑150 (MQ=255)
                                                                                              aCGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGa                                                         >  1:4286866/1‑150 (MQ=255)
                                                                                               cGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAg                                                        >  1:3955293/1‑150 (MQ=255)
                                                                                                  cATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGttg                                                     <  2:1565915/150‑1 (MQ=255)
                                                                                                  cATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGttg                                                     <  2:1565791/150‑1 (MQ=255)
                                                                                                   aTTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGttgc                                                    >  2:3635860/1‑150 (MQ=255)
                                                                                                     tAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGttgccc                                                  >  2:4595442/1‑150 (MQ=255)
                                                                                                      aGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGttgccct                                                 >  1:2256684/1‑150 (MQ=255)
                                                                                                             tACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCaa                                          >  2:582924/1‑150 (MQ=255)
                                                                                                             tACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCaa                                          <  2:2917639/150‑1 (MQ=255)
                                                                                                                  aTAATCCCTAATGCAACGCCAGCCAGCAATAACCGTCTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAAt                                     >  2:1409573/1‑150 (MQ=255)
                                                                                                                  aTAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAAt                                     <  2:37424/150‑1 (MQ=255)
                                                                                                                  aTAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAAt                                     >  2:1808448/1‑150 (MQ=255)
                                                                                                                     aTCCCTAATGCAACGCCAGCCAGCAATAACCGTCTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATAcc                                  <  1:2080366/150‑1 (MQ=255)
                                                                                                                     aTCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATAcc                                  >  1:346204/1‑150 (MQ=255)
                                                                                                                        ccTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACcgcg                               <  1:954814/150‑1 (MQ=255)
                                                                                                                        ccTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACcgcg                               <  2:4846625/150‑1 (MQ=255)
                                                                                                                         cTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCgcgg                              >  2:3157165/1‑150 (MQ=255)
                                                                                                                          tAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCgcggc                             <  2:3404496/150‑1 (MQ=255)
                                                                                                                              gCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATa                         <  2:1331532/150‑1 (MQ=255)
                                                                                                                                aaCGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAg                       >  1:2650214/1‑150 (MQ=255)
                                                                                                                                 tcGCCAGCCAGCTAGAACCGACTGGTCGATAGATGACGTCGGGCGAATCGTAAGAGTATTACAGTGATGATAAGCGCGCCAGCTATCGCACACAGCCCTAGCGCCCTGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGc                      <  2:2662046/149‑1 (MQ=255)
                                                                                                                                   gccagccaGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGccc                    <  1:896369/150‑1 (MQ=255)
                                                                                                                                    ccagccagCAATAACCGTCTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGcccc                   <  1:364113/150‑1 (MQ=255)
                                                                                                                                    ccagccagCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGcccc                   >  1:4769306/1‑150 (MQ=255)
                                                                                                                                     cagccagcAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCa                  <  2:1137664/150‑1 (MQ=255)
                                                                                                                                     cagccagcAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCa                  <  1:3324534/150‑1 (MQ=255)
                                                                                                                                         cagcAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGcc              <  1:550620/150‑1 (MQ=255)
                                                                                                                                            cAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGCCTgc           <  2:4972858/150‑1 (MQ=255)
                                                                                                                                                 aCCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGCCTGCGCCGt      <  1:3999434/150‑1 (MQ=255)
                                                                                                                                                   cGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGGTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGCCTGCGCCGTTa    <  2:2899171/150‑1 (MQ=255)
                                                                                                                                                    gACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGCCTGCGCCGTTag   <  1:4457696/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCCAGCTTTGCCGCCAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGCCTGCGCCGTTAGA  >  CP009273/1788832‑1789130

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: