Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,378,778 G→T 100% S137* (TCG→TAG)  mscM ← mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,378,7780GT98.4% 1027.0 / 3.8 315S137* (TCG→TAG) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base G (5/0);  new base T (140/170);  total (145/170)
Fisher's exact test for biased strand distribution p-value = 1.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 315 total aligned reads displayed.

GCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCATGTTG  >  CP009273/4378632‑4378924
                                                                                                                                                  |                                                                                                                                                  
gcgcCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCTGTTTGCTGTTGCGGCAGTTGATTCAGCTAAt                                                                                                                                                 <  1:4123418/150‑1 (MQ=255)
gcgcCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAAt                                                                                                                                                 <  1:3910207/150‑1 (MQ=255)
 cgcCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATc                                                                                                                                                <  1:294943/150‑1 (MQ=255)
  gcCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCg                                                                                                                                               <  1:2760000/150‑1 (MQ=255)
   ccTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCgg                                                                                                                                              <  2:2505763/150‑1 (MQ=255)
       aaGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAAt                                                                                                                                          >  2:3202664/1‑150 (MQ=255)
        aGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAAtc                                                                                                                                         <  2:115355/150‑1 (MQ=255)
        aGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAAtc                                                                                                                                         >  2:1205729/1‑150 (MQ=255)
           cGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTcg                                                                                                                                      >  1:2612670/1‑150 (MQ=255)
           cGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTcg                                                                                                                                      <  2:3829293/150‑1 (MQ=255)
             tGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTcgcg                                                                                                                                    <  1:4929807/150‑1 (MQ=255)
                agagTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGgc                                                                                                                                 >  1:4976075/1‑150 (MQ=255)
                 gagTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGgcg                                                                                                                                <  2:1646263/150‑1 (MQ=255)
                      aGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTc                                                                                                                           >  1:745943/1‑150 (MQ=255)
                      aGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTc                                                                                                                           >  2:4579613/1‑150 (MQ=255)
                       gACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTcc                                                                                                                          >  1:1721206/1‑150 (MQ=255)
                       gACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTcc                                                                                                                          >  1:4298562/1‑150 (MQ=255)
                                    aaTTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTggg                                                                                                             >  1:4596655/1‑150 (MQ=255)
                                    aaTTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGGTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTggg                                                                                                             <  2:3530388/150‑1 (MQ=255)
                                       tttGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCt                                                                                                          >  2:1557968/1‑150 (MQ=255)
                                        ttGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTg                                                                                                         <  2:3049070/150‑1 (MQ=255)
                                        ttGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTg                                                                                                         <  2:2273541/150‑1 (MQ=255)
                                            gCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACgg                                                                                                     <  1:3671285/150‑1 (MQ=255)
                                            gCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACgg                                                                                                     <  1:1222785/150‑1 (MQ=255)
                                             ccTGAGTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGc                                                                                                    <  2:1602722/150‑1 (MQ=255)
                                               tGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCtt                                                                                                  >  1:1018871/1‑150 (MQ=255)
                                                 aTTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCtttt                                                                                                <  1:1343845/150‑1 (MQ=255)
                                                  ttGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTa                                                                                               <  2:1987205/150‑1 (MQ=255)
                                                  ttGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTa                                                                                               <  1:2408699/150‑1 (MQ=255)
                                                   tGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTAt                                                                                              >  1:2518914/1‑150 (MQ=255)
                                                    gAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATc                                                                                             <  2:3191770/150‑1 (MQ=255)
                                                     aGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATcc                                                                                            >  2:5071122/1‑150 (MQ=255)
                                                      gCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTCACTGACGGCGGGCGTCGGTTTGCTGGTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCa                                                                                           <  2:3012799/150‑1 (MQ=255)
                                                      gCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCa                                                                                           <  1:4212620/150‑1 (MQ=255)
                                                      gCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCa                                                                                           <  2:750667/150‑1 (MQ=255)
                                                       cGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAg                                                                                          >  2:4253633/1‑150 (MQ=255)
                                                        ggAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGc                                                                                         <  2:2421104/150‑1 (MQ=255)
                                                        ggAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGc                                                                                         <  1:1078805/150‑1 (MQ=255)
                                                        ggAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGc                                                                                         >  1:389503/1‑150 (MQ=255)
                                                         gAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCCTTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCa                                                                                        <  1:183855/150‑1 (MQ=255)
                                                          aGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCaa                                                                                       >  1:4357691/1‑150 (MQ=255)
                                                          aGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCaa                                                                                       >  1:4454343/1‑150 (MQ=255)
                                                          aGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCaa                                                                                       <  1:1605298/150‑1 (MQ=255)
                                                            tATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACt                                                                                     >  2:3439614/1‑150 (MQ=255)
                                                             aTTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTg                                                                                    <  2:4937270/150‑1 (MQ=255)
                                                                gCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGgct                                                                                 >  1:1780749/1‑150 (MQ=255)
                                                                gCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGgct                                                                                 >  2:2279843/1‑150 (MQ=255)
                                                                 ccGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGTCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGgctg                                                                                <  1:3089502/150‑1 (MQ=255)
                                                                  cGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGgctgc                                                                               <  2:73172/150‑1 (MQ=255)
                                                                   ggTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGgctgct                                                                              >  2:3949582/1‑150 (MQ=255)
                                                                   ggTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGgctgct                                                                              <  2:2076496/150‑1 (MQ=255)
                                                                    gTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGctgctg                                                                             <  1:789667/150‑1 (MQ=255)
                                                                    gTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGctgctg                                                                             >  1:2988918/1‑150 (MQ=255)
                                                                      aaGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGAc                                                                           >  2:4855245/1‑150 (MQ=255)
                                                                        gCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCt                                                                         <  2:1020277/150‑1 (MQ=255)
                                                                         cGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTg                                                                        <  1:2702166/150‑1 (MQ=255)
                                                                          gTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTgg                                                                       >  2:279286/1‑150 (MQ=255)
                                                                            tCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGgag                                                                     <  1:1440252/150‑1 (MQ=255)
                                                                             cccAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGgaga                                                                    <  2:1805522/150‑1 (MQ=255)
                                                                              ccAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGaa                                                                   >  1:2890946/1‑150 (MQ=255)
                                                                                 ggcggcGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAAttt                                                                <  2:3626391/150‑1 (MQ=255)
                                                                                   tggcgCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTcc                                                              >  2:247750/2‑150 (MQ=255)
                                                                                   cggcgCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTcc                                                              <  1:2217815/150‑1 (MQ=255)
                                                                                   cggcgCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTcc                                                              <  2:2550548/150‑1 (MQ=255)
                                                                                    ggcgCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCt                                                             <  2:3266663/150‑1 (MQ=255)
                                                                                    ggcgCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCt                                                             <  1:2319464/150‑1 (MQ=255)
                                                                                    ggcgCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTGTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCt                                                             >  2:670500/1‑150 (MQ=255)
                                                                                    ggcgCTCGATCTCATTTAACTGACGGCGGGCGTCGGTGTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCt                                                             <  1:3469073/150‑1 (MQ=255)
                                                                                        cTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGAtt                                                         <  1:4823313/150‑1 (MQ=255)
                                                                                          cGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCa                                                       >  2:3560504/1‑150 (MQ=255)
                                                                                           gATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAg                                                      <  2:4155438/150‑1 (MQ=255)
                                                                                            aTCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAgc                                                     <  2:2860989/150‑1 (MQ=255)
                                                                                             tctcATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAgcg                                                    >  1:2208548/1‑150 (MQ=255)
                                                                                             tctcATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAgcg                                                    >  1:1980082/1‑150 (MQ=255)
                                                                                                cATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGt                                                 <  2:1390854/150‑1 (MQ=255)
                                                                                                 aTTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTc                                                >  1:2636917/1‑150 (MQ=255)
                                                                                                 aTTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTc                                                >  1:1120520/1‑150 (MQ=255)
                                                                                                  tttAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCg                                               <  2:2684432/150‑1 (MQ=255)
                                                                                                  tttAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCg                                               <  2:4308452/150‑1 (MQ=255)
                                                                                                  tttAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCg                                               <  2:4976271/150‑1 (MQ=255)
                                                                                                     aaCTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTa                                            <  1:399567/150‑1 (MQ=255)
                                                                                                     aaCTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTa                                            >  1:3943244/1‑150 (MQ=255)
                                                                                                        tGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGAc                                         >  2:819325/1‑150 (MQ=255)
                                                                                                           cggcggGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACAtt                                      >  2:596781/1‑150 (MQ=255)
                                                                                                               gggcgTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCg                                  >  2:5145901/1‑150 (MQ=255)
                                                                                                               gggcgTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCg                                  >  2:2663133/1‑150 (MQ=255)
                                                                                                                ggcgTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCgg                                 >  1:3182058/1‑150 (MQ=255)
                                                                                                                 gcgTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGc                                <  1:924148/150‑1 (MQ=255)
                                                                                                                    tCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGac                             <  2:2968429/150‑1 (MQ=255)
                                                                                                                      ggTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGacac                           >  2:5034189/1‑150 (MQ=255)
                                                                                                                       gTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACg                          >  1:4696295/1‑150 (MQ=255)
                                                                                                                                 gCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCt                <  2:3351712/150‑1 (MQ=255)
                                                                                                                                 gCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCt                <  1:4328738/150‑1 (MQ=255)
                                                                                                                                  cGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTc               >  2:4620993/1‑150 (MQ=255)
                                                                                                                                   ggCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCg              >  1:530949/1‑150 (MQ=255)
                                                                                                                                      aGTTGATTCAGCTAATCGGCAATCTCGCTGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCa           >  2:345146/1‑150 (MQ=255)
                                                                                                                                        ttGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACg         <  2:3771160/150‑1 (MQ=255)
                                                                                                                                          gATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCa       <  1:2293468/150‑1 (MQ=255)
                                                                                                                                           aTTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCAt      >  2:3811524/1‑150 (MQ=255)
                                                                                                                                            ttCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCAtg     <  1:14642/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCATGTTG  >  CP009273/4378632‑4378924

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: