Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,927,972 G→T 71.5% E279* (GAA→TAA)  rbsA → fused D‑ribose transporter subunits of ABC superfamily: ATP‑binding components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,927,9720GT71.5% 333.9 / 256.2 260E279* (GAA→TAA) rbsAfused D‑ribose transporter subunits of ABC superfamily: ATP‑binding components
Reads supporting (aligned to +/- strand):  ref base G (44/30);  new base T (99/87);  total (143/117)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

Only 100 of 263 total aligned reads displayed.

CCGAAGATTCGCTGATTGAGATGATGGTGGGTCGCAAGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGCAGGATGGCCTGGCAAACGGC  >  CP009273/3927823‑3928121
                                                                                                                                                     |                                                                                                                                                     
  gAAGATTCGCTGATTGAGATGATGGTGGGTCGCAAGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTaa                                                                                                                                                     >  2:2364074/1‑150 (MQ=255)
    aGATTCGCTGATTGAGATGATGGTGGGTCGCAAGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAAt                                                                                                                                                   <  2:37910/150‑1 (MQ=255)
       ttCGCTGATTGAGATGATGGTGGGTCGCAAGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCt                                                                                                                                                <  2:2041981/150‑1 (MQ=255)
        tctcTGATTGATATGATGCTGGGTCGCAAGATGGAAGATAAATATCCGCACCATGACAAATCGCCGGGAGATATCCGCCAGAAAGTCGATAATCTCTGCGGACATCGCGTTAACGATGTCTATTTTACTATAAGCAAAGGCTAAAATCtt                                                                                                                                               >  2:968075/4‑150 (MQ=255)
        tCGCTGATTGAGATGATGGTGGGTCGCAAGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCtt                                                                                                                                               >  2:675744/1‑150 (MQ=255)
         cGCTGATTGAGATGATGGTGGGTCGCAAGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTg                                                                                                                                              >  2:1421291/1‑150 (MQ=255)
                 gagaTGATGGTGGGTCGCAAGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTctct                                                                                                                                      <  2:4195061/150‑1 (MQ=255)
                 gagaTGATGGTGGGTCGCAAGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTctct                                                                                                                                      <  2:4194999/150‑1 (MQ=255)
                   gatgatGGTGGGTCGCAAGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTgg                                                                                                                                    >  2:3476153/1‑150 (MQ=255)
                   gatgatGGTGGGTCGCAAGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTgg                                                                                                                                    >  2:3470916/1‑150 (MQ=255)
                    atgatgGTGGGTCGCAAGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGt                                                                                                                                   <  1:1197926/150‑1 (MQ=255)
                    atgatgGTGGGTCGCAAGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGt                                                                                                                                   <  1:1197752/150‑1 (MQ=255)
                          gtggGTCGCAAGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATg                                                                                                                             <  2:1778335/150‑1 (MQ=255)
                            gggTCGCAAGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATggg                                                                                                                           >  2:2086726/1‑150 (MQ=255)
                             ggTCGCAAGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGgc                                                                                                                          >  2:107819/1‑150 (MQ=255)
                              gTCGCAAGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGgcg                                                                                                                         <  1:336212/150‑1 (MQ=255)
                                 gCAAGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAg                                                                                                                      <  1:1260972/150‑1 (MQ=255)
                                  cAAGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAgg                                                                                                                     >  1:1082799/1‑150 (MQ=255)
                                   aaGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGt                                                                                                                    >  1:3983723/1‑150 (MQ=255)
                                       tGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTa                                                                                                                <  1:618046/150‑1 (MQ=255)
                                       tGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTa                                                                                                                >  1:961555/1‑150 (MQ=255)
                                        ggAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTAc                                                                                                               >  2:4093597/1‑150 (MQ=255)
                                        ggAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTAc                                                                                                               <  1:2109297/150‑1 (MQ=255)
                                           aGATCAATATCCGCACCTGGACAACGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGa                                                                                                            <  1:3756285/150‑1 (MQ=255)
                                           aGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGa                                                                                                            <  1:3756254/150‑1 (MQ=255)
                                           aGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGa                                                                                                            >  1:4244115/1‑150 (MQ=255)
                                            gATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGaa                                                                                                           >  2:3151832/1‑150 (MQ=255)
                                            gATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGaa                                                                                                           >  1:1771852/1‑150 (MQ=255)
                                               cAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACtg                                                                                                        <  2:195187/150‑1 (MQ=255)
                                               cAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACtg                                                                                                        >  2:1916770/1‑150 (MQ=255)
                                                aaTATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACtga                                                                                                       <  1:849747/150‑1 (MQ=255)
                                                 atatCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACtgat                                                                                                      <  1:3997450/150‑1 (MQ=255)
                                                           cTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGctc                                                                                            >  2:870579/1‑150 (MQ=255)
                                                           cTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGctc                                                                                            >  1:4300681/1‑150 (MQ=255)
                                                           cTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGctc                                                                                            >  2:464514/1‑150 (MQ=255)
                                                            tGGACAAAGCGCCGGGATATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGctct                                                                                           <  1:1633740/150‑1 (MQ=255)
                                                             ggACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTa                                                                                          >  2:4388827/1‑150 (MQ=255)
                                                                  aaGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGgcg                                                                                     <  1:45719/150‑1 (MQ=255)
                                                                  aaGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGgcg                                                                                     <  1:45816/150‑1 (MQ=255)
                                                                  aaGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGgcg                                                                                     >  2:788546/1‑150 (MQ=255)
                                                                     cgcCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCAc                                                                                  >  2:216135/1‑150 (MQ=255)
                                                                      gcCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTAACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACt                                                                                 <  1:359871/150‑1 (MQ=255)
                                                                          ggAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACcg                                                                             >  1:4467094/1‑150 (MQ=255)
                                                                           gagaTATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACcgc                                                                            >  2:192658/1‑150 (MQ=255)
                                                                            agaTATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACcgcg                                                                           <  1:1640728/150‑1 (MQ=255)
                                                                            agaTATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACcgcg                                                                           <  2:4057255/150‑1 (MQ=255)
                                                                                       tGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGtt                                                                <  1:2245928/150‑1 (MQ=255)
                                                                                        gAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTa                                                               <  1:3208642/150‑1 (MQ=255)
                                                                                        gAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTa                                                               >  2:1813879/1‑150 (MQ=255)
                                                                                           aGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGt                                                            >  1:2033647/1‑150 (MQ=255)
                                                                                            gTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGAGGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTc                                                           <  1:253074/150‑1 (MQ=255)
                                                                                             tCGATAATCTCTGCGGACCTGGAGTTAACGATGTCTCTTTTACTTTACGGAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCa                                                          <  1:3604293/150‑1 (MQ=255)
                                                                                               gATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCAcc                                                        >  2:880701/1‑150 (MQ=255)
                                                                                                 tAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCt                                                      <  1:3897382/150‑1 (MQ=255)
                                                                                                 tAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCt                                                      <  1:3897390/150‑1 (MQ=255)
                                                                                                  aaTCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTg                                                     >  1:3377650/1‑150 (MQ=255)
                                                                                                  aaTCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTg                                                     <  2:110794/150‑1 (MQ=255)
                                                                                                  aaTCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTg                                                     <  1:378685/150‑1 (MQ=255)
                                                                                                  aaTCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTg                                                     >  1:4342859/1‑150 (MQ=255)
                                                                                                     ctctGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGAt                                                  >  1:2803417/1‑150 (MQ=255)
                                                                                                     ctctGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGAt                                                  >  1:2804661/1‑150 (MQ=255)
                                                                                                       ctGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATgg                                                <  1:1716257/150‑1 (MQ=255)
                                                                                                        tgtggACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATggg                                               >  2:3526874/4‑150 (MQ=255)
                                                                                                         gCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGc                                              <  1:564364/150‑1 (MQ=255)
                                                                                                            gACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATg                                           <  1:2494205/150‑1 (MQ=255)
                                                                                                            gACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATg                                           <  1:1398495/150‑1 (MQ=255)
                                                                                                             aCCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGa                                          >  1:1561759/1‑150 (MQ=255)
                                                                                                             aCCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGa                                          >  1:1561805/1‑150 (MQ=255)
                                                                                                             aCCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGa                                          >  1:1561800/1‑150 (MQ=255)
                                                                                                              ccTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGaa                                         >  2:2396732/1‑150 (MQ=255)
                                                                                                               cTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAg                                        <  2:3624542/150‑1 (MQ=255)
                                                                                                                tGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGt                                       >  2:2414453/1‑150 (MQ=255)
                                                                                                                tGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGt                                       >  1:2504117/1‑150 (MQ=255)
                                                                                                                    gTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGtt                                   >  2:1708645/1‑150 (MQ=255)
                                                                                                                    gTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGtt                                   <  2:2740121/150‑1 (MQ=255)
                                                                                                                        aCGATGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTAccc                               >  2:2985400/1‑150 (MQ=255)
                                                                                                                           aTGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGtt                            <  1:312663/150‑1 (MQ=255)
                                                                                                                            tGTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTc                           >  2:1906388/1‑150 (MQ=255)
                                                                                                                             gTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCa                          <  1:1847106/150‑1 (MQ=255)
                                                                                                                             gTCTCTTTTACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCa                          <  1:4388827/150‑1 (MQ=255)
                                                                                                                             gTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCa                          <  1:94302/150‑1 (MQ=255)
                                                                                                                             gTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCa                          <  1:94972/150‑1 (MQ=255)
                                                                                                                                 ctTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGc                      >  1:1154120/1‑150 (MQ=255)
                                                                                                                                  ttttACTTTACGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGCa                     <  2:3909719/150‑1 (MQ=255)
                                                                                                                                  ttttACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGCa                     >  2:3608612/1‑150 (MQ=255)
                                                                                                                                   tttactttacGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGCAg                    >  2:3994945/1‑150 (MQ=255)
                                                                                                                                    ttactttacGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGCAgg                   >  1:3788150/1‑150 (MQ=255)
                                                                                                                                    ttactttacGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGCAgg                   <  1:4043212/150‑1 (MQ=255)
                                                                                                                                    ttactttacGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGCAgg                   <  1:1789727/150‑1 (MQ=255)
                                                                                                                                       ctttacGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGCAGGATg                <  1:1001372/150‑1 (MQ=255)
                                                                                                                                          tacGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGCAGGATGGcc             <  1:3425691/150‑1 (MQ=255)
                                                                                                                                           acGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGCAGGATGGCCt            >  1:4117354/1‑150 (MQ=255)
                                                                                                                                           acGCAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGCAGGATGGCCt            >  1:3796626/1‑150 (MQ=255)
                                                                                                                                           acGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGCAGGATGGCCt            >  2:222026/1‑150 (MQ=255)
                                                                                                                                            cGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGCAGGATGGCCTg           >  2:3475501/1‑150 (MQ=255)
                                                                                                                                              cAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGCAGGATGGCCTGGc         <  2:3030679/150‑1 (MQ=255)
                                                                                                                                              cAAAGGCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGCAGGATGGCCTGGc         >  1:320929/1‑150 (MQ=255)
                                                                                                                                                  ggCTAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGCAGGATGGCCTGGCAAAc     <  2:1359975/150‑1 (MQ=255)
                                                                                                                                                     gAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGCAGGATGGCCTGGCAAACGGc  >  2:1429533/1‑150 (MQ=255)
                                                                                                                                                     gAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGCAGGATGGCCTGGCAAACGGc  >  2:3248274/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGAAGATTCGCTGATTGAGATGATGGTGGGTCGCAAGCTGGAAGATCAATATCCGCACCTGGACAAAGCGCCGGGAGATATCCGCCTGAAAGTCGATAATCTCTGCGGACCTGGCGTTAACGATGTCTCTTTTACTTTACGCAAAGGCGAAATTCTTGGCGTCTCTGGTTTGATGGGCGCAGGTCGTACCGAACTGATGAAAGTGCTCTACGGCGCACTACCGCGCACCAGCGGTTACGTCACCCTGGATGGGCATGAAGTCGTTACCCGTTCACCGCAGGATGGCCTGGCAAACGGC  >  CP009273/3927823‑3928121

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: