Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 815,168 C→A 47.3% R140R (CGG→AGG)  moaE → molybdopterin synthase, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273815,1680CA47.3% 150.9 / 308.7 234R140R (CGG→AGG) moaEmolybdopterin synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base C (50/73);  new base A (59/52);  total (109/125)
Fisher's exact test for biased strand distribution p-value = 6.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 236 total aligned reads displayed.

CGATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTC  >  CP009273/815023‑815317
                                                                                                                                                 |                                                                                                                                                     
cGATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGga                                                                                                                                                   <  2:395387/150‑1 (MQ=255)
 gATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGgag                                                                                                                                                  >  2:1051828/1‑150 (MQ=255)
 gATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGgag                                                                                                                                                  >  2:3327066/1‑150 (MQ=255)
    gAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGc                                                                                                                                               <  1:3149231/150‑1 (MQ=255)
         cGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATca                                                                                                                                          >  1:1008013/1‑150 (MQ=255)
         cGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATca                                                                                                                                          >  1:255185/1‑150 (MQ=255)
            ttttGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATcagca                                                                                                                                       <  1:1116573/150‑1 (MQ=255)
            ttttGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATcagca                                                                                                                                       <  1:639683/150‑1 (MQ=255)
                 tCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAggcgg                                                                                                                                  >  2:2278004/1‑150 (MQ=255)
                 tCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAggcgg                                                                                                                                  >  2:4375115/1‑150 (MQ=255)
                    gtgtCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCaa                                                                                                                               <  2:100700/150‑1 (MQ=255)
                     tgtCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCaaa                                                                                                                              >  2:3706369/1‑150 (MQ=255)
                      gtCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCACAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGATGAAGCTCGGGAGAGCGATCAGCAGGCGGCaaaa                                                                                                                             <  1:4423205/150‑1 (MQ=255)
                       tCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAAc                                                                                                                            <  1:2894442/150‑1 (MQ=255)
                            aGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTgg                                                                                                                       <  2:1761865/150‑1 (MQ=255)
                            aGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTgg                                                                                                                       >  2:1349281/1‑150 (MQ=255)
                                 gcATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGtt                                                                                                                  >  2:4244027/1‑150 (MQ=255)
                                  cATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGttt                                                                                                                 <  1:728365/150‑1 (MQ=255)
                                       cagcagTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGt                                                                                                            <  2:3196754/150‑1 (MQ=255)
                                        agcagTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGtt                                                                                                           >  2:1046619/1‑150 (MQ=255)
                                         gcagTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCCAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTa                                                                                                          <  1:2151426/150‑1 (MQ=255)
                                         gcagTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTa                                                                                                          <  1:581347/150‑1 (MQ=255)
                                          cagTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAg                                                                                                         <  1:2185377/150‑1 (MQ=255)
                                            gTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGcc                                                                                                       >  2:1304222/1‑150 (MQ=255)
                                              gagTTTGAAGCCGGGCAGTTTATTATGTATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCgg                                                                                                     >  1:566214/3‑150 (MQ=255)
                                              gCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCgg                                                                                                     <  1:1190241/150‑1 (MQ=255)
                                                gTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGAt                                                                                                   <  1:3860803/150‑1 (MQ=255)
                                                  ttGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATaa                                                                                                 <  1:3391355/150‑1 (MQ=255)
                                                   tGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAg                                                                                                <  2:1142362/150‑1 (MQ=255)
                                                   tGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAg                                                                                                <  2:903112/150‑1 (MQ=255)
                                                      aGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCa                                                                                             >  2:392616/1‑150 (MQ=255)
                                                      aGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCa                                                                                             <  2:2674868/150‑1 (MQ=255)
                                                      aGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCa                                                                                             <  2:2675145/150‑1 (MQ=255)
                                                        ccGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCAcc                                                                                           <  1:183131/150‑1 (MQ=255)
                                                             cAGTTTATTATGGATTATCTCCAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGc                                                                                      <  1:3043702/150‑1 (MQ=255)
                                                             cAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGc                                                                                      >  2:1867545/1‑150 (MQ=255)
                                                              aGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCg                                                                                     >  1:3968886/1‑150 (MQ=255)
                                                                 ttattaTGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTcc                                                                                  <  2:1544134/150‑1 (MQ=255)
                                                                 ttattaTGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTcc                                                                                  >  1:3259325/1‑150 (MQ=255)
                                                                  tattatGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCg                                                                                 >  1:1557097/1‑150 (MQ=255)
                                                                   attatGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCgg                                                                                <  2:2352931/150‑1 (MQ=255)
                                                                   attatGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCgg                                                                                >  2:3802763/1‑150 (MQ=255)
                                                                    ttatGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGc                                                                               >  1:233988/1‑150 (MQ=255)
                                                                     tatGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCa                                                                              <  2:461913/150‑1 (MQ=255)
                                                                       tGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAATCGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAAt                                                                            >  2:300787/1‑150 (MQ=255)
                                                                       tGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAAt                                                                            >  1:3944057/1‑150 (MQ=255)
                                                                        ggATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAAtt                                                                           <  2:338266/150‑1 (MQ=255)
                                                                        ggATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAAtt                                                                           >  1:4321692/1‑150 (MQ=255)
                                                                         gATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAAttt                                                                          >  1:1187779/1‑150 (MQ=255)
                                                                         gATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAAttt                                                                          >  1:2029899/1‑150 (MQ=255)
                                                                           ttATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCa                                                                        >  1:2995083/1‑150 (MQ=255)
                                                                            tATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGGGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCAc                                                                       >  1:4437216/1‑150 (MQ=255)
                                                                            tATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCAc                                                                       <  1:3705682/150‑1 (MQ=255)
                                                                             aTCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCAcc                                                                      <  2:1553247/150‑1 (MQ=255)
                                                                                tcAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt                                                                   >  1:3253609/1‑150 (MQ=255)
                                                                                     aCCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTa                                                              <  2:1499763/150‑1 (MQ=255)
                                                                                     aCCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTa                                                              >  1:3265081/1‑150 (MQ=255)
                                                                                     aCCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTa                                                              >  1:988668/1‑150 (MQ=255)
                                                                                      cccGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTAc                                                             <  2:3571778/150‑1 (MQ=255)
                                                                                       ccGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACa                                                            >  2:568850/1‑150 (MQ=255)
                                                                                        cgcgCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACAt                                                           >  2:428367/1‑150 (MQ=255)
                                                                                          cgcACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATAc                                                         <  1:1671614/150‑1 (MQ=255)
                                                                                          cgcACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATAc                                                         <  1:491737/150‑1 (MQ=255)
                                                                                             aCCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGtt                                                      >  1:3711473/1‑150 (MQ=255)
                                                                                              ccGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGttt                                                     <  2:2413450/150‑1 (MQ=255)
                                                                                                gTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTAc                                                   <  1:1236957/150‑1 (MQ=255)
                                                                                                  tCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGc                                                 <  1:3384753/150‑1 (MQ=255)
                                                                                                    tGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTg                                               <  2:3718722/150‑1 (MQ=255)
                                                                                                    tGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGATTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTg                                               <  2:3793486/150‑1 (MQ=255)
                                                                                                    tGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTg                                               >  1:2119789/1‑150 (MQ=255)
                                                                                                     ggAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGc                                              <  1:1192812/150‑1 (MQ=255)
                                                                                                       aaGCGCGAAGCCACGCCGGAAGGCGACCTCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCACAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGt                                            >  1:4046088/1‑150 (MQ=255)
                                                                                                       aaGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGt                                            <  1:1423401/150‑1 (MQ=255)
                                                                                                       aaGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGt                                            <  1:2468143/150‑1 (MQ=255)
                                                                                                         gcgcGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCg                                          <  1:3765639/150‑1 (MQ=255)
                                                                                                            cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATg                                       <  1:3967526/150‑1 (MQ=255)
                                                                                                             gAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGc                                      <  2:3695427/150‑1 (MQ=255)
                                                                                                             gAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGc                                      <  2:2782252/150‑1 (MQ=255)
                                                                                                              aaGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGca                                     >  2:676096/1‑150 (MQ=255)
                                                                                                                  cACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGc                                 >  1:116914/1‑150 (MQ=255)
                                                                                                                    cGCCGGAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCt                               <  1:1505707/150‑1 (MQ=255)
                                                                                                                        ggAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCtcatc                           >  2:2093233/1‑150 (MQ=255)
                                                                                                                         gAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCtcatca                          <  1:1420740/150‑1 (MQ=255)
                                                                                                                         gAAGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCtcatca                          <  2:1487855/150‑1 (MQ=255)
                                                                                                                           aGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCAtt                        >  1:4492331/1‑150 (MQ=255)
                                                                                                                           aGGCGACCGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCAtt                        <  2:3550473/150‑1 (MQ=255)
                                                                                                                            ggCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCAttt                       <  2:3968590/150‑1 (MQ=255)
                                                                                                                            ggCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCAttt                       <  2:3370558/150‑1 (MQ=255)
                                                                                                                                 ccGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAg                  <  1:1052025/150‑1 (MQ=255)
                                                                                                                                 ccGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAg                  <  1:1051745/150‑1 (MQ=255)
                                                                                                                                 ccGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAg                  <  1:3327066/150‑1 (MQ=255)
                                                                                                                                  cGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGt                 <  2:2866811/150‑1 (MQ=255)
                                                                                                                                  cGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGt                 >  1:3974630/1‑150 (MQ=255)
                                                                                                                                  cGCTGGGTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGATTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGt                 <  2:2803006/150‑1 (MQ=255)
                                                                                                                                     tGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATc              <  2:566214/150‑1 (MQ=255)
                                                                                                                                      gggTTGAAGCTAGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATcc             >  1:2009298/1‑150 (MQ=255)
                                                                                                                                         ttGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTa          >  2:1261490/1‑150 (MQ=255)
                                                                                                                                             aGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGat      <  2:3875021/150‑1 (MQ=255)
                                                                                                                                               cTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGatat    >  2:3391042/1‑150 (MQ=255)
                                                                                                                                                 cGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTc  >  1:3031193/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CGATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTC  >  CP009273/815023‑815317

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: