Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,594,244 T→G 12.8% L90V (TTA→GTA)  bglJ → bgl operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,594,2440TG12.8% 529.6 / 18.3 228L90V (TTA→GTA) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (83/109);  new base G (34/2);  total (117/111)
Fisher's exact test for biased strand distribution p-value = 2.21e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 235 total aligned reads displayed.

CCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGAC  >  CP009273/4594095‑4594393
                                                                                                                                                     |                                                                                                                                                     
 cTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGttt                                                                                                                                                      <  2:1938310/150‑1 (MQ=255)
  tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTa                                                                                                                                                     <  1:1529183/150‑1 (MQ=255)
  tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTa                                                                                                                                                     >  2:3118748/1‑150 (MQ=255)
     aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGtt                                                                                                                                                  >  2:3774340/1‑150 (MQ=255)
     aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGtt                                                                                                                                                  >  2:2109642/1‑150 (MQ=255)
          tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGc                                                                                                                                             >  2:280132/1‑150 (MQ=255)
           cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg                                                                                                                                            <  2:3782652/150‑1 (MQ=255)
           cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg                                                                                                                                            <  2:444792/150‑1 (MQ=255)
           cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg                                                                                                                                            <  1:4331576/150‑1 (MQ=255)
           cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTATTTTTCTGCCATGAGAAGGGAGCGCAGAGCGGGCTTATATTTCCTGACTGTACAGGCGATTAAGTTTCCGCGTACCCGGCGGTGTGTTAtttgg                                                                                                                                            >  1:2209126/1‑147 (MQ=255)
             ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTGTAGTTATTGCGga                                                                                                                                          >  2:3912995/1‑150 (MQ=255)
               aGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCCTGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatg                                                                                                                                        >  2:1488622/1‑150 (MQ=255)
               aGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatg                                                                                                                                        >  2:1116945/1‑150 (MQ=255)
               aGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatg                                                                                                                                        >  2:2285517/1‑150 (MQ=255)
               aGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGGTATTTCGgatg                                                                                                                                        >  2:3313531/1‑150 (MQ=255)
                  cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatg                                                                                                                                     <  1:2301644/150‑1 (MQ=255)
                  cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatg                                                                                                                                     <  1:2699449/150‑1 (MQ=255)
                   aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatga                                                                                                                                    >  2:151940/1‑150 (MQ=255)
                    tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatgat                                                                                                                                   <  1:1400347/150‑1 (MQ=255)
                      tCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGatat                                                                                                                                 <  1:2968954/150‑1 (MQ=255)
                       cAGGATGCCATGTCGCGGATCTCGTTTGCGGCTGTCATTTTTTCTTTTTCTGACATGAGAAGTTAGCGCATAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATAtt                                                                                                                                <  2:1944028/150‑1 (MQ=255)
                        aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTg                                                                                                                               >  2:2779846/1‑150 (MQ=255)
                        aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTg                                                                                                                               >  2:235072/1‑150 (MQ=255)
                           aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                                                                                                            <  1:1058675/150‑1 (MQ=255)
                             gCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCt                                                                                                                          <  1:264514/150‑1 (MQ=255)
                             gCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCATTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCt                                                                                                                          <  1:631313/150‑1 (MQ=255)
                                aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCgg                                                                                                                       <  1:3956829/150‑1 (MQ=255)
                                  gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCt                                                                                                                     <  1:546232/150‑1 (MQ=255)
                                     gcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGGTATTGCGGATGATGATATTGAAGCTCGGCTGAt                                                                                                                  >  2:1256016/1‑150 (MQ=255)
                                     gcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTATGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAATTTTCCGCGTACCCGGCGTGTGGTTTTTTGGGGTGGTGTTATTTAATCTCgggggtt                                                                                                                  >  1:3749602/1‑145 (MQ=255)
                                      cgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGAtt                                                                                                                 <  1:3774340/150‑1 (MQ=255)
                                       ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTg                                                                                                                <  2:2158562/150‑1 (MQ=255)
                                       ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTg                                                                                                                <  1:3052281/150‑1 (MQ=255)
                                           ctcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTTCGGATGATGATATTGAATCAAGGCTGATTTGtta                                                                                                            >  2:2594701/1‑149 (MQ=255)
                                                 tGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTc                                                                                                      <  1:3868871/150‑1 (MQ=255)
                                                 tGCGGCGGTCATTTTTTCTTTTTCAGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTc                                                                                                      <  1:1650448/150‑1 (MQ=255)
                                                    ggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGcc                                                                                                   >  2:2123242/1‑150 (MQ=255)
                                                      cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt                                                                                                 <  2:2436575/150‑1 (MQ=255)
                                                       ggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTATTTATTTCGGATGTTGATATTTAAGCTCGGCGGATTGGTTAATTTTCGCCaac                                                                                                >  1:287919/1‑148 (MQ=255)
                                                       ggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc                                                                                                <  2:4330341/150‑1 (MQ=255)
                                                       ggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGGCGCCATc                                                                                                >  2:1276197/1‑150 (MQ=255)
                                                         tCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTAGGCTGATTGGGTCATTTTCGCCATCAc                                                                                              >  2:4068794/1‑150 (MQ=255)
                                                          cATTTTGTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTGTAGTTCTTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAcc                                                                                             <  1:1897870/150‑1 (MQ=255)
                                                            ttttttCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAATCTCGGCTGTTTGGTTAATTGTCGCCATAACCGc                                                                                           >  2:3859274/1‑150 (MQ=255)
                                                             tttttctttttcTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATTATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt                                                                                          >  2:1777287/1‑150 (MQ=255)
                                                              ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTg                                                                                         <  2:84488/150‑1 (MQ=255)
                                                                ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTGTAGTTATTGCGGATGTTTATATTGAATCTCGGCTGATTGGTTCATTTTCGCCATCACCGCGGGa                                                                                       >  1:1666774/1‑150 (MQ=255)
                                                                 tctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc                                                                                      <  1:4101919/150‑1 (MQ=255)
                                                                    ttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGgt                                                                                   <  1:2853951/150‑1 (MQ=255)
                                                                       tctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCTCCATCACCGCTGGACGGTGTa                                                                                >  2:2017271/1‑150 (MQ=255)
                                                                       tctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTa                                                                                <  1:1329355/150‑1 (MQ=255)
                                                                        ctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAGTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTAt                                                                               <  2:4010560/150‑1 (MQ=255)
                                                                          gCCATGAGACGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGTTTTAGTTTCCGCGTACCCGGCGTGTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTa                                                                             <  2:24522/150‑1 (MQ=255)
                                                                          gCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTa                                                                             <  1:3118748/150‑1 (MQ=255)
                                                                                 gaAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAg                                                                      <  2:1224464/150‑1 (MQ=255)
                                                                                 gaAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAg                                                                      >  2:2351690/1‑150 (MQ=255)
                                                                                  aaGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTTCGGATGATGATATTTAAGCTCGGCTGATTTGTTCATTGTAGCCAACACCCCGGGACGGTGTATTAAGTaaaac                                                                     >  1:2445376/1‑148 (MQ=255)
                                                                                          gcAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTTCGGATGATGATATTTAAGCTCGGCTGATTTGTTGATTGTCGCCAACACCGCAGGACCGTGTATTAACTAAAACGTCAACGc                                                             >  2:2695517/1‑150 (MQ=255)
                                                                                            agagagGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTTTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTg                                                           >  2:1040570/1‑150 (MQ=255)
                                                                                            agagagGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTGTTGTTTTTTGGGGTGATGATATTTTAGCTCGGGTGGTTGGTTAATTTTCGCCCACACCCCGGGACGGGGTATTAAGTAAATAGTCaaacctg                                                           >  1:1950951/1‑145 (MQ=255)
                                                                                              agagGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTGTAGTTATTGCGGACGATGATATTGAAGCTCGGCTGATTGGGTCATTGTCGCCATCACCGCTGGACGGGGTATTAAGTAAAGCGTCAACGCTggt                                                         >  2:187763/1‑149 (MQ=255)
                                                                                               gagGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGGAGATGATATTGAATCTCGGCTGATTGGTTCATTTTCGCCATCACCGCTGGTCGGTGTATTAAGTAAAGCGTAAACGCTGgag                                                        >  2:871345/1‑150 (MQ=255)
                                                                                               gagGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGgag                                                        >  2:3977150/1‑150 (MQ=255)
                                                                                                agGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTGTAGTTATTGCGGATGTTGATATTGAAGCTCGGCGGTTTTGTTCAATGTCGCCATCAACGCTGGGCGGGGTATTTAGTAAAGCGTCAACGCgggtgt                                                       >  1:3248900/1‑147 (MQ=255)
                                                                                                 gggATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTTTCGCCATCACCGCTGGACGGTGTATGAAGTTAAGCGGCAACGCTGGAGAt                                                      >  2:4148900/1‑150 (MQ=255)
                                                                                                   gATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGAttt                                                    <  2:4167696/150‑1 (MQ=255)
                                                                                                        tCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTtcat                                               <  1:1304053/150‑1 (MQ=255)
                                                                                                           tGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTTTTTCGGGTTATGTTAATTTATCTCGGCGGGTTTGTTCAATGTGGCCAACACCGCGGGGCGGGGTAATTATTAAATCGTCAACGCTGGAGATTTTTGATAAg                                            >  2:1456068/1‑150 (MQ=255)
                                                                                                           tGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAg                                            <  1:3520823/150‑1 (MQ=255)
                                                                                                              ctgactgaACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTTCGGATGTTTATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAAAACCGCTGGACGGTGTTTTTAGTAAAGCGGACACGCTGGTGATTTTTGATAAggga                                         >  2:54391/1‑148 (MQ=255)
                                                                                                              ctgactgaACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGaa                                         <  1:483095/150‑1 (MQ=255)
                                                                                                              ctgactgaACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGaa                                         <  1:4163697/150‑1 (MQ=255)
                                                                                                                gactgaACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAAct                                       <  2:1684826/150‑1 (MQ=255)
                                                                                                                  ctgaACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAActct                                     <  2:3709869/150‑1 (MQ=255)
                                                                                                                   tgaACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGTTGATGATATTGAAGCTCGGCTGATTGGTTCCTTGTAGCCAACACCGCGGGACGGTGTATTTAGTAAAGCGTCAACGCTGGAGATTTTTCATAAGGAACTCtt                                    >  2:2641128/1‑150 (MQ=255)
                                                                                                                    gaACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCAGTGTCGCCCACACCGCTGGGCGGTGGCTTAAGGAAAGCGTCAACGCGGGAGATTTTTAATCAGGAACTCttt                                   >  2:3394358/1‑150 (MQ=255)
                                                                                                                      aCTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCttttt                                 <  2:2565687/150‑1 (MQ=255)
                                                                                                                      aCTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCttttt                                 <  2:3035296/150‑1 (MQ=255)
                                                                                                                        tGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTTCGGATGATGATATTGAAGCTCGGCTGATTGGTTCAGTTTCGCCATCACCGCGGGTCGGTGGATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATAAGGAACTCttttttt                               >  2:1836913/1‑148 (MQ=255)
                                                                                                                         ggCGATTAAGTTTCCGCGTACCCGGCGTTTTGTTATTGCGGATGATGATATTGAAGCGCGGCTGGTTGGTTCCTTGTCGCCATCACCGCGGGACGGTGTATTAAGTAAAGCGTCAACGCGGGAGATTTTTCATCAGGAACACTGTTTGTc                              >  2:2230383/1‑150 (MQ=255)
                                                                                                                          gCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCa                             <  1:1203621/150‑1 (MQ=255)
                                                                                                                          gCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCa                             <  1:228095/150‑1 (MQ=255)
                                                                                                                            gATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTTCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGGATTAAGTAAAGCGTCAAAGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCAtt                           >  2:1665982/1‑150 (MQ=255)
                                                                                                                            gATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCAtt                           <  2:368410/150‑1 (MQ=255)
                                                                                                                            gATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCAtt                           <  1:1116945/150‑1 (MQ=255)
                                                                                                                            gATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCAtt                           <  2:4281573/150‑1 (MQ=255)
                                                                                                                            gATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCAGGAGATTTTTCATCAGGAACACTTTTTGTCAtt                           >  2:3586179/1‑150 (MQ=255)
                                                                                                                            gATTAAGTTTCCGCGTACCCGGCGTTGAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCAtt                           <  1:3313531/150‑1 (MQ=255)
                                                                                                                             aTTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTa                          <  1:72500/150‑1 (MQ=255)
                                                                                                                             aTTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCGGGAGATTTTTCATCAGGAAATCTTTTTGTCATTa                          >  2:4222241/1‑150 (MQ=255)
                                                                                                                              ttAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTaa                         <  2:2445272/150‑1 (MQ=255)
                                                                                                                                aaGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAAt                       >  2:3567648/1‑150 (MQ=255)
                                                                                                                                  gTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATgg                     <  2:287919/150‑1 (MQ=255)
                                                                                                                                        gcgTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACg               <  2:1347879/150‑1 (MQ=255)
                                                                                                                                        gcgTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACg               <  1:4129693/150‑1 (MQ=255)
                                                                                                                                        gcgTACCCGGCGTTTAGTTATTGCGGATGATGAGATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACg               <  2:3844416/150‑1 (MQ=255)
                                                                                                                                           tACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCa            <  2:3937371/150‑1 (MQ=255)
                                                                                                                                           tACCCGGCGTGTAGTTATTTCGGAAGATTATAATTAATCTCGGCTGATTTGTTAATTTTCGCCCAAACCCCAGGCCGGTGTATTAAGTAAAGCGTCAAACCGGGAGATTTTTAATAAAGAACACATTTTGTAAATAAATTGTGTACGtaa            >  2:3362026/1‑148 (MQ=255)
                                                                                                                                                ggCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGcga       >  1:3658800/1‑150 (MQ=255)
                                                                                                                                                   gTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGcgaccg    >  1:1747820/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGAC  >  CP009273/4594095‑4594393

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